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INTRODUCTION 


Vascularization  plays  a  key  role  in  the  growth  and  metastasis  of  solid  tumors  [1-6],  In  two 
recent  clinical  studies,  breast  cancer  patients  whose  tumors  had  a  high  vascular  density  subsequently 
went  on  to  develop  metastases  over  a  follow  up  period  of  2.5  years  [7,  8],  Statistical  analyses  of 
these  patients  showed  that  vascular  density  was  the  single  most  important  factor  (p<0.006) 
associated  with  subsequent  formation  of  metastasis  [8];  the  other  factors  examined  were  epidermal 
growth  factor  receptor  status  (p<0.01),  node  status  (p<0.02),  estrogen  receptor  status  (p<  0.05), 
tumor  size  (p<0.06),  tumor  grade  (p<0.5),  c-erb-2  expression  (p<0.7),  p53  (p<0.8)  and  tumor  type 
(p<0.8).  Studies  correlating  vascularization  with  metastasis  have  so  far  been  performed  with 
histological  evaluation  of  excised  tissue  specimens  [7,  8]  as  a  result  of  which  information  regarding 
functioning  of  vessels  is  lost.  Similarly,  the  physiological  environment  of  these  tumors,  in  terms  of 
acidity  and  lactate  production  remains  unknown.  Thus  a  lack  of  noninvasive  methods  has  left  some 
vital  questions  about  the  precise  nature  of  the  relationship  between  vascularization  and  metastasis 
unanswered.  Tumor  neovascularization  is  induced  by  the  secretion  of  angiogenic  factors  which  act 
as  chemotactic  factors  and  mitogens  for  endothelial  cells  [1, 4-6].  One  of  the  most  potent  of  these  is 
vascular  endothelial  growth  factor  (VEGF).  VEGF  also  increases  vascular  permeability  which  in 
turn  may  allow  cancer  cells  greater  access  to  the  vasculature  [9].  In  glioblastoma  multiformae  areas 
of  necrosis  and  hypoxia  show  a  higher  expression  of  VEGF  [10, 11].  Poorly  functioning  vessels  and 
the  associated  hypoxia  and  necrosis  may  play  a  role  in  attracting  further  vascularization.  Areas  of 
hypoxia  are  also  associated  with  accumulation  of  lactate  and  low  pH.  These  two  physiological 
factors  attract  neovascularization  by  stimulating  the  secretion  of  angiogenic  factors  from 
macrophages  [12-15].  The  secretion  of  enzymes  which  degrade  the  basement  membrane  in  the 
metastatic  process  increases  at  low  pH  [16,  17].  Thus,  vascularization,  the  physiological 
environment,  and  formation  of  metastases  are  highly  interdependent.  An  understanding  of  the  role 
of  the  physiological  environment  in  vascularization  and  metastasis,  and  the  dependence  of  this 
environment  on  metastatic  phenotype  are  essential  to  delineate  the  relationship  between 
vascularization  and  metastasis.  Questions  which  are  central  to  understanding  this  relationship  are  - 
(1)  does  the  metastatic  phenotype  induce  a  higher  degree  of  vascularization  and  is  this  mediated  by 
modulation  of  the  physiological  environment  ?  (addressed  in  Specific  Aims  1  and  3)  (2)  If  so,  do 
nonmetastatic  tumors  and  preneoplastic  tissue  exhibit  proportionately  lower  vascularization  ? 
(addressed  in  Specific  Aims  1  and  2)  (3)  Which  particular  property  of  the  vascularization  e.g. 
permeability,  or  vascular  volume  is  the  dominant  factor  in  the  dependence  of  metastasis  on 
vascularization  ?  (addressed  in  Specific  Aims  1  and  2)  (4)  Is  a  significant  fraction  of  the  vessels 
observed  in  the  histological  studies  non-functional  and  does  the  resultant  unsuitable  environment 
induce  expression  of  signals  or  substances  which  prompt  and  enable  the  cells  to  metastasize  ? 
(addressed  in  Specific  Aim  1).  The  overall  goal  of  this  research  proposal  is  to  use  noninvasive 
Magnetic  Resonance  (MR)  Imaging  (I)  and  Spectroscopy  (S)  to  answer  the  questions  posed  above. 
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*  The  research  proposed  consists  of  three  closely  related  aims  designed  to  unravel  the  complex 
relationship  between  vascularization  and  metastasis.  Our  overall  goal  in  this  project  is  to  determine 
key  vascular  and  physiological  properties  which  result  in  the  close  relationship  between  vascular 
density  and  metastasis  as  this  information  may  ultimately  be  used  to  prevent  tumor  metastasis.  We 
had  proposed  the  following  three  aims: 

Aim  1:  To  investigate  the  relationship  between  metastatic  phenotype  and  the  development  of 
vascularization  and  evaluate  the  functionality  of  the  developing  vascularization  in  terms  of  vascular 
volume,  vascular  permeability  and  relative  perfusion. 

Hypothesis  #1:  More  metastatic  lines  will  exhibit  a  higher  level  of  vascularization  and  permeability 
for  similar  volumes.  A  significant  number  of  vessels  detected  by  immunoperoxidase  staining  will 
not  be  functional  and  this  number  will  increase  with  the  size  of  the  tumor. 

Aim  2:  To  investigate  the  effect  of  increasing  (a)  tumor  vascularization  and  (b)  tumor 
vascularization  and  permeability  on  the  formation  of  metastases. 

Hypothesis  #2:  Higher  vascularization  and  permeability  will  lead  to  an  increase  or  an  earlier 
incidence  of  metastases  for  all  the  lines. 

(Aims  1  and  2  are  related  to  questions  1-4  outlined  in  background) 

Aim  3:  To  determine  the  relationship  between  metastatic  phenotype  and  intra-  and  extracellular  pH 
and  lactate  production. 

Hypothesis  #3:  More  metastatic  lines  will  be  more  glycolytically  active  in  vivo,  creating  an 
environment  of  high  lactate  and  low  extracellular  pH  for  volume  matched  lesions. 

(Aim  3  is  related  to  question  1  outlined  in  background) 

In  the  progress  report  for  Year  1  we  had  presented  preliminary  data  which  demonstrated  that  there 
were  significant  differences  in  the  vascular  volume  generated  by  a  invasive  metastatic  human  breast 
cancer  line  MDA-MB-231  and  a  nonmetastatic  animal  cancer  line  RIF-1.  Studies  correlating  VEGF 
distribution  with  MRI  maps  of  vascular  volume  and  permeability  demonstrated  that  areas  around 
necrosis  showed  high  expression  of  VEGF  and  were  more  permeable.  In  Year  1  we  also  presented 
data  to  demonstrate  that  there  were  significant  differences  in  pH  regulation  and  the  phospholipid 
metabolism  in  solid  tumors  growing  in  vivo  in  SCID  mice  for  a  highly  metastatic  and  a  less 
metastatic  human  breast  cancer  line.  We  had  therefore  made  significant  progress  in  Aim  1  and  Aim 
3  by  the  end  of  the  first  year  of  the  grant. 

In  Year  2  we  had  the  following  achievements.  First  we  showed  that  the  alterations  in  phospholipid 
metabolism  observed  for  a  highly  metastatic  and  less  metastatic  tumor  could  be  generalized  to  an 
entire  panel  of  human  mammary  epithelial  cells,  ranging  from  normal  to  highly  malignant  and 
metastatic.  We  also  demonstrated  that  lactate  levels  significantly  increased  in  malignant  cells 
compared  to  normal  and  immortalized  hyperplastic  human  mammary  epithelial  cells.  In  addition  we 
developed  an  implemented  multi-slice  imaging  pulse  sequences  instead  of  the  originally  proposed 
single  slice  studies  of  tumors  to  determine  the  vascular  characteristics  over  the  entire  tumor.  In 
addition  we  developed  a  visualization  software  program  which  would  allow  us  to  interactively  relate 
MRI  information  with  histological  maps. 
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In  Year  3  we  have  used  the  significant  technical  advances  made  in  year  2  to  comprehensively 
*  characterize  vascular  patterns  for  three  human  breast  cancer  lines,  MDA-MB-435,  MDA-MB-23 1 
and  MCF-7  inoculated  in  the  mammary  fat  pad  of  SCID  mice. 

The  technical  objectives  outlined  in  our  statement  of  work  continue  to  remain  the  same  and  these 
are:  Delineate  the  role  of  vascular  volume,  permeability  and  perfusion  and  tumor  physiological 
environment  in  the  formation  of  metastasis  jfrom  human  breast  cancer  lines  with  preselected 
invasive  and  metastatic  potential  grown  in  SCID  (severe  combined  immune  deficient)  mice. 

RXPRRTMRNTAL  METHODS  USED  FOR  YEAR  3 

Imaging  studies  were  performed  on  a  GE  Omega  4.7T  instrument  widi  a  solenoidal  coil 
placed  around  the  tumor  (volumes  200-300mm3)  for  human  breast  cancer  models  MCF-7,  MDA- 
MB-23 1  and  MDA-MB-435  inoculated  in  the  mammary  fat  pad  or  flank  of  SCID  mice.  The  tail 
vein  of  the  anesthetized  ^imal  was  c^theterized  before  it  vyas  placed  in  the  magnet.  Aoimal  body 
temperature  was  maintained  at  37  ^C  by  heat  generated  from  a  warm  water  pad.  Multi-slice 
relaxation  rates  (Tr^)  in  tumors  were  obtained  by  a  saturation  recovery  method  combined  with 
fast  Ti  SNAPSHOT-FLASH  imaging  (flip  angle  of  7^,  echo  time  of  2ms).  Images  of  8  slices 
(slice  thickness  of  1mm)  acquired  with  an  in-plane  spatial  resolution  of  125  pm  (128  x  128 
matrix,  16mm  field  of  view,  NS=8)  were  obtained  for  3  relaxation  delays  (lOOms,  SOOms,  and  Is) 
for  each  of  the  slices.  Thus,  5-8  T  i  maps  (depending  upon  the  size  of  the  tumor)  were  acquired 
within  7  minutes.  An  Mq  map  with  a  recovery  delay  of  7s  was  acquired  once  at  the  beginning  of 
the  experiment.  Images  were  obtained  before  i.v.  administration  of  0.2ml  of  60  mg/ml  albumin- 
GdDlBA  in  saline  (dose  of  500mg/kg)  and  repeated  every  8  minutes,  starting  10  minutes  after  the 
injection,  up  to  32  minutes.  Relaxatfon  maps  were  r^onst^cted  from  dma  sets  for  three  different 
relaxation  times  arid  the  Mq  data  set.  At  the  end  of  the  imaging  studies,  the  animal  was  sacrificed, 
0.5  ml  of  blood  was  withdrawn  from  the  inferior  vena  cava,  and  tumors  were  marked  for 
referencing  to  the  MRI  images,  excised,  and  fixed  in  10%  formalin  for  sectioning  and  staining. 

Another  technical  advance  we  made  this  year  was  to  implement  a  coil  which  allowed  us  to 
determine  blood  Tl  concentrations  from  a  couple  of  drops  of  blood  obtained  from  the  tail  vein 
(shown  in  Figure  1).  The  advantage  of  this  development  is  that  it  is  not  necessary  to  sacrifice  the 
animal  to  determine  blood  Tl's  and  therefore  it  is  possible  to  obtain  repeated  measurements  of 
vascular  volume  and  permeability  over  a  period  of  time.  This  would  be  very  useful  for 
determining  the  effects  of  therapeutic  interventions  of  vascular  characteristics  in  the  same  tumor 
before  and  after  treatment. 

Vascular  volume  and  permeability  product  surface  area  (PS)  maps  were  generated  from 
the  ratio  of  A(1/T1)  values  in  the  images  to  that  of  blood.  The  slope  of  A(1/T1)  ratios  versus  time 
in  each  pixel  was  used  to  compute  (PS)  while  the  intercept  of  the  line  at  zero  time  was  used  to 
compute  vascular  volume.  Thus,  vascular  volumes  were  corrected  for  permeability  of  the 
vessels.  Ten  pairs  of  adjacent  5pm  thick  histological  sections  were  obtained  at  500  pm  intervals 


7 


through  the  tumor  were  stained  with  hematoxylin  eosin  or  immunostained  for  VEGF  expression. 
Sections  were  digitized  with  a  Sanyo  CCD  camera  attached  to  an  optical  microscope. 

Blood  Tj  measurements  can  be  obtained  at  4.7T  on  a  18  pi  sample 
using  drops  of  blood  from  the  tail  vein  without  sacrificing  animal 


Coil  designed  by  V.  P.  Chacko 

Tj  from  18  pi  sample  =119  msec 

Tj  from  0.3  ml  sample  taken  five  minutes  after  from  same  animal  =  124  msec 
Figure  1:  Coil  design  used  to  determine  blood  Tl's  from  a  couple  of  drops  of  blood. 

3-D  reconstruction  of  both  MRI  and  histological  sections  was  performed  using  the 
Clinical  Microscope  Visualization  software  that  we  are  developing  as  a  research  tool  for  use  in 
our  investigation.  The  software  consists  of  two  parts,  the  first  part  forms  the  volume  image  from 
2-D  images  and  the  second  part  performs  visualization  of  the  volumetric  model  in  suitable  3-D 
perspective.  The  high  resolution  image  'tiles'  representing  the  histological  sections  acquired  from 
the  microscope  are  stitched  along  the  overlapping  borders  to  produce  a  large  'virtual'  field  of  view 
image  that  represents  the  full  cross  section  which  could  not  be  acquired  otherwise  in  a  single 
actual  field  of  view  at  the  selected  high  resolution.  Such  a  virtual  field  of  view  image  can  be 
panned  around,  zoomed  in,  to  look  at  selected  regions  at  the  full  resolution  in  which  they  are 
acquired.  Each  crossection  is  then  layered  one  over  the  other,  to  compose  the  3-D  volume.  The 
volume  image  is  then  interactively  rendered  using  the  volumetric  visualization  software,  to  adjust 
the  transfer  functions  that  control  the  voxel  transparency  and  intensity  characteristics  of  various 
structures  of  interest,  to  delineate  it  from  the  surrounding  structures.  The  software  system  is 
developed  around  Silicon  Graphics  Inc.,  Workstation  systems,  taking  full  advantage  of  the 
hardware  accelerated  graphics  capabilities  such  as  2-  and  3-  D  Textures  to  provide  interactive 
rendering  results. 

RESULTS 

Since  we  are  evaluating  the  role  of  vascularization  in  metastasis  it  was  critical  to  establish  that  our 
human  breast  cancer  models  were  indeed  metastasizing.  Shown  in  Figure  2  are  examples  of  lung 
sections  obtained  from  mice  bearing  MDA-MB-435  and  MDA-MB-231  breast  tumors  in  the 
mammary  fat  pad,  demonstrating  the  presence  of  metastatic  nodules.  Metastasis  from  MDA-MB- 
231  was  less  frequent  than  from  MDA-MB-435  supporting  observations  made  by  other 
investigators  [18]. 
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Figure  2:  High  power  photomicrographs  demonstrating  metastatic  nodules  in 
5  (xm  thick  H  &  E  stained  lung  sections  obtained  from  mice  bearing  (a)  an 
MDA-MB-435  tumor  and  (b)  an  MDA-MB-235  tumor.  Metastatic  nodules 
are  marked  by  arrows. 

The  MRJ  method  used  to  generate  maps  of  vascular  volume  and  permeability  is  based  on  the 
assumption  that  the  concentration  of  the  contrast  agent  remains  constant  within  the  vascular  system 
over  the  period  of  observation.  Although  this  has  been  demonstrated  in  larger  animals  such  as  rats 
and  rabbits  [19]  this  has  not  been  demonstrated  in  mice.  We  therefore  canulated  the  mouse  carotid 
and  used  our  miniature  coil  to  estimate  Tl's  in  blood  in  the  same  animal  over  a  period  of  time.  The 
data  obtained  for  three  different  animals  are  presented  in  Figure  3  and  demonstrate  that  the 
concentration  of  albumin-GdDTPA  remains  constant  in  the  vascular  system  upto  35  minutes.  Our 
MRI  studies  are  completed  within  32  minutes. 
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Figure  3:  Blood  relaxation  rates  following  administration  of  500mg/kg  albumin-GdpTPA. 
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Figure  4:  Representative  multi-slice  data  obtained  from  a  human  breast  cancer  model. 


Typical  multi-slice  data  of  vascular  volume  and  permeability  together  with  histological  sections 
stained  with  hematoxylin  and  eosin  and  VEGF  expression,  obtained  from  an  MDA-MB-435-ip 
tumor  are  shown  in  Figure  4.  The  3D  reconstructed  maps  obtained  from  these  data  are  shown  in 
Figure  5. 


Anotlier  major  effort  of  this  yeai'  was  to  develop  a  software  program  for  volmnetric  analysis  of 
regions  of  interest  within  the  tumor  such  as  regions  of  high  vascular  volume,  or  permeability  and  the 
relationship  between  the  two  for  the  different  tumor  models.  An  example  of  this  software  analysis  is 
shown  in  Figure  6. 


3D  aod  tri  planaF  V  ievvs  of  fused  map 


Figure  5:  3D  reconstructed  maps  of  vascular  volume  and  permeability.  The  fused  map  is 
obtained  by  displaying  vascular  volume  through  a  greeen  channel  and  vascular  permeability 
through  the  red  channel.  The  absence  of  yellow  in  the  fused  images  demonstrate  that  there  is 
hardly  any  overlap  of  regions  of  high  vascular  volume  with  high  permeability. 


MO  map 


to  navigate 


Permeability  map 
to  navigate 
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T^igure  6:  Example  of  software  analysis. 


Data  obtained  from  MCF-7,  MDA-MB-23 1  and  MDA-NfB-435  tumors  are  summarized  below 
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Vascular  volume  profile 


MCF-7  MDA-MB-231  MDA-MB-435 


Tumor  type 

Figure  7:  Analysis  of  different  vascular  categories  within  the  human  breast  cancer 
models.  Approximately  five  to  eight  slices  wdre  obtained  from  each  tumor  and  five  to 
eight  frimors  were  studied  for  each  group.  **  represents  significantly  different  frorn 
MCF-7  tumors.  Statistical  analysis  was  performed  using  Analysis  of  Variance 
(ANOVA)  with  a  multi-comparison  significance  level  of  95%.  Tumor  volume?  were 
size  matched  and  were  of  the  order  of  200-300  mm  3. 


Vascular  permeability  profile 


CONCLTTSION S:  The  results  obtained  demonstrate  that  there  are  significant  differences  in 
vascular  characteristics  between  the  metastatic  and  non  or  less  metastatic  lines.  High  vascular 
volume  and  high  vascular  permeability  characterized  the  most  metastatic  breast  cancer  line.  The  non 
metastatic  human  breast  cancer  line  MCF-7  exhibited  a  lower  vascular  volume  as  well  as  lower 
vascular  permeability.  We  are  currently  in  the  process  of  determining  the  relationship  between 
vascular  volume  and  permeability  and  the  number  and  size  of  metastatic  lung  nodules  in  individual 
animals.  Also  apparent  from  the  data  are  the  observations  that  it  is  essential  to  characterize  the 
vasculature  of  the  regions  of  high  vascular  volume  or  permeability,  since  the  significance  disappears 
or  is  less  apparent  when  averaging  the  measurements  over  the  entire  tumor.  This  demonstrates  the 
necessity  of  techniques  which  can  detect  vascular  characteristics  with  spatial  information.  In 
addition,  this  study  fiirther  confirmed  our  previous  observations  that  regions  of  high  vascular  volmne 
and  high  vascular  permeability  do  not  coincide  spatially.  Regions  of  low  vascular  volume  were 
usually  associated  with  foci  of  necrosis  in  the  histological  sections.  High  permeability  was  related  to 
a  higher  expression  of  VEGF.  We  are  now  in  the  process  of  performing  studies  for  Aim  2  which  is 
to  investigate  the  effect  of  increasing  (a)  tumor  vascularization  and  (b)  tumor  vascularization  and 
permeability  on  the  formation  of  metastases. 


KEY  RESEARCH  ACCOMPLISHMENTS; 


The  major  findings  to  emerge  fi’om  the  research  studies  performed  thus  far  are : 

•  NMR  spectra  of  primary  tumors  in  SCID  mice  revealed  a  dramatic  and  consistent  difference  in  the 
phospholipid  composition  of  control  and  transgene  tumors  formed  by  derivatives  of  MDA-MB-435 
human  breast  carcinoma  cells  transfected  with  nm23  constructs.  Significant  differences  in  intra  and 
extracellular  pH  were  also  detected  for  solid  tumors  derived  from  Aese  lines.  This  was  one  of  the 
first  in  vivo  observations  to  link  the  activity  of  a  putative  metastasis  suppressor  gene  to  metabolic 
processes.  The  data  also  demonstrate  the  potential  of  noninvasive  NMR  spectroscopy  to  detect 
forms  of  gene  therapy  which  may  involve  transfection  of  cells  with  nm23. 

•  Choline  phospholipid  metabolite  levels  progressively  increased  in  cultured  HMEC  as  cells  become 
more  malignant.  We  therefore  propose  that  carcinogenesis  in  human  breast  epithelial  cells  results  in 
progressive  alteration  of  membrane  choline  phospholipid  metabolism.  This  work  is  relevant  to 
diagnosis  of  breast  cancer  and  also  provides  a  rationale  for  selective  pharmacological  intervention. 

•  Lactate  levels  increase  significantly  in  cultured  HMEC  following  malignant  transformation. 
However,  following  malignant  transformation,  there  did  not  appear  to  be  a  close  dependence 
between  lactate  levels  observed  in  malignant  cell  lines  and  the  metastatic  potential  of  these  lines.  The 
increased  lactate  production  may  result  in  an  acidic  environment  which  may  promote  invasive 
behavior  and  contribute  to  metastasis. 

•  3-dimensional  interactive  analysis  of  vascular  volume  and  permeability  and  histological 
morphology  demonstrates  that  areas  of  low  vascular  volume  are  associated  with  cell  death  and 
increasingly  permeable  vasculature.  Regions  of  high  vascular  volume  and  high  vascular 
permeability  do  not  coincide  spatially. 

•  The  more  metastatic  cell  lines  are  characterized  by  higher  vascular  volume  and  vascular 
permeability  in  vivo  for  analysis  performed  for  regions  of  high  vascular  volume  and  permeability. 
The  results  indicate  a  potential  use  of  MRI  for  evaluating  'metastatic  risk'  noninvasively. 
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ABSTRACT 

Transduction  of  mitogenic  signals  in  cells  can  be  mediated  by  molecules 
derived  from  the  synthesis  and  breakdown  of  the  major  membrane  phos¬ 
pholipid,  phosphotidylcholine.  Studies  were  performed  on  human  mam¬ 
mary  epithelial  cells  in  culture  to  understand  the  impact  of  malignant 
transformation  and  progression  on  membrane  phospholipid  metabolism. 
In  the  model  system  used  here,  phosphocholine  levels  and  total  choline- 
containing  phospholipid  metabolite  levels  increased  with  progression  from 
normal  to  immortalized  to  oncogene-transformed  to  tumor-derived  cells. 
These  changes  occurred  independently  of  cell  doubling  time.  A  “glycero- 
phosphocholine  to  phosphocholine  switch”  was  apparent  with  immortal¬ 
ization.  This  alteration  in  phenotype  of  increased  phosphocholine  relative 
to  glycerophosphocholine  was  observed  in  oncogene-transformed  and  for 
all  human  breast  tumor  cell  lines  analyzed.  The  results  demonstrate  that 
progression  of  human  mammary  epithelial  cells  from  normal  to  malignant 
phenotype  is  associated  with  altered  membrane  choline  phospholipid 
metabolism. 

INTRODUCTION 

PtC^  is  the  most  abundant  phospholipid  in  biological  membranes 
and  together  with  other  phospholipids,  such  as  phosphatidylethano- 
lamine  and  neutral  lipids,  form  the  characteristic  bilayer  structure  of 
cells  and  regulate  membrane  integrity  (1,  2).  MCPM  (Fig.  1),  i.e., 
biosynthesis  and  hydrolysis  of  PtC,  are  essential  processes  for  mito¬ 
genic  signal  transduction  events  in  cells  (3-6).  There  is  now  increas¬ 
ing  evidence  to  suggest  that  products  of  MCPM  such  as  PCho, 
diacylglycerol,  and  arachidonic  acid  metabolites  may  function  as 
second  messengers  essential  for  the  mitogenic  activity  of  growth 
factors  particularly  in  the  activation  of  the  ras-raf-l-MAPK  cascade  and 
protein  kinase  C  pathway  (3-6).  Together  with  inositol  phospholipid 
metabolism,  MCPM  can  also  provide  a  sustained  activation  of  mitogenic 
signal  transduction  via  a  positive  feedback  interaction  (4,  7). 

The  regulation  of  MCPM  can  occur  through  growth  factor  stimu¬ 
lation  (4,  5)  or  requirements  for  eicosanoid  production."^  There  is  now 
increasing  evidence  to  suggest  that  the  activity  of  key  enzymes 
involved  in  MCPM  are  regulated  by  receptor  tyrosine  kinase  cascade 
downstream  of  the  ras/raf  interaction  (8,  9).  This  assertion  is  sup¬ 
ported  by  the  fact  that  activation  of  receptor  tyrosine  kinase  growth 
factor-mediated  signal  transduction  at  the  level  of  growth  factor,  ras, 
or  rqf  produces  an  enhanced  MCPM  (8,  9).  Because  signal  transduc¬ 
tion  events  and  genetic  alterations  involving  amplification  of  onco¬ 
genes  such  as  erbB2  play  a  crucial  role  in  the  development  of  the 
normal  breast,  carcinogenesis  of  its  epithelium,  and  progression  of 
breast  cancer,  it  is  possible  that  regulation  of  the  levels  of  choline- 
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containing  metabolites  may  be  linked  to  malignant  transformation  and 
progression  of  the  breast  epithelium. 

NMR  has  been  used  to  study  choline  phospholipid  metabolism  in 
cells  or  excised  tissues,  as  well  as  noninvasively  in  vivo  (10-13). 
Depending  on  the  experimental  conditions,  NMR  methods  can 
detect  either  individual  choline  phospholipid  metabolites  or  a  peak 
corresponding  to  total  choline-containing  metabolites.  Using 
NMR,  invasive  cancer  could  be  distinguished  from  benign  breast 
lesions  by  the  high  total  choline  phospholipid  metabolite  levels  in  the 
former  (10).  In  another  study,  increased  choline  phospholipid  metab¬ 
olite  levels  characterized  two  cancer  cell  lines  (MCF-7  and  T47D) 
compared  with  that  of  a  normal  HMEC  line;  there  were  no  distinct 
differences  in  high  energy  phosphates  and  the  rates  of  glucose  con¬ 
sumption  and  aerobic  glycolysis  (14).  These  studies  support  the  ex¬ 
istence  of  differences  in  phosphatidylcholine  metabolism  between 
normal  epithelial  cells  and  cancer  cells  in  vitro  and  between  benign 
and  malignant  cells  in  vivo.  The  possibility  of  differential  regulation 
of  MCPM  in  normal  versus  tumor  cells  suggests  a  diagnostic  role  for 
enhanced  MCPM  and  has  implications  for  therapeutic  intervention. 

Despite  the  indication  of  an  altered  MCPM  in  breast  cancer  cells, 
no  attempt  has  been  made  to  systematically  relate  the  multistep 
process  of  carcinogenesis  to  altered  MCPM  in  mammary  epithelial 
cells.  To  address  this  issue,  we  have  assessed  PCho,  GPC,  and  choline 
levels  in  a  number  of  epithelial  cell  lines  derived  from  reduction 
mammoplasty  (normal)  tissues  and  neoplastic  lesions  and  also  inves¬ 
tigated  the  effects  of  immortalization  and  oncogene  transformation  on 
MCPM.  Such  a  model  has  been  used  previously  by  other  workers  to 
evaluate  the  stepwise  progression  in  mammary  epithelium  from  nor¬ 
mal  to  malignant  phenotype  (15-20).  Our  data  suggest  that  phenotypic 
changes  in  MCPM  probably  commence  early  in  carcinogenesis  and 
may,  as  with  most  other  neoplastic  phenotypes,  be  regulated  by  an 
interplay  of  cellular  immortalization  and  oncogene  transformation. 

MATERIALS  AND  METHODS 

Cell  Lines.  HMECs  used  in  this  study  include  finite  life  span  HMEC 
strains  184  and  48,  derived  from  reduction  mammoplasty  tissues;  nontumori- 
genic  immortal  cell  lines  184A1  and  184B5,  derived  from  benzo(a)pyrene- 
treated  184  cells;  and  the  lMB5-erbB2  cell  line,  derived  from  184B5  by 
transfection  with  the  erbB2  oncogene.  All  of  the  above  cell  lines  were  obtained 
from  Dr.  Martha  Stampfer  (Lawrence  Berkeley  National  Laboratory,  Berkeley, 
CA)  and  cultured  in  MCDB  170  media  supplemented  as  described  previously 
(21,  22).  MCE- 12 A,  a  spontaneously  immortalized  cell  line  established  from 
MCF-12M  mortal  cells  (23),  was  obtained  from  American  Type  Culture 
Collection  (Rockville,  MD)  and  cultured  in  DMEM-Ham’s  FI 2  medium 
supplemented  as  described  previously  (23).  All  of  the  human  breast  cancer  cell 
lines  were  derived  from  pleural  effusions  in  patients  with  breast  cancer  and 
were  obtained  from  American  Type  Culture  Collection.  The  tumor-derived  cell 
lines  were  all  cultured  in  DMEM-Ham’s  F12  medium  supplemented  with  10% 
fetal  bovine  serum. 

Growth  Rate  and  Cell  Size.  The  growth  rate  of  the  cell  lines  used  in  this 
study  were  determined  using  the  MTT  assay  (24).  Briefly,  cells  (5  X  10^)  were 
plated  in  24- well  plates  in  1  ml  of  media  and  incubated  under  normal  culture 
conditions  for  up  to  6  days.  To  estimate  cell  number,  the  cells  were  incubated 
with  MTT  (Sigma  Chemical  Co.,  St.  Louis,  MO)  for  4  h.  MTT  was  then 
removed,  and  the  resulting  formazan  crystals  were  dissolved  in  1  ml  of  DMSO 
and  125  fjd  of  glycine  buffer  (pH  10.5;  Ref.  24).  The  UV  absorbance  of  the 
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Fig,  1.  Biosynthesis  and  hydrolysis  of  PtC.  Phosphorylation  of  choline  to  PCho  by 
choline  kinase  {CK)  is  the  first  step  in  the  biosynthesis  of  PtC.  PCho  is  then  converted  to 
PtC  via  intermediates  involving  the  rate-limiting  enzyme  CTP:  phosphocholine  cytidyl- 
transferase  (CT)  and  phosphocholine  transferase  {PCT).  Hydrolysis  of  PtC  is  effected  by 
three  major  PtC-specific  enzymes,  phospholipase  C  {PLQ,  phospholipase  D  {PLD),  and 
phospholipase  A,  and  A2  {PLA^  and  PLA^.  FFA,  free  fatty  acid. 


formazan  solution  was  recorded  at  553  nm  (Amax),  Four  replicates  were  used 
to  calculate  the  cell  doubling  time  for  each  cell  line.  Because  the  cells  had 
different  morphologies  and  diameters,  the  cell  size  was  determined  for  each 
cell  line  by  trypsinizing  the  cells  and  counting  the  diameter  of  20  random  cells 
using  an  optical  microscope. 

Extraction.  To  determine  the  choline  phospholipid  metabolite  content, 
cells  growing  in  culture  were  fed  with  fresh  media  24  h  before  extraction  and 
used  at  70-80%  confluency.  Cells  (10^  to  10^)  were  trypsinized,  washed  twice 
with  normal  saline,  and  homogenized  with  ice-cold  8%  perchloric  acid  (5  ml). 
The  homogenates  were  centrifuged  (15000  rpm  for  15  min  at  4°C),  and  the 
supernatants  were  neutralized  with  3  M  K2CO3/I  m  KOH  buffer.  The  samples 
were  again  clarified  by  centrifugation,  treated  with  ~50  mg  Chelex  (Sigma)  to 
remove  divalent  ions,  lyophilized,  and  resuspended  in  0.5  ml  of  D2O  for  NMR 
analysis.  Trimethylsilyl  propionate  (5  /xl)  was  used  as  an  internal  standard. 
NMR  spectra  of  the  extracts  were  acquired  on  a  11. 7T  Bruker  NMR  spec¬ 
trometer  with  a  5-mm  probe.  Fully  relaxed  spectra  (without  saturation  effects) 
were  obtained  using  the  following  acquisition  parameters:  30°  flip  angle,  6000 
Hz  sweep  width,  4.7  s  repetition  time,  32  K  block  size,  and  512  scans.  The  data 
were  analyzed  using  an  in-house  software.  Soft  Fourier  Transform  (P.  Barker, 
The  Johns  Hopkins  University).  PCho,  GPC,  and  total  choline-containing 
(PCho  +  GPC  +  choline)  metabolite  levels  were  determined  and  normalized 
to  cell  size.  Between  three  and  five  independent  extracts  were  analyzed  per  cell 
line. 

The  reason  for  normalizing  metabolite  levels  to  cell  size  was  due  to 
differences  in  cell  size  between  the  cell  lines  used  in  this  study.  This  neces¬ 
sitated  normalization  to  either  cell  size  or  protein  concentration.  The  former 
requires  fewer  cells  and  is  therefore  suited  to  experiments  with  mortal  cells, 
which  senesce  after  5  to  18  passages.  To  determine  concentrations,  peak 
amplitudes  for  choline  PCho,  GPC,  and  total  choline-containing  metabolites 
(PCho  +  GPC  +  choline)  were  compared  with  that  of  the  internal  standard 
TSP  according  to  the  equation: 


[metabolite]  =  Amplitude(n,e.aboiitc) 


[TSP] 

Amplitude^Tsp)  X  cell  number  X  cell  volume 


Statistical  Analysis.  Statistical  analysis  of  the  data  were  performed  using 
StatView'  II  version  1.04,  1991  (Abacus  Concepts,  Inc.,  Berkeley,  CA).  The 
statistical  significance  of  differences  in  metabolite  levels  between  cell  lines 
was  determined  using  the  Mann-Whitney  U  test.  Ps  of  <  0.05  were  considered 
to  be  significant. 

RESULTS 

Identification  of  Phospholipid  Metabolites  by  *H  NMR.  'H 

NMR  of  perchloric  acid  extracts  demonstrated  the  presence  of  three 
water-soluble,  choline-containing  [-N(CH3)3]  metabolites,  Le.,  cho¬ 
line,  PCho,  and  GPC  (Fig.  2).  These  metabolites  resonate  at  ~3.2  ppm 
downfield  of  the  internal  standard  and  chemical  shift  reference  TSP. 
Peak  assignments  were  performed  with  authentic  compounds.  Ten 
epithelial  cell  lines  of  mammary  origin  were  characterized  by  this 
method;  the  phenotype  and  cell  size  of  these  cell  lines  are  indicated  in 
Table  1. 

The  “GPC  to  PCho  Switch”  in  Mammary  Epithelial  Cells. 

Analysis  of  individual  choline  metabolites  uncovered  an  early  alter¬ 
ation  in  MCPM  that  was  linked  to  immortalization  and  malignant 
transformation,  the  “GPC  to  PCho  switch”  (Fig.  3).  In  Fig.  3  it  is 
evident  that  GPC  was  the  major  choline  metabolite  in  the  finite  life 
span  HMEC  strains  48  and  184.  Thus,  these  cells  showed  a  low 
PChoiGPC  ratio  of  <  1.  Immortalization  of  cells,  however,  resulted  in 
variable  effects.  The  spontaneously  immortalized  cell  line  MCF-12A 
showed  a  similar  phenotype  (PCho:GPC)  as  the  finite  life  span  cells; 
we  do  not  have  the  finite  life  span  cells  from  which  MCF-12A  was 
established  for  comparison.  In  contrast,  the  benzo(n)pyrene-immor- 
talized  cell  lines  showed  a  GPC  to  PCho  switch,  /.e.,  PCho  was  now 
the  major  choline  metabolite.  Of  interest,  the  two  immortal  lines 
derived  from  the  184  strain  showed  variable  degrees  of  this  altered 
MCPM;  184A1  had  a  higher  PCho: GPC  level  compared  with  184B5. 
Forced  overexpression  of  normal  erbB2  gene  into  184B5  cells  dra¬ 
matically  increased  the  PCho:GPC  ratio  in  this  cell  line.  The  altered 
GPC  to  PCho  switch  was  detected  in  all  breast  cancer  cell  lines 
analyzed. 

Breast  Cancer  Cells  Have  a  High  Choline  Content.  Fig.  4  shows 
that  there  was  a  gradual  increase  in  both  PCho  levels  and  total 
choline-containing  metabolite  levels  as  cells  progressed  from  normal 
to  malignant  phenotype  (normal  <  immortal  <  oncogene-trans¬ 
formed  <  tumor-derived).  GPC  levels  also  increased,  albeit  to  a  lesser 
extent  than  PCho  levels  and  total  choline-containing  metabolite  levels. 
It  is  worth  noting  that  despite  the  GPC  to  PCho  switch,  total  choline- 


a  GPC 


where  [metabolite]  is  concentration  of  the  metabolite  expressed  as  fmolZ/xm^, 
[TSP]  is  the  molar  concentration  of  TSP  used,  and  cell  volume  was 

calculated  from  the  radius  of  the  cell  according  to  the  equation,  volume  — 
4/3  X  ttt^.  For  this  equation  to  be  valid,  it  is  necessary  that  spectra  are  fully 
relaxed,  as  was  the  case  here,  or  to  correct  for  saturation. 


Fig.  2.  Typical  NMR  spectra  obtained  from  perchloric  acid  extracts  of  MCF-12A 
(fl)  and  MDA-MB-231  {b)  brca.st  cancer  cells  grown  in  culture.  The  spectra,  expanded  to 
show  the  choline-containing  metabolite  region,  represent  qualitative  differences  between 
the  two  cell  lines,  i.e.,  not  normalized  to  display  comparable  signal-lo-noisc  levels. 
Spectral  assignments  include  GPC  (3.234  ppm),  PCho  (3.225  ppm),  and  free  choline 
(Cho:  3.207  ppm). 
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Table  1  Phenotype  and  size  of  cell  lines  used  in  this  study 


Cell  types 

Phenotype'^ 

Cell  volume^ 
(XIO^  /mm^) 

Normal  HMECs 

184 

Senescent,  ADG 

36.5  ±  6.8 

48 

Senescent,  ADG 

11.1  ±  3.1 

Spontaneously  immortalized  HMECs 

MCF-12A 

Immortal,  ADG 

9.8  ±  1.8 

Benzo{fl)pyrene-immortalized  HMECs 

184A1 

Immortal,  ADG 

7.1  ±  0.9 

184B5 

Immortal,  ADG 

7.7  ±  0.9 

Oncogene-transformed  HMECs 

mB5-erbB2 

Immortal,  AIG,  forms  low 
frequency,  high  latency  tumors 

7.0  ±  0.8 

Breast  cancer  cells 

SKBR3 

AIG,  tumorigenic,  low 
metastaticity 

21.9  ±  7.3 

MCF7 

AIG,  tumorigenic,  low 
metastaticity 

6.7  ±  1.5 

MDA-MB-231 

AIG,  tumorigenic,  highly 
metastatic 

8.4  ±  1.4 

MDA-MB-435 

AIG,  tumorigenic,  highly 
metastatic 

3.5  ±  0.4 

^  ADG,  anchorage-dependent  growth;  AIG,  anchorage-independent  growth. 
*  Mean  ±  SE. 
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Fig.  3.  PCho:GPC  ratios  in  a  panel  of  cell  lines  representing  various  stages  of  breast 
carcinogenesis.  There  was  a  statistically  significant  difference  in  PCho:GPC  ratio 
(P  <  0.05)  between  finite  life  span  versus  tumor-derived  cells,  184  strain  versus  184A1 
cell  line,  and  184B5  versus  lS4B5-erbB2  cell  lines.  The  P  for  184  strain  versus  184B5  cell 
line  was  0.1.  Bars,  SE. 


ratio,  a  low  level  of  total  choline-containing  metabolites,  and  low 
PCho  compared  with  MDA-MB-435  cells. 

DISCUSSION 

Carcinogenesis  of  the  mammary  epithelium  occurs  though  a  mul¬ 
tistep  process  involving  genetic  alterations,  amplification  of  onco¬ 
genes,  and  loss  of  tumor  suppressor  function  (27,  28).  The  stages  of 


(a) 


PCho  levels  (fmol  /p  m  ^  cell) 


(b) 


u 


GPC  levels  (fmol  /p  cell) 


containing  metabolite  levels  and  PCho  levels  in  immortalized  cell 
lines  such  as  184A1,  which  was  nontumorigenic,  and  184B5-erbB2, 
which  exhibited  low  tumorigenicity,  were  significantly  lower  than  any 
of  the  tumorigenic  breast  cancer  cells. 

Are  the  High  Choline  Phospholipid  Metabolite  Levels  in  Breast 
Cancer  Cells  a  Function  of  Their  Rate  of  Cellular  Proliferation? 
It  is  generally  thought  that  the  increase  in  phosphomonoester  (mainly 
PCho  and  phosphoethanolamine)  metabolite  levels  in  cancer  cells  is 
due  to  their  intensified  cell  membrane  synthesis  to  cope  with  rapid 
growth  and  proliferation  (12,  25).  This  assertion  is  supported,  for 
instance,  by  the  work  of  Smith  et  al  (26),  where  an  increase  in  PCho 
and  a  decrease  in  GPC  correlated  strongly  with  tumor  growth  rate. 
Thus,  we  tested  the  hypothesis  that  increased  PCho:GPC  levels,  PCho, 
GPC,  or  total-choline  containing  metabolite  levels  may  be  the  result  of 
high  proliferation.  A  systematic  measurement  of  cell  doubling  time  in 
all  of  the  cell  lines  revealed  that  there  was  no  overall  correlation 
between  cell  doubling  time  and  PCho:GPC  ratio,  PCho  levels,  GPC 
levels,  or  total  choline-containing  metabolite  levels  (Fig.  5;  r  <  0.2; 
P  >  0.1).  For  instance,  MCF-12A  cells  exhibited  a  doubling  time  in 
culture  comparable  with  the  tumorigenic  cell  line  MDA-MB-435. 
However,  MCF-12A  cells  exhibited  a  significantly  lower  PCho: GPC 


Total  choline-containing  (fmol  /p  m  ^  cell) 
metabolite  levels 

Fig.  4.  PCho  levels  (a),  GPC  levels  ib),  and  total  choline-containing  metabolite 
(PCho  +  GPC  +  choline)  levels  (c)  in  a  panel  of  cell  lines  representing  various  stages  of 
breast  carcinogenesis.  There  was  a  statistically  significant  difference  in  total  choline- 
containing  metabolite  levels  and  PCho  levels  {P  <  0.05)  for  finite  life  span  cells  versus 
tumor-derived  cells;  184  strain  versus  184A1;  184B5  versus  \S4B5-erbB2;  MDA-MB- 
435  versus  MDA-MB-231,  MCF7,  and  SKBR3;  MDA-MB-231  versus  SKBR3;  and 
MCF7  versus  SKBR3.  There  was  a  statistically  significant  difference  in  GPC  levels  for 
finite  life  span  cells,  184A1, 184B5  and  \S5B5-erbB2,  and  SKBR3  versus  MDA-MB-435, 
MDA-MB-231,  MCF7,  and  MCF-12A  cells.  Bars,  SE. 


82 


MALIGNANT  TRANSFORMATION  AND  PHOSPHOLIPID  METABOLISM 


Fig.  5.  The  relationship  between  doubling  time 
and  PCho  levels  (a),  GPC  levels  (/?),  total  choline- 
containing  metabolite  levels  (PCho  +  GPC  +  cho¬ 
line;  c),  and  PChoiGPC  ratio  (d).  Doubling  times 
were  measured  by  the  MTT  assay  (see  “Materials  and 
Methods”)  and  increased  in  the  order  MCF-12A  >  48 
(mortal)  >  MDA-MB-435  >  184B5- 

erbB2  >  18465  >  MDA-MB-23I  >  184  (mor¬ 
tal)  >  184A1  >  MCF-7  >  SKBR3;  A,  MDA-MB- 
435;  T.  MDA-MB-231;  •,  MCF-7;  SKBR3;  ♦, 
MCF-12A;  A,  184  (mortal);  V,  48  (mortal);  O, 
184AI;  □,  184B5;  O,  mB5-erhB2.  Bars,  SE. 
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carcinogenesis  can  broadly  be  classified  as  transformation  of  normal 
cells  to  benign  hyperplasia  followed  by  atypical  hyperplasia,  which 
progresses  to  carcinoma  in  situ  and  finally  to  infiltrating  carcinoma 
with  or  without  metastasis  to  distant  sites.  In  vitro  models  based  on 
immortalization  and  oncogene  transformation  of  normal  HMECs  have 
been  developed  to  study  mammary  carcinogenesis  (15-20).  These 
immortalized,  oncogene-transformed,  and  cancer  cells  show  differ¬ 
ences  in  phenotypes  that  differentiate  normal  cells  from  immortal/ 
malignant  cells,  such  as  karyotype,  vimentin/uvomorulin  expression, 
responsiveness  to  transforming  growth  factor  j3,  telomere  length/ 
telomerase  expression,  activating  protein  1  transcription  factor  activ¬ 
ity,  as  well  as  anchorage-dependent/independent  growth,  and  ability 
to  form  tumors  in  immune-deficient  mice. 

We  have  investigated  the  association  between  malignant  carci¬ 
nogenic  processes  and  MCPM  by  monitoring  the  three  choline 
phospholipid  metabolites  (choline,  PCho,  and  GPC)  in  10  cell 
lines,  which  represent  different  stages  of  malignant  progression. 
Our  findings  suggest  that  normal  human  mammary  epithelium  has 
low  steady-state  levels  of  total  choline-containing  metabolites.  In 
addition  to  their  low  total  choline-containing  metabolite  levels,  we 
also  demonstrated  that  GPC  was  the  major  metabolite  in  the 
normal  HMECs.  A  GPC  to  PCho  switch  appeared  to  be  an  early 
phenotypic  change  during  carcinogenesis,  as  observed  in  benzo- 
(<3)pyrene-immortalized  cells  and  where  instead  of  GPC,  PCho 
became  the  major  choline  phospholipid  metabolites.  However, 
despite  this  “switch,”  total  choline-containing  metabolite  levels 
remained  low  in  these  immortalized  cells.  Transformation  of 
184B5  immortal  cells  by  forced  overexpression  of  the  erbB2 
oncogene,  however,  resulted  in  a  dramatic  increase  in  both  PCho: 
GPC  ratio  and  total  choline  levels  compared  with  the  benzo- 
(a)pyrene-immortalized  cells.  However,  total  choline-containing 
metabolites  and  PCho  levels  were  still  less  than  those  of  tumor- 


derived  cells.  erhB2  is  an  important  (proto)oncogene  that  is  am¬ 
plified  in  20-30%  of  breast  cancer  cases  and  is  associated  with 
poor  prognosis;  amplification  of  this  oncogene  is  thought  to  occur 
late  in  tumor  progression  (27-30).  Transformation  of  184B5  by 
this  gene  results  in  the  ability  of  these  cells  to  form  colonies  in 
semisolid  medium  and  to  form  small,  low  frequency  tumors  with 
high  latency  in  vivo  (16).  Our  data  with  erhB2  demonstrate  a  new 
and  heretofore  unknown  metabolic  role  for  erhB2  and  support  the 
possibility  that  growth  factor-mediated  activation  of  the  tyrosine 
kinase  cascade  (involving  receptor-gr/?  2-sos-ras-raf-l-MEK- 
MAPK)  can  lead  to  an  increase  in  PCho  levels  (6-9).  In  general, 
the  levels  and  expression  of  receptors  and  proteins  involved  in  the 
growth  factor  receptor-tyrosine  kinase  pathway  tend  to  increase 
with  malignancy.  For  instance,  levels  of  epidermal  growth  factor 
receptor  are  low  in  the  184  strain,  moderately  high  in  184A1, 
184B5,  and  \MB^5-erbB2  cells,  and  very  high  in  MDA-MB-231 
cells  (19,  31).  In  addition,  Daly  et  al.  (32)  reported  up-regulation 
of  grb2  mRNA/protein  and  the  ras  signaling  pathway  in  MCF-7 
and  MDA-MB-231  cells  compared  with  normal  HMECs. 

All  of  the  breast  tumor  cell  lines  showed  the  GPC  to  PCho  switch. 
In  addition  to  this  switch,  all  breast  tumor  cells  showed  significantly 
higher  total  choline-containing  metabolite  levels  (P  <  0.05).  The 
increased  total  choline-containing  metabolite  levels  were  mainly  due 
to  an  increase  in  PCho  levels  and,  to  a  lesser  and  variable  extent,  an 
increase  in  GPC  levels.  There  was  a  gradual  increase  in  both  total 
choline-containing  metabolite  levels  and  PCho  levels  as  the  cells 
acquired  malignant  phenotype  (normal  <  immortal  <  oncogene- 
transformed  <  tumor-derived),  with  the  highly  invasive  metastatic 
cell  lines  showing  the  highest  levels.  The  high  total  choline  content  in 
the  tumorigenic  cells  may  be  related  to  the  multiple  genetic  changes 
that  are  associated  with  the  multistep  process  of  carcinogenesis  (28) 
and  may  explain  the  progressive  ability  of  these  cells  to  gain 
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alichorage-independent  growth,  form  primary  tumors  in  immune 
compromised  mice,  and  finally  to  metastasize.  Our  studies  confirm 
the  work  of  Ting  et  al.  (14),  who  showed  for  a  limited  number  of 
cell  lines  that  levels  of  choline-containing  metabolites  were  low  in 
a  normal  mammary  epithelial  strain  and  high  in  two  tumor- derived 
cell  lines.  Our  results  also  support  recent  clinical  observations  that 
the  total  choline  peak  is  higher  for  malignant  lesions  than  for 
benign  ones  (10). 

It  has  been  postulated  that  the  rapid  growth  and  proliferation  of 
cancer  cells  and  increased  membrane/fatty  acid  requirements  may  be 
responsible  for  the  high  choline  phospholipid  metabolite  levels  in 
cancer  versus  normal  tissues  (12,  25,  26);  the  same  argument  could  be 
made  for  benign  lesions  versus  invasive  cancers.  However,  the  data 
presented  here  and  that  of  Ting  et  al  (14)  show  that  choline-contain¬ 
ing  metabolite  levels  remain  low  in  normal  HMECs  in  culture  when 
the  cells  are  proliferating  at  approximately  similar  rates  as  tumor- 
derived  cells  and  suggest  that  although  proliferation-related  changes 
may  occur  (26),  the  rate  of  proliferation  per  se  cannot  completely 
account  for  the  increased  choline  phospholipid  metabolism.  In  this 
study,  we  have  demonstrated  that  an  alteration  in  MCPM  is  linked  to 
malignant  transformation  and  progression  of  mammary  epithelium. 
Presently,  the  exact  mechanisms  underlying  the  altered  metabolism 
are  unknown.  Possible  mechanisms  include  activation  of  enzymes 
involved  in  MCPM,  such  as  via  enhanced  receptor  tyrosine  kinase 
cascade  (9),  or  differential  induction  of  choline  kinase  isozymes,  as 
reported  previously  for  carcinogen-treated  rat  liver  (33).  Other  possi¬ 
ble  mechanisms  that  need  to  be  investigated  include  amplification  of 
choline  kinase,  phospholipase  C,  phospholipase  D,  and  phospholipase 
A  genes  during  carcinogenesis. 

To  conclude,  the  major  finding  to  emerge  from  the  present  study 
is  that  choline  phospholipid  metabolite  levels  progressively  in¬ 
crease  in  cultured  HMECs  as  cells  become  more  malignant.  We 
therefore  propose  that  carcinogenesis  in  human  breast  epithelial 
cells  results  in  progressive  alteration  of  membrane  choline  phos¬ 
pholipid  metabolism.  This  work  is  relevant  to  diagnosis  of  breast 
cancer  and  also  provides  a  rationale  for  selective  pharmacological 
intervention. 
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In  Vivo  Imaging  of  Extracellular  pH  Using  MRSI 
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Tumor  pH  is  physiologically  important  since  it  influences  a 
number  of  processes  relevant  to  tumorigenesis  and  therapy. 
Hence,  knowledge  of  localized  pH  within  tumors  would  contrib¬ 
ute  to  understanding  these  processes.  The  destructiveness, 
poor  spatial  resolution,  and  poor  signal-to-noise  ratio  (SNR)  of 
current  technologies  (e.g.,  microelectrodes,  magnetic  reso¬ 
nance  spectroscopy)  have  limited  such  studies.  An  extrinsic 
chemical  extracellular  pH  (pHe)  probe  is  described  that  is  used 
in  combination  with  ^H  magnetic  resonance  spectroscopic 
imaging  to  yield  pHe  maps  with  a  spatial  resolution  of  1  x  1  x  4 
mm^.  The  principle  of  the  technique  Is  demonstrated  on  a 
phantom.  Further  data  are  shown  to  demonstrate  its  application 
in  vivo,  and  results  agree  with  previously  reported  pH  values. 
The  accuracy  of  the  reported  pH  measurements  is  <0.1  pH 
units,  as  derived  from  a  detailed  analysis  of  the  errors  associ¬ 
ated  with  the  technique,  the  description  of  which  is  included. 
Magn  Reson  Med  41 :743-750, 1 999.  ©  1 999  Wiley-Liss,  Inc. 
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Since  the  discovery  of  lactic  acid  production  in  tumors 
more  than  50  years  ago  (1),  it  has  generally  been  assumed 
that  the  pH  of  tumors  is  acidic.  Indeed,  numerous  microelec¬ 
trode  measurements  have  shown  that  extracellular  tumor 
pH  (pHe)  is  acidic  (2).  This  acidic  pHg  of  tumors  has  been 
confirmed  with  less  invasive  magnetic  resonance  spec¬ 
troscopy  (MRS)  measurements  (3).  Although  the  intracellu¬ 
lar  pH  (pHi)  of  tumors  remains  neutral  to  alkaline  (4,5),  it  is 
somewhat  influenced  by  the  pHe  (6). 

An  acidic  pHg  of  tumors  is  physiologically  important 
since  it  influences  a  number  of  processes  relevant  to 
carcinogenesis  and  therapy.  Knowledge  of  localized  pH 
within  tumors,  both  intra-  and  extracellular,  would  allow 
more  detailed  study  of  these  processes  and  relate  them  to 
intratumoral  pH  heterogeneity.  For  example,  it  has  been 
found  that  low  pHe  in  vitro  causes  tumorigenic  transforma¬ 
tion  of  primary  Syrian  hamster  embryo  cells  (7)  and  can 
lead  to  chromosomal  rearrangements  in  Chinese  hamster 
embryo  cells  (8,9).  Furthermore,  culturing  cells  at  low  pH 
causes  them  to  be  more  invasive  in  vitro  (10)  and  meta¬ 
static  in  vivo  (11).  Finally,  the  orientation  of  the  pH 
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gradient  across  the  cell  membrane  may  influence  cell  drug 
resistance  (6,12) 

Previously  reported  measurements  of  extracellular  pH 
using  either  microelectrodes  or  MRS  of  3-aminopro- 
pylphosphonate  (3-APP)  (3)  have  drawbacks.  Microelec¬ 
trodes  are  invasive  and  can  destroy  the  membrane  integ¬ 
rity,  thereby  disrupting  the  mechanism  for  maintaining  the 
pHg.  ^^P  MRS  does  not  suffer  this  drawback  and  has  the 
additional  advantage  of  permitting  simultaneous  measure¬ 
ments  of  intracellular  pH.  However,  the  limited  sensitivity 
of  ^^P  MRS  allows  measurements  of  pHe  only  from  rela¬ 
tively  large  tissue  volumes.  Hence,  ^^P  MRS  provides 
measurements  of  pH  ranges  rather  than  different  pH  values 
for  discrete  spatial  locations  (13). 

The  use  of  ^H  MRS,  inherently  more  sensitive  than  ^^P 
MRS,  would  allow  measurements  of  pH  over  smaller  tissue 
volumes.  For  example,  the  imidazole  protons  of  histidine 
have  long  been  useful  as  intracellular  pH  indicators  in 
NMR  (14,15). 

Rabenstein  and  Isab  (16)  first  proposed  using  imidazoles 
as  extrinsic  pHe  indicators.  Gil  et  al  (17)  suggested  several 
modifications  of  the  basic  structure  of  the  imidazole  mol¬ 
ecule  to  improve  its  performance  as  an  extrinsic  pH  probe. 
To  date,  the  most  promising  candidate  for  a  ^H  nuclear 
magnetic  resonance  (NMR)-sensitive  pHe  indicator  is  the 
H2  resonance  of  (±)  2-imidazole-l-yl-3-ethoxycarbonyl 
propionic  acid  (lEPA),  which  has  been  shown  to  remain  in 
the  extracellular  environment  (17)  Results  from  toxicity 
studies  and  preliminary  in  vivo  data  using  lEPA  have  been 
reported  previously  (18,19). 

Here,  we  present  both  phantom  and  in  vivo  data  that 
demonstrate  the  feasibility  of  localized,  multi-voxel  pHe 
measurements  using  ^H  magnetic  resonance  spectroscopic 
imaging  (MRSI).  Furthermore,  a  detailed  discussion  of  the 
sources  and  (where  possible)  the  quantification  of  errors 
associated  with  the  proposed  measurement  technique  is 
included. 

MATERIALS  AND  METHODS 

Titration  of  lEPA 

Titration  data  for  lEPA  were  obtained  from  bovine  serum 
containing  20  mM  lEPA,  to  which  20  mM  EDTA  and  10% 
v/v  deuterium  oxide  (D2O)  were  added  for  preservation 
and  shimlock  purposes,  respectively.  As  a  chemical  shift 
reference,  10  mM  3-(trimethylsilyl)  tetradeutero  sodium 
propionate  (TSP)  was  also  added.  The  pH  of  the  solution 
was  adjusted  at  37°C  using  5  N  solutions  of  NaOH  and  HCl 
to  19  pH  values  between  4.7  and  8.0  and  the  chemical  shift 
of  the  H2  lEPA  resonance  was  measured  at  37°C  in  the 
spectrum  for  each  pH  value.  Figure  1  shows  the  chemical 
structure  and  a  high-resolution  (9.4  T)  ^H  spectrum  of  lEPA 
at  pH  7.0. 
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The  Henderson-Hasselbalch  equation  for  lEPA  is  given 

by 

pH-pK.-l„g.4^ii^),  lEq^l 

\Obase  “  ^obs/ 

where  pH  is  the  observed  pH,  pKa  is  the  negative  logio  of 
the  acid  dissociation  constant,  and  8obsi  Kcid^  and  8base 
are  the  chemical  shifts  of  the  observed,  the  fully  proto- 
nated,  and  the  fully  deprotonated  lEPA  peaks,  respectively. 
The  experimental  data  were  fitted  to  Eq.[l]  through  nonlin¬ 
ear  regression  (Sigmaplot  3.0,  Jandel),  and  the  results  are 
shown  in  Fig.  2.  The  fit  is  characterized  by  the  following 
set  of  parameters  and  their  respective  standard  errors: 

pK„  =  6.49  ±  0.02; 

=  8.92  ±  0.01; 

8be,e  =  7.77  ±  0.01. 

MRSI  Experiments 

All  MRSI  experiments  were  performed  on  a  4.7  T  GE 
Omega  scanner,  using  the  BASS  ALE  sequence  (20)  and  a 
two-turn  ^H  surface  coil.  Two-dimensional  (2D)  water- 
suppressed  MRSI  data  were  obtained  from  a  4  mm  coronal 
slice,  phase-encoding  dimensions  of  2  X  2  mm^  and  16  X 
16  steps.  The  fields  of  view  (FOVs)  in  all  imaging  and  MRSI 
studies  were  32  X  32  mm^. 

A  three-compartment  phantom  was  constructed,  consist¬ 
ing  of  two  5  mm  NMR  tubes  inside  a  larger,  12  mm  glass 
cylinder.  The  two  inner  tubes  were  filled  with  20  mM  lEPA 
at  pH  values  of  6.5  and  7.5,  respectively.  The  third,  outer 
compartment  was  filled  with  saline.  ^H  MRSI  data  were 
acquired  with  a  TE  of  100  msec,  TR  of  2  sec,  512  points, 
spectral  width  (SW)  of  6000  Hz,  4  averages.  A  separate  data 
set  was  obtained  without  water  suppression  to  correct  for 
susceptibility  effects. 


FIG.  1 .  Structure  and  9.4  T  NMR  spectrum  of  lEPA  in  solution  (pH 
7.0).  The  pH-sensitive  resonance  and  the  C2  proton  from  which  it 
originates  are  indicated  by  asterisks.  The  broad  resonance  in  the 
center  of  the  spectrum  is  due  to  residual  water. 


FIG.  2.  Chemical  shift  of  the  H2-IEPA  resonance  relative  to  TSP,  as  a 
function  of  pH  in  bovine  serum  at  37°C.  Measured  data  are  shown  as 
points,  and  the  curve  results  from  a  nonlinear  regression  of  the 
experimental  data  to  Eq.  [2]. 

In  vivo  studies  were  performed  on  human  breast  cancer 
cells  (MCF-7  and  MDA"^^-435),  grown  in  the  mammary  fat 
pad  of  severe  combined  immunodeficient  (SCID)  mice. 
Tumor  volumes  studied  were  300-600  mm^.  The  mice 
were  anesthetized  with  ketamine  (50  mg/kg)  and  aceproma- 
zine  (5  mg/kg),  and  0.15  ml  of  a  310  mM  lEPA  solution  was 
injected  intra-peritoneally.  The  animals  were  immobilized 
on  a  home-built  MR  cradle,  and  body  temperature  was 
maintained  using  a  warm  water  blanket.  ^H  MRSI  data 
were  acquired  with  a  TE  of  32  msec,  TR  of  2  sec,  128  points, 
SW  of  6000  Hz,  and  8  averages.  The  total  measurement 
time  was  34  min.  A  separate  data  set  was  obtained  without 
water  suppression  to  correct  for  susceptibility  effects. 

All  data  were  processed  using  an  MRSI  processing 
program,  developed  in-house  and  written  in  Interactive 
Data  Language  (IDL,  Research  Systems,  Boulder,  CO).  The 
raw  ^H  MRSI  data  were  Hanning  filtered  in  ic-space,  and 
zero  filled  once  in  both  k-space  and  in  the  time  domain 
prior  to  3D  Fourier  transformation. 

RESULTS  AND  DISCUSSION 

Phantom  Data 

The  chemical  shift  of  the  H2  resonance  of  lEPA  was 
determined  for  the  maximum  peak  occurring  between  7.77 
and  8.92  ppm  (the  H2  proton  titratable  range)  in  each  of  the 
MRSI  spectra.  Chemical  shifts  in  each  voxel  were  refer¬ 
enced  to  unsuppressed  water  at  4.7  ppm  to  correct  for 
susceptibility  effects,  and  converted  to  pH  values  using  Eq. 
[1].  The  pH  values  from  each  voxel  were  color-encoded. 
Peaks  lower  than  10%  of  the  maximum  lEPA  intensity 
were  discarded,  and  results  are  shown  in  Fig.  3. 

It  is  clear  from  Fig.  3b  that  the  pH  imaging  technique 
allows  the  two  compartments  of  the  phantom  to  be  distin¬ 
guished.  Furthermore,  the  reported  pH  values  are  in  good 
agreement  with  the  calibration  values  of  6.5  and  7.5  pH 
units.  Note  that  the  pH  6.5  compartment  is  spread  over 
more  voxels  than  the  pH  7.5  compartment.  Over  the  pH 
range  5-8,  the  T2  of  the  C2  proton  of  lEPA  decreases  from 
101  ±  12  msec  to  61  ±  4  msec.  In  spin-echo  sequences  like 
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FIG.  3.  a:  Coronal  image  through  the  three  compartment  phantom,  showing  the  two  tubes  with  pH  6.5  (bottom  left)  and  pH  7.5  (top  right).  The 
bright  spot  is  a  capillary  filled  with  Gd-DTPA,  used  for  orientation  purposes,  b:  pH  map  and  legend  for  the  phantom.  The  pixel  size  in  the  pH 
map  is  1  mm^. 


BASS  ALE,  resonances  with  a  longer  T2  will  have  a  higher 
signal  intensity  than  resonances  with  a  shorter  T2.  Since 
the  T2  of  the  C2  proton  of  lEPA  decreases  with  increased 
pH,  the  H2-IEPA  resonance  from  the  lower  pH  compart¬ 
ment  will  have  a  higher  intensity  than  the  H2-IEPA  reso¬ 
nance  from  the  higher  pH  compartment.  The  degree  of 
contamination  due  to  the  point  spread  function  (PSF)  is 
proportional  to  the  intensity  of  the  signal  from  which  the 
contamination  arises  (21).  As  a  result  the  contamination  for 
the  pH  6.5  compartment  extends  over  more  neighboring 
voxels  compared  with  the  pH  7.5  compartment.  Further¬ 
more,  this  T2  weighting  could  lead  to  an  underestimate  of 
the  reported  pH  value.  The  decrease  of  T2  with  increasing 
pH  is  consistent  with  the  known  increase  in  the  exchange 
rate  with  solvent  of  protons  from  primary  and  secondary 
amines  and  amides  (22).  Increasing  solvent  exchange  rate 
would  tend  to  decrease  T2  values  of  the  imidazolic  protons 
because  of  an  increased  magnetization  transfer  to  the 
surrounding  solvent.  Regardless  of  the  exact  underlying 
mechanisms  of  the  T2  dependence  on  pH,  it  is  clear  that  ^H 
MRSI  of  the  H2  proton  of  lEPA  allows  us  to  produce  a  pH 
map  of  the  phantom  and  distinguish  between  the  two 
compartments,  as  shown  in  Fig.  3. 

In  Vivo  Data 

An  illustration  of  the  quality  of  our  in  vivo  MRSI  data,  an 
array  of  7  X  7  spectra,  is  shown  in  Fig.  4.  These  data 
demonstrate  that  the  H2  resonance  of  lEPA  is  observable 
with  good  SNR  in  15-20  voxels  in  the  tumor.  Pre-injection 
spectra  contain  a  small  resonance  in  the  7-9  ppm  region 
corresponding  to  endogenous  imidazole  (data  not  shown). 
This  resonance  was  smaller  than  that  in,  for  example,  voxel 
(9,12),  and  did  not  affect  subsequent  measurements. 

Some  voxels,  in  particular  [10,8],  [10,9],  and  [10,10], 
bracketed  in  Fig.  4,  show  the  presence  of  a  complex 
H2-IEPA  resonance,  which  may  result  from  the  presence  of 
two  or  more  pools  of  distinct  pHe  values.  In  all  spectra,  the 


H2-IEPA  resonance  is  rather  broad.  The  short  T2*  of  this 
resonance  in  vivo  does  not  allow  a  finer  sampling  of  the 
spectrum  without  ramifications  for  the  SNR.  At  present, 
the  optimum  spectral  sampling  interval  is  0.23  ppm/point, 
which  is  coarse.  Clearly,  this  is  an  aspect  of  the  work  that 
needs  improvement. 

For  construction  of  the  pHe  map  from  the  in  vivo  data, 
the  chemical  shift  of  the  H2  resonance  of  lEPA  was 
determined  using  a  center  of  gravity  (COG)  calculation  (23) 
over  the  titratable  chemical  shift  range  for  this  proton.  The 
calculation  was  weighted  by  the  shape  of  the  titration 
curve  to  correct  for  its  non-linearity  (24).  Over  the  titratable 
range  of  the  lEPA  C2  proton,  the  weighting  function  is  the 
partial  derivative  of  the  observed  chemical  shift  with 
respect  to  the  pH  for  Eq.  [l],  and  is  given  by: 

^^obs  2.3(83^.}^  —  S^bsK^obs  “  ^base^  21 

apH  (8acid  -  hase)  ^ 

Outside  the  titratable  range,  the  weighting  function  is  zero. 
As  before,  chemical  shift  values  were  references  to  unsup¬ 
pressed  water  at  4.7  ppm  to  correct  for  susceptibility 
effects  and  the  calculated  pH  values  were  color-encoded. 
Peaks  lower  than  10%  of  the  maximum  lEPA  intensity 
were  discarded.  Results  are  shown  in  Fig.  5. 

The  intensity  of  the  H2-IEPA  resonance  is  the  largest  in 
the  center  of  the  tumor  (Fig.  5a).  Since  tumor  perfusion  is 
generally  worse  in  the  center  of  a  tumor  than  in  the 
periphery,  one  would  expect  to  find  a  higher  concentration 
of  lEPA  in  the  periphery  of  the  tumor  compared  with  the 
center.  However,  the  observed  intensity  distribution  is 
likely  to  be  caused  by  the  application  of  the  Hanning  filter 
to  the  raw  data.  This  filter  reduces  the  signal  intensity  at 
the  edges  of  the  FOV,  to  reduce  ringing  artifacts.  As  a  result, 
data  in  the  center  of  the  FOV  appear  to  be  at  a  higher 
intensity. 


FIG.  4.  Array  of  7  X  7  in  vivo  ‘•H  MRSI  spectra  used  for  the  determination  of  pHe  in  the  tumor.  The  H2-IEPA  resonance  is  clearly  visible  in  the 
majority  of  the  spectra.  The  displayed  chemical  shift  range  is  10-6  ppm.  Note  that  for  these  spectra,  the  voxel  size  is  2  x  2  mm^.  For 
construction  of  the  pHg  map,  data  were  interpolated  in  /c-space  to  yield  a  voxel  size  of  1  mm^.  The  boxed  spectra,  [1 0,8],  [1 0,9],  and  [1 0,1 0]  are 
discussed  in  detail  in  the  text. 


FIG.  5.  a:  Intensity  map  of  the  H2“resonance  of  lEPA  in  a  coronal  slice  through  an  MDA^M35  tumor,  b:  Corresponding  pHe  map.  In  both  the 
intensity  and  pHe  map,  the  pixel  size  is  1  mm^. 
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As  expected,  the  data  show  a  good  correlation  between 
the  signal  intensity  distribution  and  the  extent  of  the  pH^ 
map  (cf.  Fig.  5a  and  b).  Across  the  tumor,  the  observed 
extracellular  pH  in  the  tumor  is  heterogeneous  and  acidic 
(6.44  <  pHe  <6.79),  with  a  volume  average  pHg  of  6.64  ± 
0.07.  This  is  in  good  agreement  with  previous  studies  (13). 
Furthermore,  it  is  seen  that  there  are  several  discrete  areas 
of  fairly  homogeneous  pHe  within  the  tumor,  with  areas 
between  4  and  9  mm^  (Fig,  5b).  The  largest  continuous  area 
within  the  tumor  appears  to  have  an  average  pHe  of  6.6,  as 
read  from  the  legend. 

Error  Analysis 

The  pH  ranges  as  reported  in  Figs.  3b  and  5b  have 
associated  errors  resulting  from  artifacts  in  the  measure¬ 
ment  technique.  The  errors  associated  with  the  calibration 
of  the  pHe  are  relatively  small  (see  above),  so  that  the 
measurement  technique  proposed  above  has  two  remain¬ 
ing  sources  of  error  to  be  considered.  First,  the  COG 
method  used  in  determining  the  peak  position  will  have  an 
accuracy  dependent  on  the  line  width  and  SNR  of  the  peak 
under  observation.  Second,  signal  modulation  due  to  the 
PSF,  an  inherent  source  of  error  in  any  MRSI  measurement 
(21),  will  cause  signals  to  ‘'bleed”  into  adjacent  voxels. 
This  voxel  bleed  artifact  depends  on  the  intensity  and 
chemical  shift  differences  existing  between  adjacent  vox¬ 
els.  Intensity  differences  may  be  caused  by  changes  in  the 
lEPA  relaxation  times  with  pH  or  by  differences  in  lEPA 
concentration  between  adjacent  voxels.  Chemical  shift 
differences  between  voxels  may  be  caused  by  pHg  differ¬ 
ences.  These  sources  of  error  will  be  described  in  more 
detail  below. 

Accuracy  of  the  COG  Method 

To  study  the  accuracy  of  the  COG  method  for  determina¬ 
tion  of  pH  from  our  spectra,  we  simulated  C2  lEPA 
resonances  for  pH  values  6.0,  6.5,  7.0,  and  7.5,  with  a 
spectral  resolution  of  0.23  ppm/point.  Resonances  were 
simulated  with  different  linewidths  as  defined  by  the  full 
width  at  half-maximum  (FWHM).  Gaussian  noise  was 
added  to  the  lEPA  resonance  to  yield  SNRs  of  approxi¬ 
mately  20,  10,  and  5.  The  pH  as  determined  by  the  COG 
method  was  compared  with  the  pH  as  entered  into  the 
simulation  and  the  absolute  error  gives  a  measure  of  the 
accuracy  of  the  COG  method  in  determining  the  pH. 
Results  are  shown  in  Fig.  6. 

The  error  in  the  reported  pH  increases  with  decreasing 
SNR,  for  all  linewidths,  as  expected.  Furthermore,  the  pH 
is  reported  more  accurately  near  the  pKa  of  the  C2  lEPA 
proton,  which  is  6.49.  This  is  due  to  the  nonlinearity  of  the 
titration  curve. 

It  is  also  observed  that  the  error  in  the  reported  pH 
increases  as  the  H2-IEPA  resonance  linewidth  decreases 
(see  Fig.  6a  and  6c).  This  is  directly  related  to  the  spectral 
resolution  in  the  simulated  data  (0.23  ppm/point).  As 
peaks  narrow,  the  lines  are  defined  by  fewer  and  fewer 
points,  and  with  the  FWHM  approaching  the  spectral 
resolution,  the  determination  of  the  position  of  the  maxi¬ 
mum  peak  intensity  by  the  COG  method  becomes  stochas¬ 
tic,  However,  for  reasonable  SNRs  (>5)  and  broader  peaks 
(FWHM  >0.16  ppm),  as  observed  in  our  in  vivo  spectra. 


the  accuracy  of  the  COG  method  is  better  than  0.1  pH  units 
in  the  range  from  pH  6.0  to  pH  7.0. 

Modeling  PSF  Artifacts 

Wang  et  al  (25)  have  proposed  a  model  to  calculate  voxel 
bleed  artifacts  due  to  concentration  differences  alone. 
Here,  we  have  modified  their  model  to  study  the  combined 
effects  of  inhomogeneities  in  both  pHe  and  lEPA  concentra¬ 
tion. 

For  simulations,  an  8  X  8  grid  of  lEPA  free  induction 
decays  (FIDs)  was  established,  with  voxel  coordinates 
running  from  [0,0]  to  [7,7].  In  the  descriptions  below,  voxel 
[4,4]  was  used  as  the  “center  voxel”.  Initial  simulations 
demonstrated  that  the  worst  contamination  occurs  in  the 
center  voxel  when  the  signal  intensity  ratio  is  maximum, 
i.e.,  an  empty  voxel  surrounded  by  voxels  of  non-zero  lEPA 
concentration.  Also,  the  contamination  is  worst  between 
nearest  neighbor  voxels,  i.e.,  the  main  contributors  to  the 
contamination  in  the  center  voxel  are  the  8  voxels  immedi¬ 
ately  surrounding  it.  These  results  are  in  agreement  with 
those  of  Wang  et  al  (25). 

Quantitative  estimates  of  the  actual  degrees  of  lEPA 
concentration  and  pHe  inhomogeneities  in  vivo  are  re¬ 
quired  before  continuing  analysis  of  their  effects  on  the 
precision  of  the  proposed  pH  imaging  technique. 

Estimating  lEPA  Concentration  Inhomogeneity 

Since  voxel  bleed  artifacts  prevent  direct  derivation  of 
quantitative  information  on  lEPA  concentration  inhomoge¬ 
neity  from  the  MRSI  data,  studies  of  the  uptake  of  dimeglu- 
mine  gadopentetate  (Gd-DTPA,  Schering,  Germany)  in 
MCF-7  tumors  were  monitored  by  T2”Weighted  imaging  as 
a  model  for  the  lEPA  concentration  distribution.  lEPA  and 
Gd-DTPA  are  comparable  in  terms  of  partition  coefficient 
and  charge.  Both  molecules  are  about  the  same  size,  and 
their  charges  are  -2  for  Gd-DTPA  (Gd  +3,  DTPA  -5)  and  0  to 
-1  for  lEPA  (acid  group  -1,  imidazolic  nitrogen  +0  to  +1). 

One  precontrast  (t  =  0)  and  16  Gd-DTPA-enhanced 
images  (17  mm  surface  coil,  FOV  25  mm^,  128  X  128 
pixels,  slice  thickness  4  mm)  were  obtained  at  regular 
intervals  over  the  course  of  110  min  from  three  consecutive 
sagittal  slices  through  MCF-7  tumors.  Gd-DTPA-enhanced 
images  were  subtracted  from  the  precontrast  image  to  yield 
images  of  the  Gd-DTPA  enhancement  only,  which  to  a  first 
approximation  is  linear  with  Gd-DTPA  concentration  (26). 
These  subtraction  images  were  rebinned  to  16  X  16  pixels 
to  yield  a  pixel  size  of  1.56  mm^,  approximating  the  2  mm^ 
spatial  resolution  used  when  acquiring  the  lEPA  MRSI 
data.  Examples  of  a  high-resolution  (128  X  128)  and 
rebinned  (16  X  16)  images  of  Gd-DTPA  enhancement  are 
shown  in  Fig.  7. 

In  each  of  the  rebinned  subtraction  images,  6  pixels  were 
identified  that  represent  the  Gd-DTPA  distribution  in  the 
tumor  and,  for  these  6  pixels,  the  intensity  ratios  between 
nearest  neighbors  were  calculated.  The  highest  and  the 
average  nearest  neighbor  intensity  ratios  were  determined 
for  each  image  and  plotted  as  a  function  of  time,  as  shown 
in  Fig.  8.  It  can  be  seen  that  once  a  steady-state  Gd-DTPA 
concentration  has  been  reached  (after  6~8  min),  the  highest 
and  average  intensity  ratios  stabilize  at  about  1:3.0  and 
1:1.5,  respectively.  Using  these  values  as  a  measure  of  the 
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FIG.  6.  The  accuracy  in  the  pH  as  reported  by  the  COG  method  as  a 
function  of  actual  pH  and  the  signal-to-noise  ratio  (SNR)  in  simulated 
data.  Results  are  shown  for  different  linewidths  as  defined  by  the 
full-width  at  half-maximum  (FWHM)  for  (a)  FWHM  0.32  ppm,  (b) 
FWHM  0.16  ppm.  and  (c)  FWHM  0.08  ppm. 
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FIG.  7.  a:  High-resolution  (128  x 
128  pixels)  subtraction  image  of 
Gd-DTPA  uptake  in  an  MCF-7  tu¬ 
mor  in  vivo,  b:  The  same  image, 
rebinned  to  16  x  16  pixels.  The 
resulting  pixel  size  (1 .56  mm^)  ap¬ 
proximates  the  2  mm^  spatial  reso¬ 
lution  used  when  acquiring  lEPA 
MRSI  data.  In  both  images,  the 
white  box  indicates  the  location  of 
the  tumor.  The  six  pixels  entirely 
within  the  tumor  in  b  were  used  to 
calculate  concentration  differences 
between  adjacent  voxels. 
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FIG.  8.  The  nearest  neighbor  highest  (a)  and  average  (b)  intensity 
ratio  as  a  function  of  time  following  Gd-DTPA  injection  in  MCF-7 
tumors. 


lEPA  intensity  ratios  between  adjacent  voxels,  a  simulated 
lEPA  distribution  was  set  up  over  the  previously  men¬ 
tioned  grid  of  8  X  8  voxels. 

Quantifying  pHe  Inhomogeneity 

Information  on  the  maximum  spatial  pHe  gradient  ex¬ 
pected  for  adjacent  voxels  in  the  MRSI  data  was  obtained 
through  the  work  of  Helmlinger  et  al  (27),  who  reported  a 
0.7  unit  pHe  difference  between  well-oxygenated  (near  the 
blood  vessel)  and  nonoxygenated  regions  (away  from  the 
blood  vessel)  in  tumors,  over  a  distance  of  400  ixm.  This  is 
in  good  agreement  with  work  by  Raghunand  et  al  (13),  who 
found  that  85%  of  tumor  pHe  values  are  within  ±0.4  pH 
units  of  the  mean  pHg.  To  observe  a  pH  difference  of  0.7  pH 
units  between  adjacent  16  mm^  voxels  in  our  MRSI  data,  it 
is  necessary  to  have  a  fully  oxygenated  voxel  (i.e,,  full  of 
capillaries)  immediately  next  to  nonoxygenated  voxels 
(i.e.,  no  capillaries).  It  would  be  possible  to  find  com¬ 
pletely  nonoxygenated  voxels  in  vivo,  for  example,  in 
necrotic  areas.  However,  since  in  general  tumor  cells 
divide  faster  than  their  surrounding  vasculature,  the  pres¬ 
ence  of  a  voxel  consisting  entirely  of  capillaries  would  be 
less  likely.  Therefore,  we  have  chosen  to  work  with  the 
average  pHg  range  observed  by  both  Helmlinger  et  al  and 


Raghunand  et  al  (i.e.,  0.4  pH  unit),  as  the  maximum  pH 
gradient  between  adjacent  voxels. 

Combined  Effect  oflEPA  Concentration  and  pHe 
Inhomogeneity  on  COG  Precision 

The  maximum  pH  differences,  with  signal  intensity  ratios 
of  either  1:3.0  or  1:1.5,  were  inserted  into  the  simulation 
program.  The  reported  pH  for  the  center  and  surrounding 
voxels  were  compared  with  the  actual  values.  The  differ¬ 
ence  between  the  actual  pH  value  and  the  pH  value 
reported  by  the  simulation  program  expresses  the  error  due 
to  the  voxel  bleed  artifact,  as  shown  in  Fig.  9.  As  expected, 
when  the  center  voxel  is  surrounded  by  more  acidic 
voxels,  the  reported  pH  for  the  center  voxel  is  underesti¬ 
mated  and  vice  versa. 

The  error  in  the  reported  pH  for  the  center  voxel  over  the 
physiological  pH  range  (6. 8-7. 6  pH  units)  is  less  than  0.15 
pH  units  for  a  maximum  intensity  ratio  of  1:3.0,  and  less 
than  0.10  pH  units  for  an  average  intensity  ratio  of  1:1.5,  as 
shown  in  Fig.  9.  This  is  an  encouraging  result  for  the 
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FIG.  9.  The  error  due  to  the  voxel  bleed  artifact  for  the  center  voxel 
as  a  function  of  pH  and  signal  intensity  ratio  for  (a)  a  signal  intensity 
ratio  of  1 :3.0  and  (b)  1 :1 .5. 
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precision  of  the  proposed  technique.  Additionally,  it  must 
be  remembered  that  these  are  worst-case  scenarios.  In 
practice,  it  is  likely  that  the  pH  difference  between  adja¬ 
cent  voxels  is  less  than  0.4  pH  units,  thus  further  reducing 
the  error  due  to  the  voxel  bleed  artifact. 

CONCLUSIONS  AND  FUTURE  DIRECTIONS 

Using  phantom  and  in  vivo  data,  we  have  demonstrated 
the  feasibility  of  imaging  extracellular  pH  using  ^H  MRSI. 
The  possible  sources  of  error  associated  with  the  technique 
were  analyzed  in  detail.  It  was  found  that  the  COG  method 
leads  to  an  error  of  <0.1  pH  units  for  broader  lines  and 
reasonable  SNR.  For  narrower  lines,  the  poor  spectral 
resolution  increasingly  reduces  this  accuracy.  From  simu¬ 
lations,  it  was  found  that  artifacts  due  to  voxel  bleed  in 
MRSI  will  introduce  a  maximum  error  of  0.1  pH  units  in 
the  reported  pH  values.  Finally,  it  was  observed  that  the  T2 
of  the  proposed  pHe  indicator  is  pH  dependent.  This 
dependence,  its  cause,  and  the  effect  on  reported  pH  values 
are  currently  under  investigation. 
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Nm23-Transfected  MDA-MB-435  Human  Breast 
Carcinoma  Cells  Form  Tumors  With  Altered  Phospholipid 
Metaholism  and  pH:  A  Nuclear 
Magnetic  Resonance  Study  In  Vivo  and  In  Vitro 

Zaver  M.  Bhujwalla,^*  Eric  O.  Aboagye,^  Robert  J.  Gillies, ^  V.P.  Chacko,^ 

Charmaine  E.  Mendola,^  and  Joseph  M.  Backer® 


Nm23  genes  are  involved  in  the  control  of  the  metastatic 
potential  of  breast  carcinoma  cells.  To  understand  the  impact  of 
nm23  genes  on  tumor  physiology  and  metabolism,  a  nuclear 
magnetic  resonance  (NMR)  spectroscopic  study  was  performed 
on  tumors  formed  in  the  mammary  fat  pad  of  severe  combined 
immunodeficiency  mice  by  MDA-MB-435  human  breast  carci¬ 
noma  cells  transfected  with  cDNA  encoding  wild  type  nm23-H1 
and  nm23-H2  proteins.  Tumors  formed  by  MDA-MB-435  cells 
transfected  with  vector  alone  were  used  as  controls.  All  trans¬ 
gene  tumors  exhibited  significantly  higher  levels  of  phosphodi- 
ester  (PDE)  compounds  relative  to  phosphomonoester  (PME) 
compounds  in  vivo  compared  with  control  tumors.  Similar 
differences  in  PDE  and  PME  also  were  observed  for  spectra 
obtained  from  cells  growing  in  culture.  Intracellular  pH  was 
significantly  lower  and  extracellular  pH  was  significantly  higher 
for  transgene  tumors  compared  with  control  tumors.  Histologic 
analysis  of  lung  sections  confirmed  reductions  in  incidence, 
number,  and  size  of  metastatic  nodules  for  animals  bearing 
transgene  tumors.  These  results  suggest  that  nm23  genes  may 
affect  suppression  of  metastasis  through  phospholipid-medi¬ 
ated  signaling  and  cellular  pH  regulation.  Magn  Reson  Med 
41:897-903, 1999.  ©  1 999  Wiley-Liss,  Inc. 

Keywords:  human  breast  carcinoma  metastasis;  nonmetastatic 
23  transfection;  phospholipid  metabolism  and  pH;  ^^P  nuclear 
magnetic  resonance  spectroscopy 

The  ability  of  solid  tumors  to  metastasize  and  establish 
colonies  at  distant  sites  is  one  of  the  most  life-threatening 
aspects  of  cancer.  Despite  continuing  advances  in  the 
molecular  characterization  of  events  promoting  metastasis, 
little  impact  has  been  made  on  therapy  or  survival  for 
patients  with  advanced  metastatic  tumors  (1).  This  is  due 
partly  to  the  lack  of  identifiable  targets  against  which  to 
design  antineoplastic  agents  to  control  the  metastatic  spread 
of  cancer.  Multinuclear  magnetic  resonance  (multi-NMR) 
methods  have  a  unique  role  to  play  in  answering  this 
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challenge  by  providing  an  understanding  of  the  biochemi¬ 
cal  and  physiologic  mechanisms  involved  in  invasion  and 
metastasis.  Such  an  understanding  can  identify  rational 
targets  for  therapy. 

Recently  it  was  shown  that  the  nm23  (nonmetastatic) 
gene  is  related  to  suppression  of  metastasis;  the  metastasis 
suppression  function  of  the  nm23  gene  was  proposed  on 
the  basis  of  correlation  and  transfection  studies  in  murine 
and  human  systems  (2,3).  Two  highly  homologous  and 
evolutionary  conserved  nm23  genes,  nm23-l  and  nm23-2, 
have  been  identified  in  rodents  (4,5),  and  two  nm23  genes 
(nm23-Hl  and  nm23-H2)  have  been  identified  in  humans 
(6,7).  The  two  murine  nm23  genomic  DNAs  have  been 
cloned  and  sequenced  (8,9).  The  human  nm23-Hl  and 
nm23-H2  genes  have  been  localized  to  chromosome  17q21 
(10,11).  These  genes  encode  17-kDa  proteins  identified  as 
nucleoside  diphosphate  kinase  A  (NDPK  A)  and  NDPK  B, 
which  form  homomers  and  heteromers.  In  addition,  NDPK 
B  displays  an  increasing  list  of  other  activities  that  appar¬ 
ently  are  unrelated  to  its  catal3dic  functions  (4).  However, 
the  cellular  mechanisms  by  which  the  nm23  protein 
suppresses  metastatic  phenotypic  expression  is  still  un¬ 
known.  In  the  current  study,  we  have  used  ^^P  NMR 
spectroscopy  to  study  metabolic  and  physiologic  character¬ 
istics  of  tumors  induced  in  severe  combined  immunodefi¬ 
ciency  (SCID)  mice  by  MDA-MB-435  human  breast  carci¬ 
noma  cells  transfected  with  wild  type  cDNA  of  nm23-Hl 
and  nm23-H2  and  demonstrated  that  nm23  transfection 
alters  phospholipid  metabolism  and  pH  in  these  breast 
tumors.  Histologic  analysis  of  lung  sections  from  tumor¬ 
bearing  animals  was  performed  to  confirm  decreases  in 
metastatic  behavior  of  these  breast  tumors  following  trans¬ 
fection  with  nm23-Hl  and  nm23-H2.  These  studies  have 
provided  further  understanding  of  the  cellular  functions  of 
nm23  and  of  the  mechanisms  of  action  of  nm23-l  and 
nm23-2  genes  and  their  role  in  metastatic  dissemination  of 
tumor  cells. 

MATERIALS  AND  METHODS 

Coding  sequences  of  normal  nm23-Hl  and  nm23-H2  pro¬ 
teins  were  cloned  into  the  eucaryotic  expression  vector 
pPalPstNeo  under  control  of  a  constitutive  HCMV  pro¬ 
moter  (11,12).  The  vector  contains  an  neo  resistance  gene 
under  control  of  a  simian  virus  24  (SV40)  promoter. 
MDA-MB-435  breast  carcinoma  cells  were  transfected  with 
nm23  constructs  by  using  a  Lipofectin  kit  (BRL-Life  Tech¬ 
nologies,  Inc.,  Gaithersburg,  MD),  and  selection  of  trans- 
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FIG.  1.  Tumor  volumes  following  inoculation  of  identical  cell  num¬ 
bers  (10®  cells  in  0.05  ml  of  Hank’s  balanced  salt  solution)  in  the 
mammary  fat  pad  for  the  groups  of  animals  used  in  the  study.  Bars 
represent  ±  1  standard  error  of  the  mean  (S.E.M.). 

fected  clones  was  performed  in  the  presence  800  p-g/ml 
G418.  Nm23-transfected  pooled  clones  of  MDA-MB-435 
were  then  transfected  with  plZsp-Pgluc  (puro),  a  mamma¬ 
lian  expression  vector  containing  bacterial  p-glucuroni- 
dase  under  control  of  a  constitutive  HCMV  promoter  and 
puromycin-resistance  gene  under  control  of  an  SV40  pro¬ 
moter  (gift  from  Dr.  T.  Jones,  Lederle  Laboratories).  Selec¬ 
tion  of  clones  expressing  bacterial  p-glucuronidase  was 
performed  in  the  presence  of  0.375  ixg/ml  puromycin  and 
800  |xg/ml  G418.  Pooled  clones  of  double-transfected  cells 
named  MDA-MB-435-Vp,  MDA-MB-453-ip,  and  MDA-MB- 
435-2p  (for  vector/p-glucuronidase,  nm23-Hl/p-glucuroni- 
dase,  and  nm23-H2/p-glucuronidase  transfections,  respec¬ 
tively)  were  maintained  in  the  presence  of  0.375  ixg/ml 
puromycin  and  200  fjig/ml  G418  in  Dulbecco’s  minimum 
essential  medium  (Sigma  Ltd.,  St.  Louis,  MO)  containing 
10%  fetal  bovine  serum  (Sigma  Ltd.).  Presence  of  trans- 
genes  was  confirmed  by  polymerase  chain  reaction  (PGR) 
analysis.  Western  blot  analysis,  and  staining  with  5-bromo- 
4-chloro-3-indol  1  glucuronide  (X-glu).  Cell  doubling  times 
were  determined  by  counting  cells  plated  in  10  mm  Petri 
dishes  (three  plates  per  cell  line)  using  a  hemocytomter. 
Cells  were  counted  48  hr  after  plating  5X10^  cells  per  Petri 
dish.  Cell  doubling  times  were  approximately  1.5-2. 5 
days,  with  no  significant  differences  between  the  cell  lines. 
Mean  values  ±  1  standard  error  of  the  mean  (S.E.M.)  for  the 
three  lines  were  as  follows:  2,3  days  ±  0,4  (MDA-MB-435- 
Vp),  1.5  days  ±  0.1  (MDA-MB-435-ip),  and  2.4  days  ±  0.6 
(MDA-MB-435-2P). 

MDA-MB-435-Vp,  MDA-MB-435-lp,  and  MDA-MB- 
435-2p  cells  were  inoculated  in  the  upper  left  thoracic 
mammary  fat  pad  of  SCID  mice,  and  10®  cells  were 
inoculated  in  0.05  ml  of  Hank’s  balanced  salt  solution 
(Sigma  Ltd.).  The  experimental  protocol  was  approved  by 
the  Institutional  Animal  Care  and  Use  Committee.  Mice 
were  anesthetized  with  ketamine  (50  mg/kg;  Aveco  Ltd.) 
and  acepromazine  (5  mg/kg;  Aveco  Ltd.).  Tumor  volumes, 
which  were  measured  just  prior  to  performing  the  ®^P  NMR 
spectroscopic  studies,  were  calculated  from  caliper  mea¬ 
surements  of  tumor  axes  (a,b,c)  by  using  the  equation  for  an 
elliptical  volume  ('TT/6)abc.  Volumes  used  in  this  study 
were  of  the  order  of  300  mm^. 


NMR  spectroscopic  studies  were  performed  on  a  GE 
CSI  4.7  T  instrument  (General  Electric,  Fremont,  CA) 
equipped  with  shielded  gradients.  Spectra  were  obtained 
with  home-built  solenoidal  coils  fitted  around  the  tumor. 
Because  the  coil  design  allowed  the  probe  to  be  tuned  to 
the  proton  frequency  for  shimming,  proton  images  were 
acquired  with  the  coil  to  ascertain  that  only  the  tumor  (and 
skin)  was  in  the  sensitive  volume  of  the  coil.  Animal  body 
temperature  was  maintained  at  37°C  by  heat  generated 
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FIG.  2.  Representative  fully  relaxed  P  nuclear  magnetic  resonance 
(NMR)  spectra  obtained  from  an  MDA-MB-435-Vp  tumor  (a)  and 
from  an  MDA-MB-435-ip  tumor  (b).  Peak  assignments  are  1) 
3-aminopropylphosphonate  (3-APP);  2)  phosphomonoester  (PME); 
3)  inorganic  phosphate  (Pi);  4)  phosphodiester(PDE);  5)  phosphocre- 
atine  (PCr);  6)  7-nucleoside  triphosphate  (NTP);  7)  a-NTP  (set  to 
-10  ppm);  and  8)  p-NTR 
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from  a  pad  circulating  with  warm  water.  For  the  NMR 
studies,  mice  were  injected  intraperitoneally  with  a  solu- 

^  tion  of  the  extracellular  pH  marker  3-aminopropylphospho- 
nate  (3-APP;  Sigma  Ltd.)  administered  in  a  volume  of  0.2 
ml  saline  (480  mg/kg)  following  anesthetization.  Fully 
relaxed  ^^P  NMR  spectra  were  obtained  by  using  a  45°  flip 
angle  with  64  scans,  a  repetition  of  5  sec,  and  a  sweep 
width  (SW)  of  10000  Hz.  Parameters  were  determined  from 
two  spectra  obtained  per  tumor.  NMR  examinations  were 
completed  within  20  min.  Extracellular  pH  (pHe)  was 
obtained  from  the  chemical  shift  of  3-APP  (13),  and 
intracellular  pH  (pHi)  was  obtained  from  the  chemical  shift 
of  inorganic  phosphate  (Pi)  (14)  from  the  endogenous 
reference  a-nucleoside  triphosphate  (NTP)  set  to  -7.57 
ppm.  pHe  was  calculated  from  the  relationship  pH  = 
6.91  +  log  [(8  3-APP  21.11)7(24.30  ^s-APp)]*  snd  pHi  was 

calculated  from  the  relationship  pH  =  6.66  +  log  [(8pi  - 
0.65)7(3.11  —  8pi)].  Data  sets  were  processed  by  using  an 
exponential  line-broadening  factor  of  22  Hz.  Peak  areas 
were  determined  in  the  time  domain  by  using  an  in-house, 
nonlinear,  least-squares,  curve-fitting  routine  for  MR  data 
analysis  written  by  Dr.  D.C.  Shungu. 

High-resolution  spectra  of  perchloric  acid  (PGA)  extracts 
from  tumors  and  cell  lines  were  obtained  to  resolve  the 
components  of  the  peak  in  the  phosphodiester  (PDE)  and 
the  phosphomonoester  (PME)  regions.  Mice  were  anesthe¬ 
tized,  and  tumors  were  excised  and  immediately  freeze 
clamped.  Neutralized  PGA  extracts  of  tumors  were  lyophi- 
lized  and  resuspended  in  D2O.  High-resolution  ^^P  NMR 
spectra  of  tumor  extracts  were  acquired  at  11.7  T  (MSL-500 
spectrometer;  Bruker)  with  a  10-mm  high-resolution  prohe. 
Spectra  were  acquired  with  a  45°  flip  angle,  an  SW  of  8000 
Hz,  a  repetition  time  of  5  sec,  a  block  size  of  4  K,  and  3600 
scans. 

PGA  extracts  were  obtained  from  equal  numbers  of  cells 
for  all  of  the  lines  used  in  the  study.  Gell  volumes  of  the 
lines  were  identical.  Gells  from  approximately  eight  flasks 
with  similar  confluency  for  each  cell  line  were  trypsinized, 
and  the  action  of  trypsin  was  blocked  with  ice-cold  growth 
medium.  Gells  were  fed  approximately  3  hr  prior  to 
trypsinization.  Gells  were  washed  twice  with  cold  0.9% 
NaGl  solution  and  extracted  with  ice-cold  8%  (volume7 
volume)  PGA.  The  supernatant  was  neutralized  (with  3  M 
K2GO371  M  KOH),  lyophilized,  and  resuspended  in  D2O. 
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Phosphorus  spectra  of  the  extracts  were  obtained  at  11.7  T 
with  a  10-mm  high-resolution  prohe.  Spectral  acquisition 
parameters  were  45°  flip  angle,  SW  =  8000  Hz,  5  sec 
repetition  time,  4  K  block  size,  and  scans  =  32000-128000. 

Lungs  from  tumor-bearing  animals  were  excised  at  the 
end  of  the  NMR  experiments  and  fixed  in  10%  buffered 
formalin.  Three  5-|uLm-thick,  paraffin-embedded  sections 
were  obtained  from  each  pair  of  lungs  and  stained  with 
hematoxylin  and  eosin.  Lung  sections  were  examined 
under  an  optical  microscope  and  evaluated  for  incidence, 
number,  and  size  of  metastatic  nodules. 

Statistical  analysis  of  the  data  was  performed  by  using 
StatView  II  software  (version  1.04;  Abacus  Goncepts,  Inc., 
Berkeley,  GA).  One  factorial  analysis  of  variance  (ANOVA) 
was  used  to  evaluate  the  statistical  significance  of  the  in 
vivo  data. 


RESULTS 

All  of  the  cell  lines  were  tumorigenic  in  SGID  mice,  with  a 
latent  period  of  5-6  weeks.  Growth  rates  for  control  and 
transgene  tumors  were  similar,  with  a  volume- doubling 
time  of  10-14  days.  Tumor  volumes  for  groups  of  animals 
inoculated  at  the  same  time  are  shown  in  Figure  1.  These 
data  demonstrate  the  absence  of  any  significant  differences 
in  growth  rate  or  “silent  interval”  following  inoculation  of 
the  different  cell  lines  in  the  mammary  fat  pad. 

The  most  striking  difference  between  ^^P  NMR  spectra  of 
control  (Fig.  2a)  and  transgene  tumors  (Fig.  2b)  was  a 
marked  increase  of  the  peak  in  the  PDE  region  relative  to 
the  PME  region  in  spectra  of  transgene  tumors.  The  compo¬ 
nents  of  the  PME  region  were  identified  as  the  membrane 
precursors  phosphocholine  (PG)  and  phosphoethanol- 
amine  (PE),  and  those  of  the  PDE  region  were  identified  as 
the  membrane  breakdown  products  glycerophosphocho- 
line  (GPG)  and  glycerophosphoethanolamine  (GPE).  These 
components  were  identified  from  the  high-resolution  ^^P 
NMR  spectra  of  tumor  extracts  by  using  the  spectral 
assignments  of  Evanochko  et  al.  (15).  Elevation  of  the  PDE 
peak  relative  to  the  PME  peak  was  observed  consistently 
for  all  transgene  tumors  but  not  for  MDA-MB-435-Vp 
control  tumors.  An  example  of  a  high-resolution  spectrum 
obtained  from  PGA  extracts  of  MDA-MB-435-lp  tumors 

MDA-MB-435-1  p  tumor  extract 


1 


FIG.  3.  High-resolutlon  NMR  spectrum  from 
MDA-MB-435-1  tumor  extract  (obtained  from 
three  pooled  tumors).  Spectral  acquisition  pa¬ 
rameters  for  in  vivo  tumors  and  extract  are 
detailed  in  Materials  and  Methods.  GPC,  glycero- 
phosphocholine;  GPE,  glycerophosphoethanol¬ 
amine;  PC,  phosphocholine;  PE,  phosphoetha- 
nolamine. 
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Table  1 

Nuclear  Magnetic  Resonance  Parameters  and  Tumor  Volumes  Obtained  From  In  Vivo  Tumors  Derived 
From  Control  and  Transgene  MDA-MD-435  Cells 


Solid  tumor  type 

Transfection  and  nm23  levels 

Tumor  volume 
(mm^) 

PDE/PME  = 
[GPE  +  GPC] 
[PE  +  PC] 

pHi 

pHe 

Control  MDA-MB-435-Vp  (n  =  7) 

Vector/p-glucuronidase  transfection 
(low  nm23-H1,  H2) 

324  ±  42 

0.60  ±  0.05 

7.37  ±  0.07 

6.8  ±  0.11 

Transgene  MDA-MB-435-1  p  (n  =  6) 

nm23-H1/p-glucuronidase  transfection 
(high  nm23-H1) 

311  ±  45 

1.45  ±  0.24* 

7.16  ±0.05*  7.17  ±0.1* 

Transgene  MDA-MB-435-2P  (n  =  6) 

nm23-H2/p-glucuronidase  transfection 
(high  nm23-H2) 

255  ±  19 

1.36  ±  0.14* 

7.15  +  0.06*  7.09  ±0.08* 

*Ninety-five  percent  confidence  limit  (P  <  0.05),  analysis  of  variance,  Fisher’s  protected  least  significant  difference  test.  Values  represent 
mean  ±  1  standard  error  of  the  mean;  n  represents  number  of  animals  for  each  group;  PDE,  phosphodiester;  PME,  phosphomonoester;  GPE, 
glycerophosphoethanolamine;  GPC,  glycerophosphocholine;  PE,  phosphoethanolamine;  PC,  phosphocholine;  pHi,  intracellular  pH;  pHe, 
extracellular  pH. 


with  the  corresponding  peak  assignments  is  shown  in 
Figure  3. 

The  in  vivo  results  for  all  of  the  animals  in  the  study  are 
summarized  in  Table  1  and  show  that  the  PDE /PME  ratio 
was  significantly  higher  for  the  transgene  tumors  com- 
pared  with  control  tumors.  Significant  differences  in  pHi 
and  pHe  also  were  detected  for  the  transgene  tumors.  pHi 
was  significantly  lower,  whereas  pHe  was  significantly 
higher  for  nm23-Hl-  and  nm23-H2-transfected  tumors 
compared  with  those  derived  from  cells  transfected  with 
vector  only  (Table  1).  No  significant  differences  in  NTP/Pi 
were  detected  between  the  cell  lines. 

^^P  NMR  high-resolution  spectra  of  isolated  cell  extracts 
obtained  from  cells  growing  in  tissue  culture  flasks  for 
each  of  the  transfected  cell  lines  and  the  control  cell  line 
are  shown  in  Figure  4.  PDE/PME  ratios  obtained  from  cell 
extracts  are  summarized  in  Table  2.  Observations  made  in 
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FIG.  4.  HIgh-resolutlon  ^^P  NMR  spectra  of  cells  from  MDA-MB- 
435-Vp  (a),  MDA-MB-435-1  p  (b),  and  MDA-MB-435-2p  (c)  cell  lines 
maintained  in  tissue  cuiture  flasks.  Spectra  are  expanded  to  focus  on 
the  phospholipid  region.  Spectral  acquisition  parameters  are  detailed 
in  Materials  and  Methods. 


vivo  also  were  apparent  in  the  spectra  of  cell  extracts, 
suggesting  that  differences  in  PME  and  PDE  peaks  for  the 
transgene  tumors  are  due  to  intrinsic  cellular  properties 
arising  from  transfection  of  cells  with  nm23  rather  than  in 
vivo  physiologic  effects  related  to  tumor  vascularization  or 
the  fraction  of  necrosis. 

For  the  evaluation  of  metastasis,  we  analyzed  the  lungs 
of  five  animals  bearing  MDA-MB-435-Vp  tumors  and  six 
animals  with  MDA-MB-435-1  p  tumors.  Because  the  two 
nm23  genes  (nm23-Hl  and  nm23-H2)  are  88%  to  90% 
identical  in  their  amino  acid  sequences,  respectively,  we 
analyzed  the  lungs  of  two  randomly  picked  animals  with 
MDA-MB-435-2P  tumors  that  showed  no  evidence  of  lung 
metastasis.  Data  from  micrcoscopic  analysis  of  the  lung 
specimens  obtained  from  tumor  bearing  mice  are  presented 
in  Table  3.  Lungs  obtained  from  animals  bearing  transgene 
tumors  showed  a  reduction  in  the  incidence  as  well  as  the 
number  and  size  of  metastatic  nodules  compared  with 
control  tumors  (Fig.  5). 

DISCUSSION 

The  metastatic  cascade  is  a  complex  phenomenon.  To 
establish  metastatic  colonies  at  a  site  distant  from  the 
primary  tumor,  a  cancer  cell  should  pass  successfully 
through  the  following  stages:  invasion,  intravasation,  ar¬ 
rest  of  cancer  cells,  extravasation,  and,  finally,  neovascular- 


Table  2 

Phosphodiester/Phosphomonoester  Ratios  of  Control 
and  Transgene  Cells 


Tumor  cell  type 

Transfection  and 
nm23  levels 

(PDE/PME)  = 
[GPE  +  GPCF 
[PE  +  PC] 

Vector/p-glucuronidase 

Control  MDA-MB-435-Vp 

transfection 

(lownm23-H1,2) 

0.53 

Transgene  MDA-MB- 
435-1  p 

nm23-H1/p-glucuroni- 
dase  transfection 
(high  nm23-H1) 

1.96 

Transgene  MDA-MB- 
435-2P 

nm23-H2/p-glucuroni- 
dase  transfection 
(high  nm23-H2) 

1.7 

^Values  obtained  from  eight  flasks  per  cell  line.  Cells  were  obtained 
from  a  total  of  eight  T-1 50  flasks  for  each  cell  line. 


Noninvasive  In  Vivo  NMR  Following  Nm23  Transfection 


901 


c  d 


FIG.  5.  Photomicrographs  demonstrating  the  marked  differences  in  metastatic  nodules  in  lung  sections  obtained  from  mice  bearing 
l\/IDA-MB-435-V(i  and  MDA-MB-435-1  (3  tumors.  a,b:  Low-  and  high-power  photomicrographs,  respectively,  from  lung  sections  obtained  from 
a  mouse  with  an  MDA-MB-435-Vp  tumor.  c,d:  Low-  and  high-power  photomicrographs,  respectively,  from  lung  sections  obtained  from  a 
mouse  with  an  MDA-MB-435-1  p  tumor.  Metastatic  nodules  are  marked  by  arrows,  and  n  denotes  a  region  of  necrosis  in  the  metastatic  nodule 
shown  in  b.  No  nodules  were  observed  in  lung  sections  that  were  examined  for  mice  bearing  MDA-MB-435-2P  tumors. 


ization.  Attenuation  of  metastasis  may  occur  at  any  of  these 
stages.  Histologic  analyses  of  lung  sections  from  tumor- 
bearing  animals  demonstrated  that  control  tumors  and 
transgene  tumors  formed  by  derivatives  of  MDA-MB-435 
human  breast  carcinoma  cells  transfected  with  nm23  con¬ 
structs  showed  different  metastatic  potential.  NMR 
spectra  of  primary  tumors  in  SCID  mice  revealed  signifi¬ 
cant  differences  in  the  phospholipid  composition,  pHi, 
and  pHe  of  control  and  transgene  tumors.  Because  trans¬ 
gene  tumors  were  formed  by  pooled,  transfected  cells, 
these  alterations  were  not  due  to  properties  of  individual 
clones.  These  differences  in  NMR  spectra  have  pro¬ 
vided  an  insight  into  how  the  loss  of  the  nm23  gene 
induces  metabolic  and  physiologic  alterations  that  may 
facilitate  metastatic  dissemination. 

Transgene  tumors  in  vivo  formed  by  wild  type  forms  of 
nm23-Hl  and  nm23-H2  exhibited  a  significantly  higher 
amount  of  PDE  relative  to  PME  levels  compared  with 
control  tumors.  These  differences  in  phospholipid  metabo- 


Table  3 

Incidence,  Number,  and  Size  of  Metastatic  Nodules  Detected  in 
Histologic  Sections  Obtained  from  Lungs  of  Tumor-Bearing  Animals 


Tumor  type 

Incidence  of 
metastatic 
nodules^ 

Number  of 
metastatic 
nodules 
averaged  for 
all  animals 

Range  of 
sizes  of 
metastatic 
nodules 

Control  MDA-MB- 
435-Vp 

5/5 

20 

50-3000  pm^ 

Transgene  MDA-MB- 
435-1  p 

4/6 

2 

50-700  pm'^ 

Transgene  MDA-MB- 
435-2P 

0/2 

0 

— 

umber  of  animals  with  lungs  containing  metastatic  nodules/ 
number  of  animals  with  lungs  examined  (three  histologic  sections 
were  analyzed  for  each  lung  specimen). 

^Mean  nodule  diameter  averaged  over  all  the  sections  was  800  pm. 
^Only  one  nodule  was  700  pm  in  diameter.  The  remaining  nodules 
were  less  than  100  pm  in  diameter. 
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lism  also  were  detected  in  cultured  cell  extracts.  The 
higher  PDE  levels,  including  an  elevation  of  GPC,  detected 
for  the  transgene  tumors  therefore  appeared  to  be  due  to 
cellular  characteristics  rather  than  to  in  vivo  physiologic 
characteristics,  such  as  vascularization  or  necrosis.  GPC  is 
formed  by  degradation  of  membrane  phospholipids  by  the 
enzyme  phospholipase  Aa  (PLAa),  and  this  mechanism  of 
PLAa  on  membrane  phospholipids  provides  precursors  for 
eicosanoids  when  the  cleaved  fatty  acid  is  arachidonic  acid 
(16).  Arachidonic  acid-derived  eicosanoids  include  impor¬ 
tant  lipid  mediators  in  signal  transduction  and  tumor 
promotion,  such  as  prostaglandins,  thromboxanes,  and 
leukotrienes  (17).  These  oxygenated  fatty  acids  can  have 
diverse  effects  on  cancer  cell  motility,  invasion,  and  metas¬ 
tasis  (17-20),  and  the  changes  in  phospholipid  metabolism 
detected  following  nm23  transfection  may  reflect  the  in¬ 
volvement  of  the  nm23  gene  product  in  phospholipid- 
mediated  signaling  mechanisms  regulating  cell  motility. 
This  possibility  is  consistent  with  the  observation  by 
Kantor  et  al.  (21)  that  murine  melanoma  and  human  breast 
carcinoma  cells  that  are  transfected  stably  with  nm23-Hl 
do  not  migrate  in  response  to  chemoattractants. 

The  second  significant  metabolic  effect  was  the  differ¬ 
ence  in  both  pHi  and  pHe  between  the  control  and 
transgene  tumors.  Transgene  tumors  formed  by  MDA-MB- 
435— Ip  and  MDA-MB-435— 2p  cells  showed  significantly 
lower  steady-state  pHi  values  compared  with  control  tu¬ 
mors.  Studies  of  cells  within  a  bioreactor  system  have 
shown  that  low  pHi  can  result  in  an  increase  of  GPC  and  a 
decrease  of  PC,  which  is  most  likely  due  to  inhibition  of 
GPC/GPE  phosphodiesterase  activity  (22).  Inhibition  of 
GPC/GPE  phosphodiesterase  activity  may  explain  in  part 
the  differences  in  PDE  and  PME  levels  for  the  control  and 
transgene  cancer  cells  (23,24).  In  the  bioreactor  system, 
however,  inhibition  of  GPC/GPE  phosphodiesterase  activ¬ 
ity  occurred  only  for  pHi  lower  than  7.  Because  pHi  values 
in  all  of  the  tumors  in  our  study  were  greater  than  7, 
inhibition  of  GPC/GPE  phosphodiesterase  activity  cannot 
explain  entirely  the  differences  in  PDE  and  PME  observed. 
The  lower  pHe  values  of  the  control  tumors  and  the  higher 
pHe  values  of  transgene  tumors  formed  by  MDA-MB- 
435— Ip  and  MDA-MB-435— 2p  cells  are  particularly  interest¬ 
ing,  because  low  pHe  may  enhance  the  invasive  behavior 
of  human  breast  carcinoma  cells.  An  acidic  pericellular 
pH,  for  instance,  was  found  to  increase  the  secretion  of  the 
active  form  of  the  lysosomal  protease  cathepsin  B  over  time 
for  human  breast  carcinoma  cells  (25).  Human  melanoma 
cells  have  been  observed  to  secrete  a  higher  level  of  90-kDa 
gelatinase  (a  Type  IV  collagenase)  at  a  pHe  of  6.8  compared 
with  pH  7.3  (26).  This  acid-induced  secretion  of  gelatinase 
was  blocked  by  cycloheximide,  indicating  that  the  enzyme 
induction  was  due  to  de  novo  synthesis.  Thus,  pH-related 
mechanisms  may  also  play  a  role  in  facilitating  invasion 
and  metastasis  in  vivo  for  breast  carcinoma  cells  that  lack 
the  nm23  gene. 

Although  the  physiologic  and  metabolic  effects  detected 
in  our  study  may  alter  invasive  and  metastatic  behavior 
independently,  it  is  also  possible  that  they  may  be  synergis¬ 
tic  or  necessary  but  not  sufficient  without  the  other.  The 
results  obtained  here  may  open  new  opportunities  for 
diagnosis/prognosis  of  metastatic  dissemination  as  well  as 
potential  targeting  of  phospholipid  metabolism  and  pH 


regulation  for  antimetastatic  therapeutics.  This  is  the  first 
in  vivo  observation  that  links  activity  of  a  putative  metasta¬ 
sis  suppressor  nm23  gene  to  metabolic  processes.  The 
results  also  demonstrate  the  potential  of  noninvasive  NMR 
to  detect  forms  of  gene  therapy  for  suppression  of  metasta¬ 
sis  that  may  involve  transfection  of  cells  with  nm23. 
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Abstract 

The  monitoring  of  antibody-directed  enzyme-prodrug  therapies  re¬ 
quires  evaluation  of  drug  activation  within  the  tissues  of  interest.  We  have 
demonstrated  the  feasibility  of  noninvasive  magnetic  resonance  spectros¬ 
copy  and  spectroscopic  imaging  (chemical  shift  imaging)  to  detect  activa¬ 
tion  of  the  prodrug  5-fluorocytosine  (5-FCyt)  to  the  cytotoxic  species 
5-fluorouraci]  (5-FU)  by  monoclonal  antibody-cytosine  deaminase  (CD) 
conjugates.  In  vitro,  L6-CD  but  not  1F5-CD  selectively  metabolized 
5-FCyt  to  S-FU  on  H2981  human  lung  adenocarcinoma  cells  because  of 
the  presence  and  absence  of  cel!  surface  L6  and  CD20  antigens,  respec¬ 
tively.  After  pretreatment  of  H2981  tumor-bearing  mice  with  L6-CD,  in 
vivo  metabolism  of  5-FCyt  to  5-FU  within  the  tumors  was  detected  by  *^F 
magnetic  resonance  spectroscopy;  the  chemical  shift  separation  between 
5-FCyt  and  5-FU  resonances  was  —1.2  ppm.  5-FU  levels  were  50-100%  of 
5-FCyt  levels  in  tumors  10-60  min  after  5-FCyt  administration.  Whole 
body  ^^F  chemical  shift  imaging  (6x6  mm  in-plane  resolution)  of 
tumor-bearing  mice  demonstrated  the  highest  signal  intensity  of  5-FU 
within  the  tumor  region.  This  study  supports  further  development  of 
noninvasive  magnetic  resonance  methods  for  preclinical  and  clinical  mon¬ 
itoring  of  CD  enzyme-prodrug  therapies. 

Introduction 

Present  therapeutic  modalities  for  treating  cancer  are  frequently 
ineffective  because  of  the  lack  of  specificity  of  the  agents  for  cancer 
cells.  The  ensuing  toxicity  to  normal  tissues  can  severely  limit  the 
dose  of  drug  administered.  A  promising  new  approach  for  increasing 
therapeutic  index  of  anticancer  therapy,  ADEPT, ^  is  based  on  the 
delivery  to  cancer  cells  of  a  nontoxic  antibody-enzyme  conjugate. 
After  a  suitable  time  to  allow  for  clearance  of  the  unbound  conjugate, 
a  nontoxic  prodrug  is  administered  that  is  selectively  converted  to  a 
toxic  compound.  This  paradigm  is  exemplified  by  systemic  adminis¬ 
tration  of  the  monoclonal  antibody-CD  conjugate,  L6-CD,  which 
renders  tumor  cells  susceptible  to  the  nontoxic  prodrug  5-FCyt  (1,2). 
The  L6-CD  contains  a  yeast  enzyme,  CD,  which  converts  5-FCyt  to  a 
potent  cytotoxic  agent,  5-FU.  The  specificity  of  this  therapeutic  ap¬ 
proach  to  cancer  cells  is  based  on  selective  binding  of  the  L6  mono¬ 
clonal  antibody  portion  of  the  conjugate  to  a  cell  surface  antigen 
strongly  expressed  on  most  human  carcinomas  (3,  4),  A  similar 
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approach  that  involves  delivery  of  a  gene  that  codes  for  bacterial  CD 
has  also  been  reported  by  other  investigators  (5-9). 

Unlike  traditional  chemotherapeutic  approaches  for  cancer  that  can 
be  monitored  by  analysis  of  the  drug  and  metabolite  levels  in  blood, 
clinical  monitoring  of  enzyme-prodrug  therapies  requires  evaluation 
of  drug  activation  within  the  tissue  of  interest,  tumor.  A  clinically 
applicable  noninvasive  method  for  monitoring  drug  activation  in 
target  tissue  will  allow  evaluation  of  drug  activation  profiles  and 
selectivity  for  tumors  versus  normal  tissues.  Such  methods  may  be 
more  acceptable  to  patients  than  multiple  biopsies,  leading  to  in¬ 
creased  compliance  to  the  therapeutic  protocols.  In  this  study,  we  have 
demonstrated  the  feasibility  of  noninvasive  MRS  and  spectroscopic 
imaging,  CSI,  for  monitoring  L6-CD/5-FCyt  therapy  in  vitro  and  in 
mice  bearing  human  lung  cancer  (H2981)  xenografts.  To  our  knowl¬ 
edge  this  constitutes  the  first  report  of  its  kind  and  is  also  applicable 
to  the  alternative  approach  involving  delivery  of  a  “suicide”  gene  that 
codes  for  CD  (5-9).  MRS  of  natural  abundant  fluorine,  ^^F  MRS,  has 
been  used  successfully  to  monitor  the  pharmacokinetics  and  metabo¬ 
lism  of  5-FU  in  both  experimental  animals  and  patients  (10-14).  ^^F 
MRS  offers  several  advantages  such  as  a  spin  of  1/2,  100%  natural 
abundance,  high  detection  sensitivity  (0.83%  that  of  protons),  and  low 
background  signal.  The  technique  is,  therefore,  suitable  for  monitor¬ 
ing  in  vivo  levels  of  fluorinated  compounds  and  their  metabolites, 
particularly  when  such  compounds  are  administered  at  high  doses 
(>100  mg/kg  body  weight)  as  is  the  case  for  5-FCyt  (1,  5).  In  this 
study,  we  have  shown  that  ^^F  MRS  and  CSI  can  be  used  to  monitor 
signal  intensities  of  5-FCyt  and  5-FU  in  tumors  as  well  as  their  whole 
body  biodistribution  in  mice.  Issues  relating  to  the  use  of  these 
techniques  to  evaluate  patient  tumors  are  also  discussed. 

Materials  and  Methods 

Materials.  Antibody  conjugates  L6-CD  and  1F5-CD  and  13B  anti-idio- 
typic  antibody  were  obtained  from  Bristol-Myers  Squibb  (Seattle,  WA).  L6 
(IgG2a)  binds  to  a  ganglioside  antigen  expressed  on  the  surface  of  human 
carcinomas  (4);  1F5  (IgG2a)  binds  to  the  CD20  antigen  expressed  on  both 
normal  and  neoplastic  B  cells  but  not  on  carcinomas  (2, 15);  13B  anti-idiotypic 
antibody  recognizes  and  binds  to  the  L6  portion  of  circulating  L6-CD  but  not 
of  bound  L6-CD  and  enhances  the  clearance  of  unbound  conjugate  (1, 16).  The 
human  lung  adenocarcinoma  cell  line  H2981  as  well  as  tumor  sections  derived 
from  this  cell  line  and  maintained  by  in  vivo  passage  were  also  obtained  from 
Bristol-Myers  Squibb  (3).  5-FU  and  5-FCyt  were  purchased  from  Sigma 
Chemical  Co.  (St.  Louis,  MO). 

In  Vitro  Studies.  H2981  ceils  were  cultured  in  Iscove’s  modified  Dulbec- 
co’s  medium  (Sigma)  supplemented  with  10%  fetal  bovine  serum,  200  units/ml 
penicillin,  and  0.1  mg/ml  streptomycin.  To  demonstrate  the  ability  of  MRS  to 
detect  selective  drug  activation  in  vitro,  subconfluent  H2981  cells  (~  2  X  10^ 
cells;  4  repeats)  were  incubated  for  30  min  with  growth  medium  alone  or  with 
medium  containing  10  /Jig/ml  L6-CD  or  1F5-CD.  Cells  were  then  washed  three 
times  with  PBS  and  incubated  for  1  h  with  medium  containing  5-FCyt.  Cells 
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Fig.  1.  Metabolism  of  5-FCyt  by  H298I  cells  pretreated  with  either 
L6-CD  (a  and  b)  or  1F5-CD  (c  and  d).  The  vertical  scale  is  propor¬ 
tional  to  the  signal  intensity,  which  was  normalized  to  display  com¬ 
parable  signal-to-noise  levels  for  cell  extracts  {a  and  c)  or  media  {b 
and  d).  Cells  were  incubated  with  conjugate  for  30  min  followed  by 
three  cycles  of  washing  and,  finally,  were  incubated  with  5-FCyt  for 
1  h.  Extraction  and  analysis  of  samples  were  performed  as  reported  in 
“Materials  and  Methods.” 
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Fig.  2.  Typical  in  vivo  *®F  spectrum  obtained  from  a  H2981  tumor  60  min  after  5-FCyt 
injection  (400  mg/kg  i.p.).  Mice  were  pretreated  with  L6-CD  and  13B  anti-idiotypic 
antibody  at  48  and  24  h,  respectively,  before  the  MRS  study.  Data  were  processed  with 
a  line  broadening  of  25  Hz. 

incubated  with  medium  alone  were  treated  with  either  5-FCyt  or  5-FU.  Media 
containing  5-FCyt  and  5-FU  were  collected  after  incubation  and  frozen  in 
liquid  nitrogen.  Cells  were  washed  three  times,  trypsinized,  and  frozen.  The 
cells  were  homogenized,  centrifuged  to  remove  cell  debris  (15,000  rpm;  4°C; 
15  min),  and  lyophilized.  Both  incubation  media  (0,5  ml)  and  lyophilized  cell 
extracts  dissolved  in  0.5  ml  of  distilled  water  were  analyzed  for  fluorine- 
related  compounds  on  a  Bruker  500-MHz  NMR  spectrometer.  Spectra  were 
acquired  using  a  one-pulse  sequence  with  a  flip  angle  of  45°,  5  s  repetition 
time,  16,384-block  size,  22,000-Hz  sweep  width,  and  16/64  scans. 

Tumors.  H2981  tumor  sections  of  approximately  32  mm^  were  implanted 
s.c.  in  the  right  flank  of  female  B  ALB/cAnNCr-nu  nude  mice  (National  Cancer 
Institute,  Frederick,  MD).  Mice  were  given  laboratory  chow  and  water  ad 
libitum.  Tumor  sizes  of  500-600  mm^  were  used  for  the  MR  experiments.  The 
experimental  protocol  was  approved  by  the  Institutional  Animal  Care  and  Use 
Committee.  Mice  were  anesthetized  by  i.p,  injection  of  a  mixture  of  ketamine 
(25  mg/kg;  Aveco  Ltd.,  Fort  Dodge,  lA),  acepromazine  (2.5  mg/kg;  Aveco), 
and  0.9%  NaCl  solution  (1:1:2  by  volume). 

In  Vivo  Studies.  In  vivo  experiments  were  performed  to  demonstrate:  (o) 
activation  of  5-FCyt  within  tumors  by  MRS;  and  (b)  whole  body  distribution 
of  5-FCyt  and  5-FU  by  CSI,  H2981  tumor-bearing  mice  were  treated  i.v.  with 
L6-CD  (300  jLig/25  g  mouse)  48  h  before  the  experiment.  To  eliminate 
unbound  circulating  conjugate,  13B  anti-idiotypic  antibody  (200  /jig/25  g 
mouse)  was  injected  i.p.  in  the  same  mice  24  h  before  the  experiment.  Control 
animals  did  not  receive  antibody  treatment.  MRS  experiments  were  performed 
on  GE  CSI  4.7T  and  9.4T  NMR  spectrometers  with  a  two-turn  solenoidal  coil 
tunable  to  or  ^^F  frequency.  Mice  were  anesthetized  and  a  catheter  to 
deliver  5-FCyt  was  placed  i.p.  before  commencing  the  MRS  studies.  Body 
temperature  was  maintained  within  the  magnet  by  a  thermostat-regulated 


heating  pad.  For  MRS  studies,  5-FCyt  was  injected  as  a  single  bolus  dose  of 
400  mg/kg  after  acquisition  of  baseline  spectra.  To  evaluate  drug  activation  in 
tumors,  serial  ‘^F  spectra  were  acquired  every  5  min  for  1  h  using  a  one-pulse 
sequence  (60°  flip  angle,  3-s  repetition  time,  1,024-block  size,  10,000-Hz 
sweep  width  and  32/64  scans). 
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Fig.  3.  Profile  of  5-FCyt  and  5-FU  in  H2981  tumors  following  treatment  with  400 
mg/kg  5-FCyt  (i.p.).  o,  control  mice  not  receiving  conjugate  (n  =  3),  in  which  the 
conversion  of  5-FCyt  (O)  to  5-FU  was  not  detected  by  MRS.  h,  L6-CD  and  13B 
anti-idiotypic  antibody-pretreated  mice  (n  =  5),  in  which  the  conversion  of  5-FCyt  (•)  to 
5-FU  (■)  was  detected  by  MRS.  Mice  were  pretreated  with  L6-CD  and  13B  anti-idiotypic 
antibody  at  48  and  24  h,  respectively,  before  5-FCyt  injection.  Data  arc  means  ±  SE 
{error  bars). 
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Fig.  4.  a,  typical  proton  image  of  a  H2981 
tumor-bearing  mouse  showing  the  anatomy  of  the 
mouse,  b,  typical  CSI  images  showing  distribution 
of  5-FCyt  and  (c)  of  5-FU.  Mice  were  pretreated 
with  L6-CD  and  13B  anti-idiotypic  antibody  at  48 
and  24  h,  respectively,  before  5-FCyt  injection  (600 
mg/kg  i.p.).  *®F  CSI  was  commenced  5  min  after 
5-FCyt  injection  (total  imaging  time,  35  min;  in¬ 
plane  resolution,  6.25  X  6.25  mm).  The  octagon 
indicates  the  tumor  region. 
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Biodistribution  studies  were  performed  on  a  4.7T  spectrometer  with  a 
home-built  whole  body  loop-gap  resonator  constmcted  to  fit  the  mice. 

Mice  were  treated  as  above  except  with  600  mg/kg  5-FCyt.  As  a  prelude  to  CSI 
experiments,  multislice  proton  images  from  2-mm  slices  were  acquired  with  a 
two-dimensional  spin  echo  imaging  pulse  sequence  to  define  mouse  anat¬ 
omy  (100-mm  field  of  view,  1-s  repetition  time,  18-ms  echo  time,  4  scans,  and 
256  X  128  matrix  size).  *^F  CSI  experiments  were  started  5  min  after  drug 
injection  with  a  total  imaging  time  of  35  min.  A  two-dimensional  spin  echo 
CSI  pulse  sequence  (17)  was  used  (thick  slice,  100-mm  field  of  view,  4  scans, 
16  X  16  phase-encoding  steps,  2-s  repetition  time,  9-ms  echo  time).  Intensity 
distribution  of  each  ‘^F  peak  (5-FCyt  or  5-FU)  was  obtained  by  image 
reconstruction  using  in-house  software  (dOs)  written  by  Dr.  D.  C.  Shungu. 

Results 

Both  5-FCyt  and  5-FU  resonances,  which  have  a  chemical  shift 
difference  of  1.2  ppm,  were  detected  by  ^^F  NMR  spectroscopy. 
Using  this  technique,  the  selectivity  of  L6-CD  binding  to  H2981  cells 
was  demonstrated  in  vitro  (Fig.  la-d).  After  incubation  of  H2981  cells 
with  L6-CD  and  three  washing  cycles,  the  cells  were  able  to  convert 
5-FCyt  rapidly  to  5-FU.  Levels  of  5-FU  in  cells  and  media 
(mean  ±  SD)  after  1  h  incubation  with  3  him  5-FCyt  were  0.3  ±  0.03 
/xmol/2  X  10^  cells  and  1.62  ±  0.03  mM,  respectively.  5-FCyt  was 
only  observed  in  the  medium  at  this  time  point  (0.4  ±  0,05  mM).  In 
contrast,  conversion  of  5-FCyt  to  5-FU  did  not  occur  in  cells  treated 
with  1F5-CD.  Levels  of  5-FCyt  in  cells  and  media  after  1-h  incubation 
with  3  niM  5-FCyt  were  0.42  ±  0.02  /xmol/2  X  10^  cells  and 
2.50  ±  0.08  mM,  respectively.  When  incubated  directly  with  5-FCyt, 
1F5-CD  was  equally  capable  of  activating  5-FCyt  to  5-FU  (data  not 
shown).  This  means  that  1F5-CD  did  not  bind  to  antigens  on  H2981 
cells  and  was,  therefore,  eliminated  during  the  washing  cycles.  In 
control  experiments  in  which  cells  were  not  incubated  with  conjugate, 
the  levels  of  5-FCyt  in  cells  and  medium  after  incubation  with  3  mM 
5-FCyt  for  1  h  were  0.43  ±  0.02  /xmol/2  X  10^  cells  and  2.95  ±  0.42 
mM,  respectively.  Corresponding  control  experiments  with  3  mM  5-FU 
gave  values  of  0.36  ±  0.02  /xmol/2  X  10^  cells  and  1.96  ±  0.16  mM 
for  cells  and  media,  respectively.  Of  interest  was  the  finding  that  no 
other  signals  from  5-FU  catabolites  or  anabolites  were  observed  in 
these  experiments  within  the  experimental  time  used. 


Both  5-FCyt  and  5-FU  resonances  were  detected  nonin vasively  in 
tumors  of  mice  treated  with  L6-CD/13B/5-FCyt  with  adequate  spec¬ 
tral  resolution  (Fig.  2).  Only  the  5-FCyt  peak  was  detected  in  control 
mice.  For  each  mouse,  resonance  intensities  of  5-FCyt  or  5-FU  were 
normalized  to  the  highest  5-FCyt  intensity.  The  mean  intensities  are 
illustrated  in  Fig.  3,  a  and  b.  In  control  mice,  maximum  tumor  levels 
of  5-FCyt  occurred  early  at  10  min  after  5-FCyt  injection  and  de¬ 
creased  to  50%  of  maximum  levels  at  60  min  {n  =  3).  In  L6-CD/ 
13B-treated  mice,  maximum  5-FCyt  tumor  levels  occurred  at  25  min 
and  were  60%  of  maximum  levels  at  60  min  {n  =  5).  Maximum  5-FU 
tumor  levels  occurred  early  at  10  min,  and  the  5-FU/5-FCyt  ratio  (an 
index  of  CD  activity)  ranged  between  0.5  and  1.0  within  the  study 
time  (60  min).  On  the  basis  of  previous  studies  using  the  L6-CD/13B/ 
5-FCyt  approach  (1),  these  results  imply  that  tumor  concentrations  of 
5-FU  will  be  0.5-1. 0  molar  equivalent  of  intratumoral  5-FCyt  when 
—8.7  jLxg  of  L6-CD  (specific  activity,  80.6  units/mg)  are  present  in  1  g 
of  a  H2981  tumor.  Again,  no  other  signals  from  5-FU  catabolites  or 
anabolites  were  observed  in  these  experiments. 

Initial  “phantom”  studies  demonstrated  that  it  was  possible  to 
spatially  resolve  peaks  derived  from  5-FCyt  and  5-FU  resonances  by 
CSL"^  We  used  this  technique  to  study  the  biodistribution  of  5-FCyt 
and  5-FU  in  L6-CD/13B -treated  mice.  Fig.  4a  shows  a  typical  proton 
image  (0.4  X  0.8  mm  in-plane  resolution)  illustrating  the  anatomy  of 
a  mouse.  Corresponding  CSI  data  from  5-FCyt  and  5-FU  resonances 
in  the  same  mouse  are  displayed  in  Fig.  4,  b  and  c.  Signals  from 
5-FCyt  were  mainly  observed  in  the  abdominal  and  thoracic  regions 
including  tumor  (Fig.  4b).  Importantly,  the  highest  5-FU  intensity  was 
observed  within  the  tumor  region  (Fig.  4c). 

Discussion 

5-FU  is  frequently  used  in  combination  therapy  for  the  treatment  of 
adenocarcinomas  of  the  colon,  breast,  and  lung.  Both  antibody-CD 
conjugates  and  genes  coding  for  CD  have  been  previously  used  to 
generate  5-FU  from  5-FCyt  (1,2,  5-9).  These  novel  approaches  are 
capable  of  generating  high  and  sustained  levels  of  5-FU  selectively 
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within  mouse  tumors  and  are,  therefore,  associated  with  dramatic 
cytotoxicity,  tumor  regression,  and  even  cures  (2,  5).  A  critical  issue 
that  may  influence  the  clinical  applicability  of  these  approaches  is  the 
ability  to  detect  drug  activation  and  target  tissue  selectivity.  To 
facilitate  this,  we  have  applied  MR  methods  to  detect  selective  pro¬ 
drug  activation  in  vitro  as  well  as  to  define  the  magnitude  and  location 
of  CD  enzyme  activity  in  mouse  tissues  after  systemic  administration 
of  an  antibody-enzyme  conjugate. 

In  this  study,  we  showed  that  MRS  and  CSI  methods  are  capable  of 
detecting  conversion  of  5-FCyt  to  5-FU  in  vitro  and  in  vivo,  5-FCyt 
and  5-FU  resonances  were  adequately  resolved  (~1.2  ppm).  This 
implies  that  one  can  monitor  the  levels  of  these  compounds  directly 
without  resorting  to  methods  that  measure  total  signal  after  washout  of 
the  original  compound.  This  is  particularly  important  since  the  me¬ 
tabolism  of  5-FCyt  by  CD  was  rapid,  and  both  prodrug  and  metabolite 
had  similar  elimination  profiles  within  the  time  interval  studied  (1  h). 
We  demonstrated  the  selectivity  of  using  L6-CD  in  vitro.  L6-CD  but 
not  1F5-CD  activated  5-FCyt  to  5-FU  when  H2981  human  lung 
adenocarcinoma  cells  were  incubated  with  5-FCyt.  1F5-CD  binds  to 
CD20  antigen  expressed  on  both  normal  and  neoplastic  B  cells  but  not 
to  L6  antigens  (2,  15).  Thus,  after  three  cycles  of  washing,  no  enzyme 
conjugate  was  available  for  metabolic  activation  of  5-FCyt.  In  con¬ 
trast,  binding  of  L6-CD  to  L6  antigens  enabled  metabolic  activation 
even  after  three  washing  cycles.  It  is  possible  that  5-FU  is  further 
metabolized  in  H2981  cells,  but  despite  a  reduction  in  5-FU  levels  of 
>30%  after  1-h  incubation  of  cells  with  5-FU,  no  metabolites  were 
detected  within  20  ppm  up-  or  down-field  of  the  5-FU  resonance.  This 
may  be  due  to  the  low  levels  of  these  multiple  putative  5-FU  metab¬ 
olites. 

In  vivo  studies  were  performed  by  pretreating  tumor-bearing  mice 
with  L6-CD  followed  by  13B  anti-idiotypic  antibody  24  h  after  having 
received  L6-CD.  13B  binds  to  and  clears  circulating  L6~CD  but  not 
the  conjugate  that  is  tumor-associated.  This  has  been  shown  previ¬ 
ously  (1)  to  lead  to  high  tumor:blood  conjugate  ratios  (42:1)  and  has 
enabled  large  doses  of  5-FCyt  to  be  subsequently  administered.  In 
H2981  tumors  pretreated  with  L6-CD/13B,  metabolic  activation  of 
5-FCyt  to  5-FU  was  detected  and  monitored  for  60  min.  Metabolic 
activation  was  rapid  and  5-FU  levels  were  between  50  and  100%  of 
5-FCyt  tumor  levels.  No  activation  was  observed  in  control  mice.  CSI 
experiments  demonstrated  that  5-FCyt  localized  mainly  in  the  abdom¬ 
inal  and  lower  thoracic  regions  of  the  mice  after  i.p.  injection.  The 
highest  5-FU  levels  were  observed  in  the  tumor  region,  supporting  the 
selectivity  of  this  ADEPT  approach  for  tumors.  However,  other  re¬ 
gions  showed  signal  from  5-FU.  This  may  be  caused  by  the  diffusion 
of  5-FU  out  of  the  tumor  region  and/or  partial  volume  averaging, 
although  we  cannot  completely  rule  out  contribution  from  nonspecific 
activation.  The  absence  of  the  hepatic  catabolite  of  5-FU,  fluoro-^- 
alanine,  in  in  vivo  spectra  suggests  that  there  were  low  levels  of 
circulating  5-FU  in  the  blood.  The  absence  of  anabolites  of  5-FU, 
however,  may  be  caused  by  the  short  observation  time  of  1  h  and/or 
by  low  levels  of  such  metabolites.  The  results  obtained  in  these 
experiments  are  in  general  agreement  with  those  obtained  by  remov¬ 
ing  the  tumors  from  the  animals  and  analyzing  the  drug  concentrations 
by  high-performance  liquid  chromatography  (1). 

These  MR  methods  for  detecting  conversion  of  5-FCyt  to  5-FU  will 
aid  in  the  clinical  testing  of  CD-based  prodrug  therapies.  Clinical 
applicability  at  the  present  widely  available  magnetic  field  strength  of 
1.5  tesla  may  be  hampered  by  low  spectral  resolution.  This  problem 
will  be  solved  by  employing  high-field  magnets  such  as  3,  4,  and  9 


tesla  presently  available  in  some  institutions.  MR  methods  are  inher¬ 
ently  insensitive.  However,  at  the  5-FCyt  dose  (400  mg/kg)  and  the  ^ 
magnetic  field  strength  (4.7  tesla)  used,  the  compounds  are  easily  ^ 
detected.  Lower  doses  of  5-FCyt  (100  mg/kg)  can  be  detected  at  this 
field  strength,  and  sensitivity/spatial  resolution  can  be  further  en-  * 
hanced  by  increasing  the  number  of  scans  or  imaging  time.  In  con¬ 
clusion,  we  have  demonstrated  feasibility  for  using  MRS  and  CSI 
methods  to  noninvasively  detect  the  activation  of  5-FCyt  to  5-FU  by 
ADEPT.  These  MR  methods  can  be  applied  to  the  preclinical  and 
early  clinical  testing  of  CD/5-FCyt  therapy.  In  addition  MRS  and  CSI 
may  have  applications  in  a  variety  of  other  enzyme-prodrug  combi¬ 
nations  of  clinical  interest. 
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Summary;  Angiogenesis,  the  process  by  which  new  blood  vessels  are  generated,  occurs 
during  wound  healing,  in  the  female  reproductive  system  during  ovulation  and  gestation, 
and  during  embryonic  development  The  process  is  carefully  controlled  with  positive  and 
negative  regulators,  because  several  vital  physiological  functions  require  angiogenesis. 
The  consequences  of  abnormal  angiogenesis  are  either  excessive  or  insufficient  blood 
vessel  growth.  Ulcers,  strokes,  and  heart  attacks  can  result  from  the  absence  of  angio¬ 
genesis  normally  required  for  natural  healing,  whereas  excessive  blood  vessel  prolifera¬ 
tion  may  favor  tumor  growth  and  dissemination,  blindness,  and  arthritis.  In  this  review, 
the  process  of  angiogenesis  and  the  characteristics  of  the  resulting  tumor  vasculature  are 
outlined.  Contrast-enhanced  magnetic  resonance  imaging  techniques  that  currently  are 
available  for  basic  research  and  clinical  applications  to  study  various  aspects  of  tumor 
angiogenesis  and  neovascularization  are  discussed.  Key  Words:  Angiogenesis — 
Vascularization — Cancer — Contrast-enhanced  magnetic  resonance  imaging. 


BRIEF  HISTORICAL  PERSPECTIVE  OF 
TUMOR  VASCULATURE 

Of  all  the  complexities  and  protean  facets  that  cancer 
displays,  the  ability  of  a  colony  of  cancer  cells  to  recruit, 
form,  establish,  and  maintain  a  vascular  network  is  one  of 
the  most  fascinating.  Interest  in  tumor  vasculature  has  ex¬ 
isted  since  the  1920s  and  earlier.  For  instance,  in  1907, 
Goldmann  (1)  noted  that  the  development  of  a  tumor  in 
the  liver,  stomach,  or  any  other  organ  resulted  in  “chaotic 
irregularity”  of  the  blood  vessels  and  that  a  growing  tumor 
exhibited  an  extensive  formation  of  blood  vessels  that  was 
most  apparent  in  the  zone  of  proliferation  that  for  infil¬ 
trating  tumors,  is  at  the  periphery.  He  also  observed  that 
necrosis  occurred  in  highly  vascularized  growths,  and  he 
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made  the  unique  suggestion  that  areas  of  necrosis  were 
"the  battlefield  on  which  assailant  and  defendant  both  per¬ 
ished."  Similarly,  it  was  recognized  as  early  at  1945  by 
Algire  ct  al.  (2)  that  malignant  cells  provoked  a  continued 
vascular  proliferation. 

Research  through  the  subsequent  decades  has  only 
served  to  further  reveal  the  undeniable  importance  of  tu¬ 
mor  vasculature  not  only  for  tumor  growth  and  hemato- 
geneous  dissemination,  but  also  for  cancer  therapy.  Ob¬ 
servations  that  had  a  significant  impact  on  the  field  of 
radiation  therapy  were  made  by  Thomlinson  and  Gray  (3) 
on  human  specimens  of  bronchogenic  carcinoma  in  1955. 
In  histologic  sections  of  bronchogenic  carcinoma,  necrosis 
was  found  to  occur  at  a  distance  of  1 60  pim  or  greater  from 
the  nearest  vessel.  This  distance  corresponded  closely  to 
the  diffusion  distance  calculated  by  them,  where  the  con¬ 
centration  of  oxygen  would  approach  zero.  Their  obser¬ 
vations  provided  some  of  the  earliest  indirect  evidence  of 
the  existence  of  hypoxia  and  hence  radioresistant  cells  in 
the  tumors  (4). 

Characteristic  “hallmarks"  of  tumor  vasculature  have 
emerged  from  numerous  studies,  which  now  span  almost 
a  century  (for  reviews  see  references  5-7).  Notable  among 
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these  characteristics  of  tumor  vasculature  are  (1)  spatial 
heterogeneity  and  chaotic  structures;  (2)  arteriovenous 
shunts;  (3)  acutely  collapsed  vessels  and  transiently  col¬ 
lapsing  vessels;  (4)  poorly  differentiated,  fragile,  and 
leaky  vessels  lacking  in  smooth  muscle  cell  lining;  and  (5) 
vasculature  that  frequently  is  unable  to  match  the  rapid 
growth  of  cancer  cells,  resulting  in  areas  of  hypoxia  and 

©necrosis.  In  an  influential  review,  Rubin  and  Ca/saretl  (8) 
pointed  out  that  although  the  ability  of  a  tumor  cell  to 
survive  depends  on  its  distance  from  the  capillary  from 
which  essential  substrates  and  oxygen  diffuse,  necrosis  is 
not  always  due  to  absence  of  vascularization  or  due  to  the 
tumor  outgrowing  its  blood  supply,  a  sentiment  first  ex¬ 
pressed  by  Goldmann  (1)  in  1907.  Acute  vascular  collapse 
due  to  high  tumor  interstitial  pressure  (6,9,10),  the  im¬ 
mune-related  destruction  of  cells,  and  the  release  of  toxic 
agents  (11)  have  been  cited  as  some  of  the  other  causes. 

Vascular  properties  of  tumors  have  always  been  con¬ 
sidered  important  because  of  their  role  in  the  delivery  of 
therapeutic  agents,  and  because  radioresistant  hypoxic 
cells  may  be  a  potential  source  of  treatment  failure  fol¬ 
lowing  radiation  therapy.  The  recent  overwhelming  inter¬ 
est  in  tumor  vascularization  and  angiogenesis,  however, 
has  arisen  due  to  two  findings.  One  of  these  is  that  high 
vascular  density  in  histologic  specimens  may  be  predictive 
of  the  disposition  of  the  tumor  to  metastasize  (12-16). 
Because  angiogenesis  is  critical  for  tumor  growth,  antian- 
giogenic  treatment  as  well  as  antineovasculature  therapies 
have  been  proposed  and  investigated  for  cancer  treatment 
(17,18).  The  second  finding,  related  to  antiangiogenic 
therapy,  is  that  repeated  treatment  cycles  with  the  angio¬ 
genesis  inhibitor  endostatin  induced  tumor  dormancy 
without  inducing  tumor  drug  resistance  in  at  least  three 
experimental  tumor  models  (19). 


tance;  however,  once  a  tumor  has  reached  this  threshold 
stage  of  growdi,  vascularization  is  essential  for  additional 
growth  and  progression  (17,22,23).  The  normal  prolifera¬ 
tion  rate  of  endothelial  cells  is  very  low  compared  to  many 
other  cell  types  in  the  body,  and  the  turnover  time  of  these 
cells  can  exceed  1,(XX)  days.  In  contrast,  tumor  capillary 
endothelial  cells  have  a  potential  doubling  time  of  2.14  to 
13  days  (24,25). 

Tumor  cells,  as  mentioned  earlier,  have  the  ability  to 
induce  neovascularization  through  tlie  release  of  angio¬ 
genesis  factors  (17,22,23).  However,  not  only  cancer  cells 
but  also  inflammatory  cells  that  infiltrate  the  tumor  and 
the  extracellular  matrix  can  be  a  source  of  angiogenesis 
factors  (26).  Also  as  mentioned  earlier,  there  are  positive 
and  negative  regulators  of  angiogenesis.  Several  factors 
have  now  been  identified  as  angiogenic  promoters  and 
angiogenic  inhibitors  (26),  and  several  of  these  promoters 
and  inhibitors  are  listed  in  Table  1.  The  phenomenon  of 
tumor  angiogenesis  begins  with  the  degradation  of  the 
basement  membrane  of  the  parent  vessel  by  proteases  se¬ 
creted  by  activated  endothelial  cells  that  migrate  and  pro¬ 
liferate,  leading  to  the  formation  of  solid  endothelial  cell 
sprouts  into  the  stromal  space  in  response  to  the  angio¬ 
genic  factors.  Vascular  loops  are  formed  and  capillary 
tubes  develop,  with  formation  of  tight  junctions  and  de¬ 
position  of  new  basement  membrane  (17,22,23,26). 

An  angiogenic  factor  that  has  attracted  significant  in¬ 
terest  is  vascular  endothelial  growth  factor  (VEGF).  Of  all 
the  various  angiogenesis  promoters  and  inhibitors,  VEGF 
appears  to  be  the  most  responsive  to  the  abnormal  physi¬ 
ological  environments  of  hypoxia,  extracellular  acidosis, 
and  substrate  deprivation  that  occur  in  solid  tumors  (27). 
Current  information  on  VEGF  has  been  reviewed  exten¬ 
sively  by  Neufeld  et  al.  (28).  VEGF  induces  angiogenesis 
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(5^ 


Although  it  was  recognized  from  very  early  on  that 
tumor  cells  secrete  substances  to  induce  vascularization 
(20),  Folkman  and  his  colleagues  (21)  were  the  first  to 
isolate  a  tumor  angiogenic  factor  in  1971.  Studies  using 
transgenic  mice  now  suggest  a  model  of  tumorigenesis  in 
which  mutated  oncogenes  and  tumor  repressor  genes  in¬ 
duce  cell  proliferation  and  hyperplasia.  There  appears  to 
be  a  secondary  set  of  events  after  cell  proliferation  leading 
from  hyperplasia  to  neoplasia  (17,22,23).  The  induction  of 
angiogenesis  is  one  of  the  critical  events  or  “switches”  in 
this  progression,  given  that  the  prevascular  phase  can  per¬ 
sist  for  years  with  limited  growth  of  the  lesion  (23). 

Tumor  growth  to  about  1  mm  can  be  accomplished  in 
the  absence  of  neovascularization,  because  all  essential 
nutrients  and  waste  products  can  diffuse  across  this  dis- 


TABLE  1. 

Angiogenesis  promoters  Angiogenesis  inhibitors 


Fibroblast  growth  factors 
Placental  growth  factor 
Vascular  endothelial  growth  factor 
Transfonning  growth  factor  a 
Traasforming  growth  factor  p 
Angiogenin 
Interleukin  8 

Hepatocyte  growth  factor 
Platelet-derived  endothelial 
ceil  growth  factor 
Angiopoietin  1 

Granulocyte  colony-stimulating 
factor 

Wound  fluid 
Epidermal  growth  factor 
Tumor  necrosis  factor  a 
Prostaglandins 


Thrombospondin  t 

Angiostatin 

Endostatin 

Prolactin  16-kD  fragment 
Metalloproteinase  inhibitors 
Platelet  factor  4 
Genistein 

Transfonning  growth  factor  ^  (?) 
Interferon  a 

Placental  proliferin-related  protein 

Fumagiltol 

Retinoids 

Protamine 
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(29,30)  and  is  a  potent  vascular  permeability  factor  (31). 
There  are  five  isoforms  of  VEGF  (with  121,  145,  165, 
189,  and  206  amino  acids)  generated  from  alternative 
splicing  of  the  VEGF  gene  that  are  distinguishable  by  their 
heparin  and  heparin-surface  binding  properties,  water 
solubility,  and  target  receptors  (28).  The  various  forms  of 
VEGF  bind  to  two  tyrosine-kinase  receptors  {flt-I)  and 
(KDR/flk-I),  which  are  expressed  almost  exclusively  on 
endothelial  cells.  In  addition,  endothelial  cells  express 
neuropilin-1  and  neuropilin-2  corcceptors,  which  bind  se¬ 
lectively  to  the  165  amino  acid  form  of  VEGF.  Of  the 
various  isoforms  of  VEGF,  the  121  and  165  forms  pre¬ 
dominate,  although  the  189  form  also  is  found  in  most 
VEGF-producing  cell  types  (32).  Expression  of  VEGF145<:— 
is  restricted  to  cells  derived  from  reproductive  organs  (28). 
Hypoxia-mediated  transcription  of  VEGF  is  thought  to 
occur  through  the  binding  of  the  hypoxia- inducible  factor 
1  (HIF-la)  to  an  HEF-l  binding  site  located  in  the  VEGF 
promoter  (33).  Hypoxia  also  promotes  the  stability  of 
VEGF  mRNA  through  binding  proteins  (34).  There  ap¬ 
pears  to  be  a  dependence  of  VEGF  expression  on  COX-2 
and  iNOS  that  is  mediated  through  prostaglandin  E2  and 
nitric  oxide  (35).  Under  hypoxic  conditions,  which  are 
typical  for  most  solid  tumors,  VEGF  expression  is  regu¬ 
lated  by  pH  and  glucose  concentration  (36,37).  VEGF 
expression  also  can  be  modulated  by  oncogenes  (28)  and 
by  the  tumor  suppressor  gene  p53  (39,40)  although  results 
to  the  contrary  have  been  observed  (41). 

Non  tWACc  NBTM^ 
SOME  ASPECTS  OF  NOIS^WAGNETIC 

RESONANCE  BIOASSAYS  FOR  MEASURING 
ANGIOGENESIS  AND  VASCULARIZATION 

The  most  commonly  used  bioassays  for  angiogenesU 
are  (1)  the  developing  chick  chorioallantoic  membran^^ 
(2)  the  comeal  pocket  assay,  and  (3)  the  transparent  Mil- 
lipore  chamber  model  grown  in  the  rabbit  ear  or  hamster 
cheek  pouch  (20,42,43).  The  CAM  bioassay  is  carried  out 
in  fertilized  chick  eggs  in  which  a  window  is  prepared  by 
removal  of  a  section  of  egg  shell  or  in  which  the  shell  is 
totally  removed  and  the  egg  cultured  in  a  Petri  dish.  Test 
substances  are  placed  onto  the  chick  chorionic  membrane 
of  8-  to  10-day-old  fertilized  chick  egg  on  filters  or  cov- 
erslips.  Neovascularization  is  monitored  visually,  1-2 
days  after  implantation,  and  the  membranes  are  fixed  and 
evaluated  histologically.  In  the  rabbit  comeal  assay,  test 
substances  are  incorporated  into  sustained  release  polymer 
pellets  that  are  implanted  into  a  surgical  pocket  prepared 
in  the  cornea.  With  a  positive  response  there  is  directional 
growth  of  new  capillaries  from  the  limbal  blood  vessel 
toward  the  implant.  The  bioassay  can  be  quantified  by 
measuring  the  rate  of  vessel  growth,  as  well  as  the  length 


and  the  density  of  the  new  blood  vessels.  Angiogenesis  in 
the  transparent  Millipore  chamber  is  evaluated  similarly. 
These  techniques  have  been  used  for  pioneering  studies  in 
angiogenesis,  but  some  of  the  disadvantages  are  the  limi¬ 
tations  on  the  number  of  samples  that  can  be  assayed  and 
the  difficulty  in  quantifying  the  results. 

The  initial  angiogenic  stimulus  results  in  the  formation 
of  tumor  neovasculature  necessaiy  to  establish  the  solid 
tumor  and  deliver  substrates  and  nutrients.  Angiogenesis, 
continues  to  occur  through  the  life  time  of  the  tumor, 
because  new  vasculature  has  to  be  continually  generated 
with  tumor  growth.  An  example  of  this  is  shown  in  Fig.  I. 
Thus,  at  any  given  time  in  the  history  of  a  growing  solid 
tumor,  there  will  be  established  vasculature  and  vascula¬ 
ture  in  the  process  of  being  formed.  Three  quantities,  tu¬ 
mor  blood  flow,  tumor  vascular  volume,  and  vascular  per¬ 
meability,  have  mainly  been  used  to  characterize  tumor 
vasculature.  Several  methods  arc  available  to  measure  tu¬ 
mor  blood  flow  quantiatively;  however,  most  of  these 
methods,  with  the  exception  of  positron  emission  tomog¬ 
raphy  are  either  invasive  or  lack  spatial  information.  The 
reader  is  referred  to  three  excellent  reviews  on  measure¬ 
ment  and  characterization  of  tumor  blood  flow  (7,27,44). 

Some  of  the  earliest  measurements  of  percent  vascular 
volume  were  obtained  by  morphometric  estimation  of  ves¬ 
sels  in  histologic  sections  (2,45).  The  functionality  of  ves¬ 
sels  cannot  be  evaluated  in  histologic  ^sections,  and  it  may 
be  that  a  fraction  of  the  vessels  nonfunctional.  This 

may  result  in  a  discrepancy  between  results  obtained  using 
morphometric  analysis  and  results  obtained  using  mag¬ 
netic  resonance  imaging  (MRI)  (discussed  later),  which 
measures  vasculature  into  which  the  contrast  agent  can  be 
delivered.  This  issue  becomes  increasingly  important  as 
immunoperoxidase  staining  of  molecules  such  as  CD3 1 
and  CD34  (46),  which  are  present  on  the  surface  of  en¬ 
dothelial  cells,  are  used  to  identify  vessels.  It  also  should 
be  mentioned  here  that  in  studies  that  have  shown  a  de¬ 
pendence  of  metastasis  on  microvessel  density,  the  analy¬ 
sis  usually  is  performed  by  identifying  the  area  of  highest 
neovascularization  within  the  section  and  obtaining  mi¬ 
crovessel  counts  in  those  areas.  In  the  study  by  Weidner  et 
al.  (12),  for  patients  with  invasive  breast  carcinomas  with 
any  metastases,  the  mean  microvessel  counts  for  the  areas 
of  highest  neovascularization  were  101  (at  200x  field) 
compared  to  45  in  carcinoma  from  patients  without  me¬ 
tastases.  Thus,  the  number  of  microvessels  in  the  regions 
of  highest  neovascularization  determines  risk  of  metasta¬ 
sis.  Several  macrosopic  and  microscopic  methods  of  mea¬ 
suring  vascular  volume  have  been  developed  (471|.that 
mthods  usually  require  administration  of  radioactive  trac¬ 
ers  and  excision  of  the  tumor  (48). 

Different  mathematical  formulae  have  been  used  to 
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(a)  (b) 

FIG.  1.  Htgh>power  micrograi^  (x400)  of  sprouting  endothelial  cells  in  two  adjacent  S*p.m  thick  sections  stained  (A)  with  hematoxylin  and  eosin  and 
(B)  immunostained  for  VEGF  expression.  Sections  were  obtained  from  a  lO-mm-diameter  MDA-MB-231  human  breast  cancer  model  growing  in  the 
flank  of  a  SCID  mouM. 


compute  vessel  permeability  (5).  Usually  an  effective 
“permeability-surface  area”  product  (PS),  also  referred  to 
as  “capillary-diffusion  capacity,”  is  obtained  by  determin¬ 
ing  the  rate  of  extravasation  of  the  test  substance  into  the 
tumor.  PS  values  quoted  in  literature  depend  on  the  mo¬ 
lecular  size  of  the  test  substance  used;  the  reader  is  re¬ 
ferred  to  excellent  reviews  on  tumor  vascular  permeability 
(5,48). 

MRI  OF  ANGIOGENESIS  AND 
VASCULARIZATION 

Angiogenesis  is  a  dynamic  process.  A  method  to  quan¬ 
titate  neovascularization  in  situ  without  perturbing  or  ex¬ 
cising  the  tumor  would  be  useful  to  examine  the  temporal 
relationship  between  neovascularization  and  tumor  growth 
and  determine  the  effects  of  antiangiogenic  therapy.  A 
noninvasive  method  of  obtaining  an  index  of  vasculariza¬ 
tion  or  microvesscl  density  also  would  be  advantageous 
clinically  to  (1)  determine  the  effects  of  preoperative  hor¬ 
monal  deprivation  on  vascularization,  (2)  determine  the 
aggressiveness  of  the  lesion,  and  (3)  play  a  role  in  the 
selection  and  monitoring  of  patients  treated  with  antian¬ 
giogenic  drugs.  MRI  is  suitable  for  these  purposes  because 
MR  images  can  be  obtained  noninvasivety  with  an  inplane 
resolution  on  the  order  of  100  p.m.  Vascularization  can  be 
evaluated  using  intrinsic  endogenous  contrast  (49)  or  with 
exogenous  paramagnetic  contrast  agents.  Gadolinium  che¬ 
lates  are  the  contrast  agents  frequently  used  for  MRI. 
These  agents  are  tightly  bound  complexes  of  the  rare  earth 


element  gadolinium  (Gd)  and  various  chelating  agents. 
The  seven  unpaired  electrons  of  Gd  produce  a  large  mag¬ 
netic  moment  that  results  in  paramagnetic  properties  when 
the  ion  is  placed  in  a  magnetic  field  and  creates  contrast  in 
an  MR  image.  Paramagnetic  agents  shorten  both  T1  (spin- 
lattice  relaxation  time )  and  T2  (spin-spin  relaxation  time), 
but  because  tissue  T2  values  are  intrinsically  short,  the  T1 
effect  of  the  contrast  agent  predominates;  tissues  that  take 
up  a  paramagnetic  agent  are  brightened  (positive  enhance¬ 
ment).  Susceptibility  effects  of  Gd-based  contrast  agents 
resulting  in  enhancement  of  T2  and  T2*  relaxation  have 
been  used  to  measure  vascular  volume  and  flow  in  clinical 
studies  in  brain  (50).  These  studies  have  been  reviewed 
elsewhere  (51)  and  are  not  discussed  here  because  their 
application  to  tumors  is  limited  due  to  leaky  vasculature 
and  relatively  slow  perfusion  rate  in  tumors. 

The  tissue  concentration  of  Gd-based  agents  can  be 
calculated  from  changes  in  T1  obtained  from  quantitative 
T1  maps  or  Tl-weighted  images.  Several  Gd  complexes 
are  under  development  or  in  use.  Low-molecular-weight 
Gd-DTPA  compounds  (0.57  kD)  are  used  in  clinics  for 
contrast  enhancement  of  various  lesions  including  malig¬ 
nant  tumors.  Macromolecular  contrast  agents,  such  as  al- 
bumin-(Gd-DTPA),  remain  in  the  intravascular  space  with 
a  half-life  of  several  hours  because  of  their  large  molecu¬ 
lar  weight  (approximately  65  kD).  Thus,  analysis  of  re- 
laxivity  changes  induced  by  an  intravascular  agent  can  be 
used  to  determine  blood  or  vascular  volume  and  vascular 
permeability  (52,53).  A  novel  and  promising  area  is  the 
design  of  receptor- targeted  contrast  agents  where  the  re- 
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suiting  contrast  is  dependent  on  the  density  of  receptor 
expression  (54).  In  the  following  sections  we  present  brief 
descriptions  of  current  MR  methods  used  to  characterize 
tumor  vascularization  based  on  the  physical  properties  of 
the  contrast  agent  used.  Of  the  three  parameters  mainly 
used  to  characterize  tumor  vasculature,  currently  only  vas¬ 
cular  volume  and  vascular  permeability  can  be  quantified 
with  contrast-enhanced  proton  imaging.  Quantitative 
blood  flow  measurements  have  been  obtained  using  deu¬ 
terium  MRI  with  D2O  as  the  blood  flow  tracer  (55,56).  A 
detailed  discussion  of  this  technique  is  beyond  the  scope 
of  the  present  review  which  is  restricted  to  proton  imag¬ 
ing. 

Assessment  of  angiogenesis  with  gradient-echo  MRI 
Detection  of  tumor  vasculature  using  the  intrinsic  T2* 
contrast  produced  by  deoxygenated  hemoglobin  in  tumor 
capillaries  (49)  is  based  on  the  BOLD  technique  originally 
proposed  by  Ogawa  ct  al.  (57)  for  brain  imaging.  The 
functional  dependence  of  T2*  is  given  by  Equation  1 : 

1/T2*  a  (1  -  Y)b,  (7) 

where  Y  is  the  fraction  of  oxygenated  blood  and  b  is  the 
partial  vascular  or  blood  volume  fraction.  In  hypoxic  tu¬ 
mors  where  0  <  y  <  0.2,  the  contrast  produced  by  the 
method  is  primarily  dependent  on  6,  which  permits  rapid 
in  vivo  determination  of  vessel  densities  in  tumors.  The 
method  works  best  in  poorly  oxygenated  tumors  such  as 
subcutaneous  models  and  human  xenografts  with  random 
orientation  of  sprouting  capillaries.  The  technique  has 
been  used  to  detect  changes  in  tumor  vascularization  fol¬ 
lowing  induction  of  angiogenesis  by  exogenous  angio¬ 
genic  compounds  (49)  and  to  obtain  '’functional  vascular'* 
maps  of  genetically  modified  HIF-1  (+/+  and  -/-)  animal 
tumor  models  (58).  The  method  provides  fast  and  truly 
nonin  vast  ve  of  measurement  of  tumor  vascularization  be¬ 
cause  it  does  not  require  injection  of  contrast  agent  How¬ 
ever,  it  cannot  provide  quantitative  measurements  of  tu¬ 
mor  vascular  volume,  and  it  cannot  measure  vascular 
permeability  or  blood  flow. 

MRI  detection  of  tumor  angiogenesis  by 
receptor-targeted  contrast  agents 
The  feasibility  of  imaging  newly  formed  vasculature  in 
squamous  cell  carcinoma  in  rabbits  was  demonstrated  re¬ 
cently  by  Sipkins  et  al.  (54)  using  a  paramagnetic  contrast 
agent  targeted  to  the  endothelial  integrin.  The  con¬ 
trast  agent  was  composed  of  Gd-labeled  polymerized  li¬ 
posomes  conjugated  with  biotinylated  antibodies  targeted 
against  ayP3  receptors  on  endothelial  cells.  A  selective 
enhancement  of  the  tumors  was  demonstrated  without  any 
significant  binding  of  the  contrast  agent  to  normal  muscle 


vasculature  adjacent  to  the  tumor.  The  large  particle  di¬ 
ameter  of  300-350  nm  provided  a  prolonged  life  time  of 
the  contrast  agent  in  plasma  lasting  approximately  8  h, 
allowing  accumulation  of  the  contrast  agent  at  the  otyp3 
receptors.  Image  intensity  proportional  to  the  density  of 
the  receptors  was  obtained  from  standard  T1 -weighted 
images.  The  method  provides  information  about  the  den¬ 
sity  of  neovasculature,  but  dynamic  vascular  characteris¬ 
tics  such  as  permeability  or  flow  cannot  be  measured. 

Low-molecular-weight  contrast  agents 

Because  low-molecular-weight  contrast  agents  are  the 
only  class  of  paramagnetic  agents  approved  for  routine 
clinical  use,  there  are  multiple  reports  describing  applica¬ 
tions  of  these  agents  to  image  a  variety  of  tumors  includ¬ 
ing  breast  tumors  (59-61),  brain  tumors  (51),  and  uterine 
cervical  carcinomas  (62).  Similarly,  several  theroretical 
models  have  been  developed  to  derive  and  quantify  tumor 
vascular  parameters  for  these  agents. 

A  frequently  used  theoretical  model  of  analysis  for  low- 
molecular-weight  contrast  agents  is  that  of  Tofts  (63), 
which  mainly  relies  on  the  existence  of  vascular  and  in¬ 
terstitial  compartments  with  uniform  tracer  concentrations 
within  the  compartments.  The  model  also  assumes  that 
transport  of  the  agent  across  the  capillaaiy  wall  is  not  flow 
limited.  The  model  ultimately  derives  influx  and  outflux 
transfer  constants  and  extravascular  extracellular  space. 
The  arterial  input  function  required  by  the  model  can  be 
measured  separately  or  defined  in  real  time  using  voxels 
localized  within  large  blood  vessels.  The  analytical  solu¬ 
tion  of  the  model  can  be  derived  by  approximating  the 
arterial  input  function  to  a  multiexponential  decay  as  first 
described  by  Ohno  et  al.  (64).  The  concentration  of  Gd- 
DTPA  is  measured  from  changes  in  the  T1  relaxation  rate 
assuming  that  water  is  in  fast  exchange. 

If  the  flow  rate  is  low  and  unable  to  replenish  the  leak 
of  contrast  agent  from  the  vasculature  and  the  extraction 
fraction  of  the  agent  is  high,  i.e.,  the  fraction  of  agent 
leaving  the  vasculature  during  a  single  pass  of  blood 
through  the  tissue  is  close  to  1,  which  is  likely  for  low- 
molecular-weight  compounds  and  highly  permeable  tu¬ 
mors,  then  only  the  product  of  {blood  flow  •  extraction 
fraction)  can  be  defined  using  the  model  (65).  However,  if 
flow  can  be  measured  independently  then  the  extraction 
fraction  (E)  and  PS  can  be  defined  and  quantified.  Kovar 
et  al.  (66)  recently  used  this  approach  where  the  freely 
diffusable  tracer  D2O  was  used  to  measure  flow  indepen¬ 
dently  of  Gd-DTPA-dimeglumine,  which  was  used  to  de¬ 
termine  the  extraction  fraction.  The  method  is  somewhat 
restricted  by  the  low  sensitivity  of  deuterium  MRI  result¬ 
ing  in  low  spatial  resolution  and  because  D2O  cannot  be 
used  clinically.  Combining  MRI  with  an  imaging  modality 
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such  as  positron  emission  tomography  can  provide  an  al¬ 
ternative  approach  for  quantitative  measurements  of  PS  in 
human  tumors. 

Another  approach,  developed  by  Degani  et  al.  (61),  uses 
a  three-time-point  method  to  obtain  wash-in  and  wash-out 
patterns  of  the  contrast  agent.  The  patterns  arc  displayed 
as  an  artificial  color  map  with  the  intensity  and  hue  of  the 
colors  proportional  to  the  patterns.  The  method  is  based  on 
the  two-compartment  flow-independent  model  proposed 
by  the  authors  previously  (60).  In  the  current  version  the 
method  is  model-free  and  can  be  used  to  differentiate  be¬ 
tween  malignant  and  benign  breast  lesions,  even  for  situ¬ 
ations  where  the  original  kinetic  model  fails  to  character¬ 
ize  the  contrast  uptake  curves.  The  method  requires 
obtaining  three  T1 -weighted  maps  (one  precontrast  and 
two  postcontrast)  that  can  be  acquired  with  high  spatial 
resolution.  However,  using  the  hue  and  intensity  of  the 
color  maps  to  independently  quantify  vascular  permeabil¬ 
ity,  vascular  volume  or  flow  rate  currently  is  not  feasible, 

Macromolecular  contrast  agents 

Quantitative  determination  of  parameters  of  tumor  vas¬ 
culature  with  low-molecular- weight  contrast  agent  is  com¬ 
plicated  by  fast  extravasation  of  the  contrast  agent  from 
leaky  capillaries  in  tumors.  The  availability  of  high- 
molecular-weight  contrast  agents  such  as  albumin-Gd- 
DTPA  complexes  or  synthetic  compounds  such  as  poly- 
lysine-Gd-DTPA  and  gadomer-17  provides  a  unique 
opportunity  for  quantitative  determination  of  tumor  vas¬ 
cular  volume  and  vascular  PS  for  molecules  of  compa¬ 
rable  sizes  (52,67)  The  relatively  slow  leakage  of  these 
agents  from  the  vasculature  results  in  a  long  half-life  time 
and  complete  equilibration  of  plasma  concentrations 
within  the  tumor  independently  of  blood  flow,  which  in 
this  case  is  large  compared  with  PS.  Assuming  fast  ex¬ 
change  of  water  between  all  the  compartments  in  the  tu¬ 
mor  (plasma,  interstitium,  and  cells),  the  concentration  of 
the  contrast  agent  within  any  given  voxel  is  proportional 
to  changes  in  relaxation  rate  (I/Tl)  before  and  after  ad¬ 
ministration  of  the  contrast  Relaxation  rates  can  be  mea¬ 
sured  either  directly  using  fast  single-shot  TI  methods 
(53)  or  from  Tl -weighted  steady-state  experiments  (68) 
that  provide  better  temporal  resolution  but  are  susceptible 
to  experimental  artifacts  caused  by  variations  in  T2  and 
T2*  relaxation  times.  Pixel-by-pixel  maps  can  be  recon¬ 
structed  from  the  acquired  data  and  processed  with  an 
appropriate  model  to  obtain  spatial  maps  of  tumor  vascu¬ 
lar  volume  and  vascular  permeability  surface  area  product 

A  simple  linear  compartment  model  describing  uptake 
of  the  contrast  agent  from  plasma  postulates  a  negligible 
reflux  of  contrast  agent  from  the  interstitium  back  to  the 
blood  compartment  Blood  concentrations  of  the  contrast 


agent  can  be  approximated  to  be  constant  for  the  duration 
of  the  MR  experiment  and,  under  these  conditions,  con¬ 
trast  uptake  is  a  linear  function  of  time  (48,69,70).  On  a 
plot  of  contrast  agent  concentration  versus  time,  the  slope 
of  the  line  provides  the  parameter  PS,  and  the  intersect  of 
the  line  with  the  vertical  axis  at  time  zero  provides  vas¬ 
cular  volume.  For  quantitative  determination  of  these  pa¬ 
rameters,  the  change  in  relaxation  rate  of  the  blood  also 
should  be  quantified.  Changes  in  blood  TI  can  be  obtained 
separately  for  blood  samples  taken  before  injection  of  the 
contrast  agent  and  at  the  end  of  the  experiment  The 
method  is  anaologous  to  the  measurement  of  vascular  vol¬ 
ume  using  radiolabeled  proteins,  usually  *^*I-albumin  (48, 
69,70).  However,  as  pointed  out  by  Tozer  and  Morris  (71), 
the  invasive  method  requires  excision  of  the  tumor,  and 
there  may  be  a  possibility  of  overestimating  vascular  vol¬ 
ume  because  the  protein  may  already  have  extravasaied 
within  2  min  usually  required  for  equilibration  of  the  pro¬ 
tein  within  plasma. 

A  more  complex  two-compartment  model,  which  takes 
into  account  reflux  of  the  contrast  agent  and  changing 
concentrations  of  the  contrast  agent  in  plasma,  was  devel¬ 
oped  by  Su  et  al.  (72)  and  Brasch  et  al,  (68).  The  model 
used  by  Su  et  al.  (72)  derives  three  parameters — vascular 
volume,  vascular  permeability,  and  a  reflux  constant. 
Brasch  et  al.  (68)  use  a  simple  exponential  approximation 
of  clearance  of  the  contrast  agent  from  blood,  and  their 
model  does  not  require  determining  the  arterial  input  func¬ 
tion.  Five  parameters  (vascular  volume,  fractional  leakage, 
reflux  rates,  the  initial  concentration  of  contrast  agent  in 
plasma,  and  the  fractional  turnover  rate  or  clearance  of  the 
contrast  agent  from  blood)  can  be  derived  from  the  non¬ 
linear  fitting  of  tissue  concentrations  of  the  contrast  agent 
obtained  experimentally.  PS  can  be  obtained  from  the 
product  of  {vascular  volume  ^  fractional  leakage  rate). 
Blood  concentrations  of  the  contrast  agent  can  be  defined 
from  large  blood  vessels  within  the  field  of  view  of  the 
image  (68),  serial  blood  samples  collected  from  an  arterial 
line  at  different  time  points  (72),  or  from  signal  enhance¬ 
ment  of  reference  tissue  such  as  liver  (73). 

For  large  macromolecular  agents  such  as  albumin  Gd- 
DTPA  (molecular  weight  approximately  65  kD),  blood 
concentrations  equilibrate  within  2-3  min  and  do  not 
change  for  at  least  60  min  after  intravenous  bolus  injec¬ 
tion.  Tissue  concentrations  of  the  agent  for  a  time  period 
starting  5  min  after  the  bolus  required  for  plasma  equili¬ 
bration  ai^U|^Jo^40  min  after  the  bolus  increase  linearly 
with  time.^e  simple  linear  model  is  preferable  for  analy-  ^ 
sis  of  intrinsically  noisy  relaxation  data,  because  this 
model  is  much  more  more  stable  in  comparison  with  non¬ 
linear  fitting  algorithms  required  for  the  two-compartment 
models  discussed  earlier.  Accuracy  of  measurement  of  the 
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pre-^albGdDTPA 


Time  point  2 
post-albGdDTPA 


Time  point  3 
post-aIbGdDTPA 


FIG.  2.  Saturation  rccovcfy  spin-echo  images  showing  the  uptake  and  distribution  of  albumin  Cd-DTPA  in  a  subcutaneous  flank  MDA-MB-23 1  tumor 
(volume  154  mm^)  before  and  at  18,  34,  and  50  min  after  injection  (IV)  of  albumin-Gd-DTPA.  A  small  glass  capillary  containing  water  doped  with 
Gd-DTPA  can  be  observed  in  the  images.  Images  are  displayed  for  thiw  saturation  recovery  lime  intervals  (1(X)  ms,  500  ms,  1  sec)  for  a  7(X)-fim  sagittal 
slice  through  the  center  of  the  tumor  (Held  of  view  =  16  mm;  128  x  128  matrix;  number  of  averages  =  2;  sweep  width  =  20,000;  in*plane  resolution 
—  125  pm;  repetition  time  =  100  ms;  echo  time  =»  16  ms).  Images  are  displayed  scaled  so  that  the  capillary  intensity  is  constant  in  alt  the  images  to 
allow  easier  comparison  of  the  image  intensities. 


tissue  vascular  volume  does  depend,  however,  on  the  wa¬ 
ter  exchange  rate  between  the  vascular  and  extracellular 
compartments.  Using  a  simplified  model  of  fast  exchange 
where  there  may  be  intermediate-to-slow  exchange  can 
lead  to  significant  underestimation  of  vascular  volume. 
Experimental  approaches  to  minimize  these  errors  are 
based  on  observations  that  the  initial  slope  of  the  relax¬ 
ation  curve  is  independent  of  the  exchange  rate  (74.75). 

Large-molecular- weight  contrast  agents  potentially 
may  be  used  to  measure  tumor  blood  flow  by  detecting  the 
first  pass  of  the  agent  through  tumor  vasculature,  similar 
to  the  method  described  previously  (50,51),  although  this 
approach  may  not  be  feasible  when  the  heart  beat  is  very 


rapid,  as  for  rodents.  Studies  exploring  the  feasibility  of 
such  an  approach  are  yet  to  be  performed. 

In  our  laboratory,  we  used  the  dynamic  distribution  of 
the  intravascular  paramagnetic  contrast  agent  albumin- 
Gd-DTPA  to  quantitate  vascular  volume  and  permeability 
and  related  the  images  obtained  to  the  distribution  of 
VEGF  and  cell  morphology  in  histologic  sections  ob¬ 
tained  from  the  imaged  slice.  Vascular  volume  and  per¬ 
meability  were  characterized  for  two  different  tumor  mod¬ 
els  (MDA-MB-23 1  human  breast  cancer  line  in  severe 
combined  immune  deficient  (SCID)  mice  and  RIF-1  in 
C3H  mice)  with  the  primary  aim  of  determining  the  ability 
of  the  technique  to  detect  differences  in  these  parameters 
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Tt  Maps  in  msec 

3000  msec 


pre-aIbGdDTPA 


post-albGdDTPA 
(18  min) 


post-albGdDTPA 
(34  min) 


post-aIbGdDTPA 
(50  min) 


FIG.  3.  TI  maps  and  maps  of  vascular  volume  and  permeability  generated  from  the  images  shown  in  Rg.  1.  Vascular  volume  and  permeability  product  ^ 
surface  area  (PS)  maps  were  generated  from  the  ratio  of  A(  l/Tl )  values  in  the  images  to  that  of  blood  using  an  IDL  (Interactive  Data  Language,  Research  ^ 
Systems,  Inc.)  program.  The  slope  of  A(1/TI)  ratios  versus  time  in  each  pixel  was  used  to  compute  PS,  and  the  intercept  of  the  line  at  zero  lime  was  ^ 
used  to  compute  vascular  volume.  The  gray  scale  is  labeled  with  the  minimum  and  maximum  values  of  vascular  volume  and  permeability. 


for  these  two  models.  Such  a  technique  then  may  be  ap¬ 
plied  to  determine  the  correlation  between  vascular  vol¬ 
ume  and  permeability  for  tumors  derived  from  clones  with 
different  metastatic  capacity.  We  detected  significant  dif¬ 
ferences  in  vascular  volume  generated  by  the  MDA-MB- 
231  breast  cancer  line  compared  to  the  RIF-1  cell  line.  We 
also  detected  significant  differences  in  vascular  perme¬ 
ability  for  MDA-MB-231  tumors  growing  in  the  mam¬ 
mary  fat  pad  compared  to  MDA-MB-231  tumors  growing 
in  the  flank.  These  latter  results  confirm  previous  obser¬ 
vations  made  for  brain  tumor  cells  (76)  that  the  tumor 
microenvironment  can  dictate  the  permeability  of  the  vas¬ 
culature. 

Tl  maps  were  obtained  for  a  7CX)-iJLm  sagittal  slice 
through  the  center  of  the  tumor  using  a  progressive  satu¬ 
ration  recovery  technique  for  saturation  recovery  time  in¬ 
tervals  of  100  ms,  500  ms,  and  1  sec  followed  by  spin- 
echo  imaging  for  each  phase  encoding  step.  A  0.2-ml 
bolus  of  a  solution  of  60  mg/ml  albumin-Gd-DTPA  made 
up  in  saline  was  injected  through  the  tail  vein  (dose  of  500 
mg/kg).  Tl  maps  were  obtained  before  and  for  three  to 


five  time  points  after  the  intravenous  injection.  At  the  end 
of  the  imaging  studies,  the  animal  was  sacrificed,  0.5  ml  of 
blood  was  withdrawn  from  the  inferior  vena  cava  to  mea¬ 
sure  blood  Tl,  and  the  tumors  were  excised  and  fixed  in 
10%  buffered  formalin  for  sectioning  and  staining.  Pre- 
contrast  blood  Tl  values  were  obtained  from  a  separate 
tumor-bearing  group  of  four  SCID  and  three  C3H  mice. 
Tl  values  of  blood  collected  in  a  tube  were  measured 
using  inversion  recovery  with  10  recovery  time  points. 
Five  5-p.m-thick,  paraffin-embedded  histologic  sections 
were  obtained  from  the  imaged  slice  and  immunohisto- 
chemically  stained  for  VEGF  expression  (Santa  Cruz  Bio¬ 
technology,  Santa  Cruz,  CA,  U.S.A.)  (16).  Adjacent  sec¬ 
tions  were  stained  with  hematoxylin  and  eosin  for 
morphologic  information. 

A  representative  dataset  demonstrating  the  raw  images, 
obtained  in  this  case  for  an  MDA-MB-231  flank  tumor,  is 
shown  in  Fig.  2.  The  corresponding  Tl  maps  as  well  as 
maps  of  vascular  volume  and  permeability  generated  from 
these  images  are  shown  in  Fig.  3.  Both  vascular  volume 
and  permeability  were  spatially  heterogeneous,  and  this 
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H  &  E  VEGF 


Vascular  volume  ^rmeability  product  surface  area 

no.  4.  Adjacent  5-fun-thick.  histologic  sccUons  obtained  from  the  imaged  slice  in  Figs.  /an?/.  The  arrow  in  both  sectiorw  indicates  an  area  containing 
necrotic  cells.  The  necrotic  area  is  shown  at  a  higher  magnification  within  the  inset  of  each  section.  One  section  is  stained  with  hematoxylin  and  eosin 
(section  marked  H  &  E  in  figure);  the  other  section  is  stained  for  vascular  endothelial  growth  factor  (VEGF)  expression  using  immunoperoxidasc  staining 
(section  marked  VEGF  in  figure).  Regions  with  an  increased  brown  stain  indicate  higher  levels  of  VEGF  expression.  Maps  of  vascular  volume  and 
permeability  are  included  for  comparison. 


was  observed  for  all  the  tumors  examined.  Most  tumors 
showed  a  higher  vascular  volume  around  the  periphery, 
which  was  consistent  with  the  presence  of  higher  vessel 
density  in  the  histologic  sections  and  consistent  with  pre¬ 
vious  observations  that  peripheral  vascularization  is  one  of 
two  basic  patterns  of  tumor  vascularization,  the  other  be¬ 
ing  central  vascularization  (8).  Regions  with  the  highest 
vascular  volume  usually  did  not  exhibit  high  permeability. 
Two  adjacent  tumq^ sections  obtained  from  the  imaged 
slice  (shown  in  Figs,  and  stained  with  hematoxy- 

lin  and  eosin  and  immunostained  for  VEGF  expression  are 
.  displayed  in  Fig.  4  with  the  corresponding  maps  of  vas- 
c— ^  cular  volume  and  permeability  included  for  comparison. 


Regions  of  low  vascular  volume  in  the  MRI  maps  fre¬ 
quently  contained  necrotic  foci  (see  magnified  inset  in  Fig. 
4)  in  the  histologically  stained  sections:  the  regions  cor¬ 
responding  to  low  vascular  volume  usually  were  more 
permeable,  as  can  be  observed  by  comparing  the  regions 
marked  by  an  arrow  in  Fig.  4.  A  possible  explanation, 
consistent  with  the  necrotic  foci  observed  in  these  areas,  is 
that  areas  with  low  vascular  volume  may  contain  the  most 
leaky  and  nonfunctional  vasculature.  VEGF  staining  was 
most  intense  around  necrotic  areas  (see  arrow  in  Fig.  4).  A 
coarse  spatial  agreement  was  observed  between  areas  with 
high  VEGF  expression  in  the  histologic  section  and  areas 
of  high  permeability  observed  in  the  MR  images.  Vascular 
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volumes  of  MDA-MB-23I  mammary  fat  pad  tumors  (46.6 
±  9  tJLi/g,  n  =  6)  were  significandy  higher  (p  <  0.05)  than 
RIF-1  tumors  (18.7  ±  4  p,l/g,  n  =  4);  values  represent 
mean  ±  1  SEM.  MDA-MB-231  tumor  vasculature  was 
significantly  (p  <  0.016)  more  permeable  in  the  flank  (1.64 
±  0.26  |il/g  •  min,  n  =  4)  compared  to  the  mammary  fat 
pad  (0.85  ±  0.13  p.l/g  •  min,  n  =  6).  Despite  the  lower 
vascular  volume,  the  volume  doubling  time  of  the  RIF-1 
tumor  model  was  3  days  compared  to  14  days  for  the 
MDA-MB-231  model,  yet  RIF-1  tumors  at  the  volumes 
studied  here  do  not  exhibit  necrosis.  These  data  demon¬ 
strate  the  complex  balance  that  exists  among  vascular  vol¬ 
ume,  tumor  growth,  and  necrosis.  Cell  cycle  time,  cell  loss 
factor,  the  rate  of  dead  cell  clearance,  the  ability  of  cells  to 
generate  an  angiogenic  response,  the  endothelial  cell  pro¬ 
liferation  rate,  the  growth  fraction,  and  the  energy  and 
oxygen  requirements  of  cells  all  a  role  in  this  balance 

CONCLUSION 

The  undoubted  advantage  of  contrast-enhanced  MRI 
lies  in  its  ability  to  obtain  vascular  parameters  noninva- 
sively  with  in-plane  spatial  resolutions  on  the  order  of  200 
pim  or  less.  Thus,  it  is  possible  to  determine  the  relation¬ 
ship  between  vascular  volume  and  vascular  permeability 
within  a  single  pixel  and  to  study  the  dynamics  of  this 
relationship  as  well  as  other  vascular  properties  with  tu¬ 
mor  growth  and  therapy.  Using  MRI  to  study  transgenic 
tumors  with  specific  genetic  alterations,  in  combination 
with  immunostaining  of  specific  molecules  not  only  on 
endothelial  cells  but  targeted  at  various  stages  of  the  an¬ 
giogenic  pathway,  will  provide  an  even  greater  versatility 
for  understanding  tumor  angiogenesis  and  vascularization 
especially  when  combined  with  three-dimensional  inter¬ 
active  visualization  software  (78).  With  three-dimensional 
MRI  maps,  or  multislice  imaging,  it  should  be  possible  to 
measure  regions  with  the  highest  vascular  volume  and 
permeability,  as  well  as  the  total  tumor  areas  or  volumes 
associated  with  these  high  vascular  volume  and  perme¬ 
ability  to  understand  the  relationship  between  vasculariza¬ 
tion  and  metastasis.  Clinically,  some  of  the  major  foresee¬ 
able  applications  of  contrast-enhanced  MRI  are  to  increase 
the  specificity  for  detecting  malignant  lesions  and  for  the 
prognosis  of  malignant  lesions,  to  detect  the  efficiency  of 
antiangiogenic  treatments  as  they  become  clinically  avail¬ 
able,  and  equally  importantly  to  use  the  uptake  and  dis¬ 
tribution  of  contrast  agents  to  predict  the  efficiency  of 
delivery  of  therapeutic  agents  of  similar  size  to  the  tumor 
(79). 
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INTRODU  CTION ;  Magnetic  Resonance  Imaging  and 
spectroscopy  techniques  are  continually  evolving  to  address  critical 
issues  in  tumor  physiology  such  as  parameters  of  tumor 
vascularization.  Consequently  there  is  a  necessity  to  develop 
visualization  software  able  to  analyze  the  MR  information  and  relate 
it  to  data  obtained  from  the  same  tumor  with  other  technologies 
such  as  immunostaining,  autoradiography  and/or  PET.  Such  a 
multi-modal  approach  will  allow  a  high  level  of  versatility  in 
probing  the  disease  and  may  result  in  revealing  of  new  information 
hidden  in  the  original  data.  Here  we  present  results  obtained  with 
our  novel  visualization  software  which  we  are  developing  for  these 
purposes. 

METHODS:  3D  maps  of  blood  volume  and  vascular 

permeability  obtained  from  a  tumor  as  described  in  (1)  were 
compared  with  three  dimensional  histological  maps  reconstructed 
form  tissue  slices  stained  with  hematoxylin  eosin  (H&E)  and  for 
VEGF  expression.  Briefly,  multislice  Ti  images  were  obtained 
from  prostate  cancer  MLL  xenografts  in  SCID  mouse  with  Turbo- 
FLASH  saturation  recovery  technique  (64x64  resolution,  8  slices, 
100ms,  500ms,  and  Is  recovery  delays)  with  an  experimental  time 
of  7  min  per  multislice  map.  An  Mq  image  used  for  pixel  by  pixel 
reconstruction  of  Tj  maps  was  acquired  once  at  the  beginning  of  the 
experiment  with  a  recovery  delay  of  7  s.  Images  were  obtained 
before  and  for  32  min  after  i.y.  administration  of  0.2ml  of  60 
mg/ml  albumin-GdDTPA  in  saline  (dose  of  500mg/kg),  starting  10 
min  after  the  injection.  At  the  end  of  the  imaging  studies,  the 
animal  was  sacrificed,  0.5  ml  of  blood  was  withdrawn  from  the 
inferior  vena  cava,  and  tumors  mariced  for  referencing  to  the  MRI 
images  were  excised,  and  fixed  in  10%  formalin  for  sectioning  and 
staining.  Permeability  product  surface  area  (PS)  and  vascular 
volume  [W]  maps  were  generated  on  a  pixel  by  pixel  basis  from 
the  slope  and  intersect  at  zero  time  of  the  linear  plots  of  ratios  of 
A(1/T1)  values  in  the  images  to  that  of  blood  vs.  time  (1).  Ten 
pairs  of  adjacent  5p,m  histological  sections  obtained  at  500  \im 
intervals  through  the  tumor  were  stained  with  H&E  and 
immunostained  for  VEGF  expression  (2).  Sections  were  digitized 
with  a  Sanyo  CCD  camera  attached  to  an  optical  microscope.  3-D 
reconstruction  and  coregistration  of  MRI  and  histological  sections 
was  performed  using  die  Clinical  Microscope  Visualization  software 
that  we  are  developing  as  a  research  tool  for  our  investigation.  The 
software  consists  of  two  sets  of  routines.  The  first  module  builds 
volume  images  from  2-D  multislice  histology  images  and 
optionally  provides  stitching  of  high  resolution  image  'tiles' 
acquired  from  the  histological  sections  with  the  microscope  at  high 
magnification.  The  tiles  stitched  along  the  overlapping  borders 
produce  a  large  'virtuar  field  of  view  that  represents  the  full  cross 
section  which  could  not  be  acquired  otherwise  in  a  single  actual  field 
of  view  at  the  selected  high  resolution.  This  virtual  image  can  be 
panned  around,  zoomed  in,  to  look  at  selected  regions  at  the  full 
resolution  in  which  they  were  acquired.  The  second  module  of  the 
software  performs  visualization  of  the  volumetric  model  in  suitable 
3-D  perspective.  The  volume  images  can  be  interactively  lendaed 
using  the  volumetric  visualization  software  The  transfer  functions 
that  control  the  voxel  transparency  and  intensity  characteristics  of 
various  structures  of  interest  can  be  interactively  adjusted  to 
delineate  them  from  the  surrounding  structures.  Interactive 
registration  of  the  MR  data  with  histology  data  is  done  by  rendering 
the  two  volume  data  sets  together  using  individual  transfer  function 
to  control  the  overlapping  of  the  voxels  from  each  volume  so  that 
the  two  volunies  can  be  visualized  together.  Each  volume  can  be 
independenUy  manipulated  in  the  3b  space  to  achieve  better 
registration.  The  software  system  is  developed  around  Silicon 
Graphics  Inc.,  Workstation  systems,  taking  full  advantage  of  the 
hardware  accelerated  graphics  capabilities  such  as  2-  and  3-  D 
Textures  to  provide  interactive  rendering  results. 


RESULTS:  Figures  1  and  2  demonstrate  3-dimensional 
coregistration  of  3D  images  of  blood  volume  and  vascular 
permeability  with  3D  histological  maps  stained  with  H&E  and  for 
VEGF  expression  respectively. 


Figure  1:  3D  superimposed  maps  of  vascular  volume  and  H&E 
staining  obtained  from  a  tumor. 


Figure  2:  3D  map  of  vascular  permeability  and  VEGF  expression 
obtained  from  the  tumor  shown  in  Figure  1. 

DISCUSSION;  The  3D  multi-modal  volume  rendering  provides 
unique  information  about  the  relationship  between  vascular  volume, 
permeability  and  histology  maps  The  rendering  algorithm  preserves 
the  color  encoding  in  the  histology  maps  and  uses  gray  scales  for 
MR  maps,  which  allows  us  to  simultaneously  display  both  images 
and  navigate  through  them  for  direct  comparison  of  functional 
parameters  of  tumor  vasculature  with  cellular  characteristics  of  the 
tumor.  We  are  currently  investigating  other  approaches  to  combine 
multiple  volumes,  such  as  using  the  separate  color  channels  for  MR 
maps  and  histology  intensity  data  with  suitable  color  space 
transformations  to  identify  regions  of  spatial  correspondence. 
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INTRODUCTION:  Our  ongoing  studies  are  directed  towards 
understanding  the  role  of  vascularization  and  permeability  in  metastatic 
dissemination.  The  nm  (nonmetastatic)23  gene  is  related  to 
suppression  of  metastasis  [1].  Two  highly  homologous  and 
evolutionarily  conserved  nm23  genes,  nm23‘Hl  and  nm23-H2  have 
been  identified  in  humans.  These  genes  encode  17kDa  proteins 
identified  as  nucleoside  diphosphate  kinases  (NDPK)  A  and  B.  We 
recently  detected  significant  differences  in  the  metabolic  and 
physiological  characteristics  of  metastatic  MDA-MB-435-Vp  tumors 
containing  the  vector  only,  and  less  metastatic  MDA-MB-435-lp 
Iransgcne  tumors  containing  the  nm23"Hl  gene  (2].  Here  we  have  used 
MRI  to  quantitate  vascular  volume  and  permeability  of  metastatic  MDA- 
MB-435-Vpand  less  metastatic  MDA*MB-435-ip  tumors,  to  further 
understand  the  impact  of  the  loss  of  the  nm23  gene  on  tumor  vascular 
characteristics.  We  have  related  multi-slice  MRI  data  of  vascular  volume 
and  permeability  to  maps  of  VEGF  expression  and  cellular  morphology. 
METHODS:  Multi-slice  experiments,  relaxation  rates  (Ti‘^)  were 
obtained  by  a  saturation  recovery  method  combined  with  fast  Tj 
SNAPSHOT-FLASH  imaging  (flip  angle  of  5°,  echo  time  of  2  msec). 
Images  of  8  slices  (slice  thickness  of  1mm)  acquired  with  an  in-plane 
spatial  resolution  of  250p.m  (64x64  matrix,  16mm  field  of  view, 
NS=16)  were  obtained  for  3  relaxation  delays  (100msec,  500mscc,  and 
Isec)  for  each  of  the  slices.  Thus,  64x64x8  Tj  maps  were  acquired 
within  7  minutes.  An  Mq  map  with  a  recovery  delay  of  7s  was  acquired 
once  at  the  beginning  of  the  experiment.  Images  were  obtained  before 
i.v.  administration  of  0.2ml  of  60  mg/ml  albumin-GdDTPA  in  saline 
(dose  of  500mg/kg)  and  repeated  every  8  minutes,  starting  10  minutes 
after  the  injection,  up  to  32  minutes.  Relaxation  maps  were 
reconstructed  from  data  sets  for  three  different  relaxation  times  and  the 
Mq  data  set  on  a  pixel  by  pixel  basis.  At  the  end  of  the  imaging  studies, 
the  animal  was  sacrificed,  0.5  ml  of  blood  was  withdrawn  from  the 
inferior  vena  cava,  and  tumors  were  marked  for  referencing  to  the  MRI 
images,  excised,  and  fixed  in  10%  formalin  for  sectioning  and  staining. 
Vascular  volume  and  permeability  product  surface  area  (PS)  maps  were 
generated  from  the  ratio  of  A(1/T1)  values  in  the  images  to  that  of 
blood.  The  slope  of  A(1/TI)  ratios  versus  time  in  each  pixel  was  used  to 
compute  (PS)  while  the  intercept  of  the  line  at  zero  time  was  used  to 
compute  vascular  volume  (1).  Thus,  vascular  volumes  were  corrected  for 
permeability  of  the  vessels.  Adjacent  54m  thick  histological  sections 
obtained  at  500  \lm  intervals  through  the  tumor  were  stained  with 
hematoxylin  eosin  or  for  VEGF  distribution  using  a  rabbit  polyclonal 
anti- VEGF  antibody  (Santa  Cruz  Biotechnology,  Santa  Cruz,  CA)  . 
Sections  were  digitized  with  a  Sanyo  CCD  camera  attached  to  an  optical 
microscope.  3-D  reconstruction  of  both  MRI  and  histological  sections 
was  performed  using  the  Clinical  Microscope  Visualization  software  we 
are  developing  as  a  research  tool  for  our  studies.  The  software  system  is 
developed  around  our  Silicon  Graphics  Inc.  Octane  Workstation,  taking 
full  advantage  of  hardware  accelerated  graphics  capabilities  such  as  2- 
and  3-  D  Textures  to  provide  interactive  rendering  results. 

RESULTS:  7^“ 


Typical  multi-slice  data  of  vascular  volume  and  permeability  together 
with  histological  sections  stained  with  haematoxylin  and  eosin  and 
VEGF  expression,  obtained  from  an  MDA-MB-435-lp  tumor  (tumor 
volume  =  370  mm^)  are  shown  in  Figure  1.  The  corresponding  3D 
•  reconstructed  maps  of  (a)  vascular  volume  (b)  vascular  permeability  (c) 
Haematoxylin  Eosin  stained  map  and  (d)  VEGF  expression  are  shown  in 
Figure  2. 


The  table  below  shows  Mean  ±  S.E.M.  values  of  vascular  volume  and 
permeability  obtained  for  a  group  of  tumors;  n  =  number  of  animals. 


Tumor  Type 

Tumor  Volume 
in  mm3 

Vascular 
Volume  (jil/gm) 

Vascular 

Permeability 

(pl/gm-min) 

VP(n=6) 

349 

6.7  +  1.2 

0.39  +  0.15 

IP  (n=6) 

460 

9.04  +  1.8 

1.03  +  0.3* 

*  p  <  0.05,  unpaired  t-test 


No  significant  differences  in  vascular  volumes  averaged  over  the  whole 
tumor  were  detected  between  the  metastatic  MDA-MB-435-Vp  and  the 
less  metastatic  MDA-MB-435-lp  transgene  tumor.  Thus  average 
vascular  volume  does  not  appear  to  be  affected  by  the  loss  of  the  nm23 
gene.  Significant  differences,  however,  were  detected  for  vascular 
permeability;  the  less  metastatic  MDA-MB-435-lp  tumors  were 
significantly  more  permeable  compared  to  the  Vp  tumors.  These  values, 
however,  include  regions  of  necrosis.  We  are  currently  using  our 
interactive  software  to  analyze  vascular  volume  and  permeability  of 
viable  regions  only.  A  visual  comparison  of  VEGF  expression  in 
histological  sections  obtained  from  a  tumor  in  the  Vp  group  (0.3  nl/gm- 
min)  and  the  most  permeable  tumor  in  the  Ip  group  (1.6  |il/gm-min) 
showed  no  significant  differences  in  VEGF  staining  intensity.  This 
suggests  that  the  presence  of  nm23  may  alter  vascular  permeability 
independently  of  VEGF  expression.  Most  of  these  tumors  (Vp  and  Ip) 


exhibit  necrosis  (see  histological  sections  in 
Fig.  1  and  2)  and  therefore  must  contain  stressed 
environments.  Since  the  stress  response  of 
normal  tissues  is  also  characterized  by  increased 
permeability  (3),  our  working  hypothesis  is  that 
loss  of  the  nm23  gene  in  metastatic  breast 
cancer  may  result  in  a  reduction  of  stress  induced 
increase  of  vascular  permeability. 
REFERENCES:  1.  Steeg,  P.  et  al.,  JNCI,  80, 
1988;  2.  Bhujwalla  et  al.,  ISMRM  Proceedings, 
1998;  3.  Maeda,  H.  et  al,  Biochem.,  67,  1998. 
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The  Anti-Inflammatory  Agent  Indomethacin  Reduces  the  Malignant  Phospholipid  Phenotype 
of  Metastatic  Human  Breast  Cancer  Cells 

Z.  M.  Bhujwalla,  K,  Natarajan,  E.  0.  Aboagye,  N.  Mori,  D.  Artemov,  U.  Pilatus  and  V.  P.  Chacko. 

Oncology  Section  -  Division  of  MR  Research,  Dept,  of  Radiology,  The  Johns  Hopkins  University  School  of 

Medicine,  Baltimore,  MD  21205. 


INTRODUCTION:  We  recentiy  demonstrated  that  human  mammary 
epithelial  cells  exhibit  alterations  in  choline  phospholipids  with 
progression  to  the  malignant  phenotype  (1).  The  most  striking 
changes  detected  in  proton  spectra  of  cell  extracts  were  an  Increase  of 
glycerophosphocholine  (GPC)  relative  to  phosphocholine  (PC)  for  the 
non-malignant  lines  and  a  progressive  increase  in  PC  levels  with 
increasingly  malignant  phenotype.  Further  evidence  of  the  role  of  the 
choline  phospholipid  compounds  in  breast  cancer  progression  also 
comes  from  our  studies  of  the  metastatic  human  breast  cancer  line  MDA- 
MB-453  following  transfection  with  the  non  metastatic  gene  nm23; 
lowly  metastatic  tumors  derived  from  nm23  transfected  cells  exhibited  a 
significantly  higher  phosphodiesters  (PDE)  relative  to 
phosphomonoesters  (PME)  compared  to  control  metastatic  tumors 
derived  from  cells  transfected  with  vector  only  (2).  Several  studies  have 
shown  that  the  nonsteroidal  anti-inflammatory  agent  indomethacin 
(INDO)  can  inhibit  the  invasive  and  metastatic  behaviour  of  human 
breast  cancer  cells.  Our  question  here  was:  Does  treatment  of  highly 
malignant  breast  cancer  cells  with  INDO  alter  the  PDE/PME  ratios 
consistent  with  predictions  from  our  observations  of  the 
'invasive/metastatic  phospholipid  phenotype’  ?  We  studied  intact 
human  MDA-MB-435  (wild  type)  breast  cancer  cells  immobilized  on 
microcarrier  beads  in  our  NMR  cell  perfusion  system  and  also  obtained 
proton  spectra  from  extracts  of  cells  growing  in  tissue  culture  flasks. 
METHODS;  For  the  ^^P  NMR  experiments  of  intact  cells,  2  -  3  ml  of 
lO'^  MDA-MB-435  breast  cancer  cells/ml  of  microcarrier  beads  were 
transferred  into  a  10  mm  screw  cap  NMR  tube.  The  cells  were 
continuously  perfused  at  a  rate  of  i.5ml/min  with  growth  medium 
consisting  of  RPMI  supplemented  with  10%  FCS  and  containing  15mM 
HEPES.  The  medium  was  equilibrialed  with  95%  O2  and  5%  CO2.  ^*P 
NMR  spectra  (45®  flip  angle;  1  sec  repetition  time;  na=4000)  were 
obtained  at  80  min  intervals  from  these  intact  isolated  perfused  cells, 
prior  to,  and  following  treatment  with  100  pM  INDO.  Parallel  control 
experiments  of  untreated  cells  were  performed  over  the  same  time  period. 
For  the  1 H  NMR  extract  studies,  cells  were  cultured  in  supplemented 
DMEM  without  HEPES  but  with  bicarbonate  buffer.  Cells  were  fed  with 


Proton  spectra  obtained  from  cell  extracts  of  MDA-MB-435  cells  (Figure 
2)  also  showed  a  significant  increase  of  GPC  relative  to  PC  following 
incubation  with  280  pM  INDO  for  3  hrs.  There  were  no  significant 
differences  in  cell  viability  for  the  treated  and  control  flasks. 


medium  alone  or  medium  plus  280  pM  INDO  3  hours  prior  to  extraction  DISCUSSION:  Treatment  with  INDO  altered  the  phospholipid  profile 
and  used  at  70-80%  confluency.  Perchloric  acid  (PCA)  extracts  were  of  highly  malignant  breast  cancer  cells  towards  a  less  malignant 
obtained  from  approximately  I08  cells.  Cells  were  trypsinized  and  the  phenotype.  Recently,  we  have  also  detected  similar  changes  following 
action  of  trypsin  was  blocked  by  cold  serum  containing  growth  medium.  INDO  for  MDA-MB-231  cells.  Since  INDO  is  a  cyclooxygenase  (COX) 
Cells  were  then  washed  twice  with  normal  saline  and  extracted  with  ice  inhibitor,  these  results  suggest  that  the  altered  phospholipid  profile  of 
cold  8%  (v/v)  PCA.  Supernatants  were  neutralized  with  3M  K2CO3/  IM  highly  malignant  cells  may  partly  be  related  to  alterations  of  the  COX- 
KOH,  lyophilized,  and  resuspended  in  D2O.  Trimethylsilyl  propionate  l/COX-2  pathway.  COX-1  and  COX -2  arc  rate-limiting  enzymes  which 
(TSP)  was  used  as  internal  standard.  Fully-relaxed  water-suppressed  convert  arachidonic  acid  to  prostanoids;  COX-2  is  an  inducible  enzyme 
proton  spectra  of  the  extracts  were  acquired  on  a  1 1.7T  Bruker  5(X)MHz  and  its  constitutive  expression  has  been  shown  to  enhance  the  invasive 
NMR  spectrometer  with  a  5  mm  probe.  Spectral  acquisition  parameters  and  metastatic  behaviour  of  colorectal  cancer  cells  (3).  Our  results 
were  SO®  flip  angle,  SW  =  6000  Hz.  4.7  sec  repetition  time.  32  K  block  Provide  direct  evidence  for  the  first  time  that  the  phospholipid  peaks  in 
size,  and  number  of  scans  =  512.  spectra  of  breast  cancer  cells  can  be  modulated  by  a 

RESULTS;  Treatment  with  INDO  resulted  in  a  significant  increase  of  'y^looxygenase  inhibitor  and  indicate  the  utility  of  NMR  in 
GPC  and  glycerophosphoethanolamine  (OPE);  control  spectra  obtained  “"^eptanding  the  role  of  anti-inflammatory  agents  in  cancer  treatment, 
over  a  similar  time  period  showed  no  changes  in  the  PDE  region.  The 

increase  in  GPC  (shown  at  9  hrs  in  Figure  1)  was  evident  within  80  j  REFERENCES:  1.  Aboagye,  E.  and 

minutes  following  addition  of  100  liM  INDO  in  the  medium.  .  Bhujwalla,  Z.  M.  Cancer  Research,  1999, 

2.  Bhujwalla  et  al.,  ISMRM  Proceedings, 
1998;  3.  Tsujii,  M.  et  al,  PNAS,  1997. 
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Characterization  of  Lactate  Levels  in  Normal  and  Malignant  Human  Mammary 

Epithelial  Cells 

Eric  0.  Aboagye  and  Zaver  M.  Bhujwalla 

Oncology  Section  -  Division  of  MR  Research,  Department  of  Radiology, 

The  Johns  Hopkins  University  School  of  Medicine,  Baltimore,  MD  21205,  USA 


INTRODUCTION ;  Interest  in  lactate  production  by  cancer  cells  and 
the  role  of  lactate  in  tumors  has  existed  since  studies  by  Warburg  in  the 
1930's  [1].  Recently  it  was  shown  that  patients  with  cervical  cancers 
containing  high  lactate  had  a  high  risk  of  metastases  [2].  Two  possible 
explanations  for  the  findings  by  Schwickert  et  al  [2]  are  (i)  A  high 
glycolytic  rate  is  part  of  the  malignant  invasive  phenotype  or  (ii)  High 
lactate  is  indicative  of  hypoxia  and  substrate  deprivation  which 
stimulates  or  promotes  invasion  and  metastasis.  Supporting  the  latter 
possibility  are  observations  that  glucose  starvation,  acidosis  and 
hypoxia  enhance  the  metastatic  potential  of  murine  tumor  cells  [3,  4]. 
In  this  study  we  have  explored  the  first  possibility  for  human  breast 
cells  by  determining  lactate  levels  in  a  series  of  human  mammary 
epithelial  cells  (HMECs)  including  finite  lifespan  cells  (normal, 
senescent),  spontaneously  or  benzo(a)pyrene  immortalized  cells  and 
tumor-derived  cells.  Our  results  indicate  that  malignant  breast  cancer 
cells  produce  a  significantly  higher  amount  of  lactate  compared  to 
normal  or  immortalized  HMECs. 

METHODS;  The  epithelial  cell  lines  used  here  have  been  previously 
described  [5],  They  include  184  and  48  finite  lifespan  cells,  184A1  and 
184B5  which  were  derived  from  184  by  benzo(a)pyrene 
immortalization,  spontaneously  immortalized  MCF12A  cells,  185B5- 
erbB2  cells  which  were  derived  from  184B5  by  oncogene  transformation 
with  erbB2,  and  the  breast,  cancer  cell  lines  SKBR3,  MCF7j  MDA-MB- 
231  and  MDA-MB-435.  Cells  were  cultured  in  MCDB  170  or 
DMEM/F12  with  or  without  supplements  as  previously  described  [5]. 
For  NMR  measurements,  cells  were  fed  24  hours  prior  to  extraction  and 
used  at  70-80%  confluency.  Perchloric  acid  (PCA)  extracts  were 
obtained  from  10^  to  10^  cells.  Cells  were  trypsinized  and  the  action  of 
trypsin  was  blocked  with  cold  serum  containing  growth  medium.  Cells 
were  then  washed  twice  with  normal  saline  and  extracted  with  ice  cold 
8%  (v/v)  PCA.  Supernatants  were  neutralized  (with  3M  K2CO3/  IM 
KOH),  lyophilized  and  resuspended  in  D2O.  Trimethylsilyl  propionate 
(TSP)  was  used  as  an  internal  standard.  Fully-relaxed  water-suppressed 
proton  spectra  of  the  extracts  were  acquired  on  a  1L7T  Bruker  500MHz 
NMR  spectrometer  with  a  5  mm  probe.  Spectral  acquisition  parameters 
were  30^  flip  angle,  SW  =  6000  Hz,  4.7  sec  repetition  time,  32  K  block 
size,  and  number  of  scans  =  512. 

RESULTS:  The  lactate  doublet  was  detected  by  IH  NMR  spectroscopy 
of  PCA  extracts  of  cells  at  a  resonance  frequency  of  1.2ppm  downfield  of 
TSP.  High  resolution  IH  NMR  spectra  with  comparable  signal  to  noise, 
of  extracts  obtained  from  immortalized  MCF12A  cells  and  malignant 
MDA-MB-231  cells  are  shown  in  Figure  1.  A  large  difference  in  lactate 
levels  is  apparent  between  these  cell  lines.  A  summary  of  lactate  levels 
calculated  for  each  cell  line  and  normalized  to  cell  number  and  cell  size  is 
shown  in  Fig  2.  Unlike  the  changes  reported  for  choline  metabolites 
[5],  no  gradual  increase  in  lactate  levels  was  observed  in  our  model 
system.  All  tumor-derived  cells,  however,  showed  significantly  higher 
levels  of  lactate  compared  to  the  normal,  immortalized  and  oncogene 
transformed  cells  (p  ^  0.05).  There  was  no  correlation  between  lactate 
levels  and  cell  doubling  time  measured  by  the  MTT  assay.  Lactate  levels 
in  tumors  derived  from  MCF7,  MDA-MB-231  and  MDA-MB-435 
;  correlated  with  the  levels  measured  in  cells.  In  addition  to  lactate  levels, 
;  we  also  analyzed  total  creatine  (at  3.04  ppm)  for  these  lines  (shown  in 
i  Fig.  3).  With  the  exception  of  the  most  mali^ant  cell  line,  MDA-MB- 
’  435,  aJl  cell  lines  showed  low  levels  of  total  creatine. 
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.  Figure  2:  Lactate  levels  for  cell  lines  representing  various  stages  of 
malignant  progression  of  HMECs.  Bars  represent  ±  S.E.M.  There  was  a 
statistically  significant  difference  in  lactate  (p  ^  0.05)  for  finite 

lifespan  and  immortalized  versus  tumor-derived  cells  (Ifmol/fim^  = 
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;  Figure  3;  Creatine  levels  in  a  panel  of  cell  lines  representing  various 
stages  of  malignant  progression  of  HMECs.  Bars  represent  ±  S.E.M. 
DISCUSSION:  Our  results  demonstrate  that,  unlike  choline  levels 
which  gradually  increased  with  malignant  progression  in  this  model 
system  [5],  lactate  levels  remained  low  in  all  normal,  immortalized  and 
oncogene,  transformed  cells.  Consistent  with  published  data  for  human 
breast  tumors  versus  uninvolved  breast  tissues  [6],  lactate  levels  were 
significantly  higher  for  the  malignant  cell  lines.  Following  malignant 
transformation,  however,  there  appeared  to  be  no  further  dependence  of 
Uactate  levels  with  increasing  invasive/metastatic  behaviour  (e.g.  MCF- 
7  vs  MDA-MB-231).  Increased  lactate  levels  may  be  due  to  increased 
expression  of  glycolytic  enzymes  such  as  lactate  dehydrogenase  and/or 
of  glucose  transporters.  These  data  suggest  that  the  increase  in  lactate 
;  levels,  characteristic  of  most  tumors,  occurs  late  in  carcinogenesis.  Our 
data  also  demonstrate  that,  with  the  exception  of  the  MDA-MB-435 
line,  mammary  epithelial  cells  do  not  show  an  increase  in 
creatine/phosphocreatine  levels  with  malignant  transformation.  Given 
the  importance  of  lactate  levels  in  the  diagnosis  of  cancer  and  prediction 
of  tumor  response,  the  results  obtained  here  support  further  studies  to 
understand  the  mechanism  of  regulation  of  lactate  in  breast  epithelial 
cells. 
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